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Abstract 

These two accompanying papers are concerned with two mode entanglement 
for systems of identical massive bosons and the relationship to spin squeezing 
and other quantum correlation effects. Entanglement is a key quantum feature 
of composite systems in which the probabilities for joint measurements on the 
composite sub-systems are no longer determined from measurement probabili¬ 
ties on the separate sub-systems. There are many aspects of entanglement that 
can be studied. This two-part review focuses on the meaning of entanglement, 
the quantum paradoxes associated with entangled states, and the important 
tests that allow an experimentalist to determine whether a quantum state - in 
particular, one for massive bosons is entangled. An overall outcome of the re¬ 
view is to distinguish criteria (and hence experiments) for entanglement that 
fully utilise the symmetrisation principle and the super-selection rules that can 
be applied to bosonic massive particles. 

In the hrst paper (I), the background is given for the meaning of entan¬ 
glement in the context of systems of identical particles. For such systems, the 
requirement is that the relevant quantum density operators must satisfy the 
symmetrisation principle and that global and local super-selection rules pro¬ 
hibit states in which there are coherences between differing particle numbers. 
The justification for these requirements is fully discussed. In the second quan¬ 
tisation approach that is used, both the system and the sub-systems are modes 
(or sets of modes) rather than particles, particles being associated with different 
occupancies of the modes. The definition of entangled states is based on hrst 
dehning the non-entangled states - after specifying which modes constitute the 
sub-systems. This work mainly focuses on the two mode entanglement for mas¬ 
sive bosons, but is put in the context of tests of local hidden variable theories, 
where one may not be able to make the above restrictions. The review provides 
the detailed arguments necessary for the conclusions of a recent paper, where 
the question of how to rigorously demonstrate the entanglement of a two-mode 
Bose-Einstein condensate (BEC) has been examined. 

In the accompanying review paper (II), we consider spin squeezing and other 
tests for entanglement that have been proposed for two-mode bosonic systems. 
We apply the approach of review (I) to determine which tests, and which mod¬ 
ifications of the tests, are useful for detecting entanglement in massive bosonic 
(BEC), as opposed to photonic, systems. Several new inequalities are derived, 
a theory for the required two-mode interferometry is presented, and key exper¬ 
iments to date are analysed. 


PACS Numbers 03.65 Ud, 03.67 Bg, 03.67 Mn, 03.75 Gg 


2 



Contents 


1. Introduction 

1.1 Entanglement: Definition and Historical Context 

1.2 Measures and Tests for Entanglement 

1.3 Particle versus Mode Entanglement 

1.4 Symmetrization and Super-Selection Rules 

1.5 Entanglement Tests and Experiments - Paper II 

1.6 Outline of Papers I and II 

2. Entanglement - General Features 

2.1 Quantum States 

2.2 Entangled and Non-Entangled States 

2.2.1 General Considerations 

2.2.2 Local Systems and Operations 

2.2.3 Constraints on Sub-System Density Operators 

2.2.4 Classical Entanglement 

2.3 Separate, Joint Measurements. Reduced Density Operator 

2.3.1 Joint Measurements on Sub-Systems 

2.3.2 Single Sub-System Measurements. Reduced Density Operator 

2.3.3 Mean Value and Variance 

2 . 3.4 Conditional Probabilities 

2.3.5 Conditional Mean and Variance 

2.4 Non-Entangled States 

2 . 4.1 Non-Entangled States - Joint Measurements on Sub-Systems 

2 . 4.2 Non-Entangled States - Single Sub-System Measurements 

2 . 4.3 Non-Entangled States - Conditional Probability 

2 . 4.4 Non-Entangled States - Mean Values and Correlations 

2.5 Local Hidden Variable Theory (LHV) 

2.5.1 LHV - Mean Values and Correlations 

2.5.2 LHV - GHZ State 

2.6 Paradoxes 

2.6.1 EPR Paradox 

2.6.2 Schrodinger Cat Paradox 

2.7 Bell Inequalities 

2.7.1 LHV Result 

2.7.2 Non-Entangled State Result 

2.7.3 Bell Ineguality Violation and Entanglement 

2.8 Non-Local Correlations 

2.8.1 LHV Result 

2.8.2 Non-Entangled State Result 

2.8.3 Correlation Violation and Entanglement 

3. Identical Particles and Entanglement 

3.1 Symmetrisation Principle 


3 



3.1.1 Sub-Systems - Particles or Modes ? 

3.1.2 Multi-Mode Sub-Systems 

3.2 Super-Selection Rules - General 

3.2.1 Global Particle Number SSR 

3.2.2 Examples of Global SSR Compliant States 

3.2.3 SSR and Conservation Laws 

3 . 2.4 Global SSR Compliant States and Quantum Correlation Func¬ 
tions 

3.2.5 Testing the SSR 

3.2.6 SSR Justification - No Suitable Phase Reference 

3.2.7 SSR Justification - Physics Considerations 

3.2.8 SSR Justification - Galilean Frames ? 

3.2.9 SSR and Photons 

3.3 Reference Frames and SSR Violation 

3.3.1 Linking SSR and Reference Frames 

3.3.2 Can Coherent Superpositions of an Atom and a Molecule 

Occur ? 

3.3.3 Detection of SSR Violating States 

3.4 Super-Selection Rules - Separable States 

3 . 4.1 Local Particle Number SSR 

3 . 4.2 Criteria for Local and Global SSR in Separable States 

3 . 4.3 States that are Global but not Local SSR Compliant 

3 . 4.4 Particle Entanglement Measure 

3 . 4.5 General Form of Non-Entangled States 

3.5 Bipartite Systems 

3.5.1 Two Single Modes - Coherence Terms 

3.5.2 Two Pairs of Modes - Coherence Terms 

4. Discussion and Summary of Key Results 

References 

Acknowledgements 

Appendices 

A. Projective Measurements and Conditional Proba¬ 
bilities 

A.r Projective Measurements 
A. 2 Conditional Probabilities 

A. 3 Conditional Mean and Variance 

B. Inequalities 

B. r Integral Inequality 

B. 2 Sum Inequality 

C. EPR Spin Paradox 

C. l Local Spin Operators 


4 



C. 2 Conditioned Variances 

D. Extracting Entanglement due to Symmetrisation 

D. l Three Particle Case - Bosons 

D.2 Two Particle Case - Bosons 

D. 3 Two Particle Case - Fermions 

E. Reference Erames and Super-Selection Rules 

E. l Two Observers with Different Reference Frames 

E.2 Symmetry Groups 

E.3 Relationships - Situation A 

E.4 Relationships - Situation B 

E.5 Dynamical and Measurement Considerations 

E.6 Nature of Reference Frames 

E.7 Relational Description of Phase References 

E.8 Irreducible Matrix Representations and SSR 

E. 9 Non-Entangled States 

E. Super-Selection Rule Violations ? 

F. l Preparation of a Coherent Superposition of an Atom and 
a Molecule ? 

F.1.1 Hamiltonian 

F.l.2 Initial State 

F.1.3 Implicated Reference Frame 

F. 1.4 Process - Alice and Charlie Descriptions 

F.l.5 Interference Effects Without SSR Violation 

F.1.6 Conclusion 

F.2 Detection of Coherent Superposition of a Vacuum and 
One-Boson State ? 


5 



1 Introduction 


Since the paradoxes of Einstein-Podolski-Rosen (EPR) [T] , Schrodinger [2] and 
Bell [3], entanglement has been recognised as a key feature that distinguishes 
quantum physics from classical physics. Entangled quantum states underlie the 
EPR and Bell paradoxes, which reveal the conflict between quantum mechanics 
and local realism, and the famous Schrodinger cat paradox, where a cat is appar¬ 
ently prepared in a state simultaneously both dead and alive. Entanglement not 
only provides a way to rigorously test quantum principles, but is the basis for 
the many quantum information tasks like quantum cryptography. Despite the 
fundamental interest, there have been only a few experimental tests of entan¬ 
glement for systems of massive particles. Yet, the substantial recent progress in 
cooling atomic systems, in particular to form Bose-Einstein condensates, makes 
such entanglement tests much more feasible. 

In this review (I), we explain the meaning of entanglement, and examine 
how to verify entanglement, for systems of identical boson particles. This leads 
us to focus on symmetrisation and superselection rules, and to consider their 
implication for entanglement criteria when applied to massive bosonic particles. 
This paper provides the theoretical background for a recent paper [4] and a 
subsequent review (II) that analyses the suitability of specific criteria, new 
and old, to detect entanglement in bosonic systems, and applies the criteria to 
interpret experimental findings. 

As well as reviewing the topic and presenting some new results in review 
II, these two articles are intended as comprehensive papers for post docs and 
postgraduates who are changing field or starting work in a new one and need to 
gain a thorough understanding of the present state of knowledge. With this aim 
in mind we have not followed the conventional approach in review articles of 
merely quoting formulae and giving references, but instead have presented full 
proofs of the key results. To really understand a field, we believe it is necessary 
to work through the derivations. However, in order to shorten the main body 
of the articles, we have included many of the details in Appendices. 

1.1 Entanglement: Definitions and Historical Context 

Entanglement arises in the context of composite quantum systems composed of 
distinct components or sub-systems and is distinct from other features of quan¬ 
tum physics such as quantization for measured values of physical quantities, 
prohahilistic outcomes for such measurements, uncertainty principles involving 
pairs of physical quantities and so on. Such sub-systems are usually associ¬ 
ated with sub-sets of the physical quantities applying to the overall system, 
and in general more than one choice of sub-systems can be made. The formal¬ 
ism of quantum theory treats pure states for systems made up of two or more 
distinct sub-systems via tensor products of sub-system states, and since these 
product states exist in a Hilbert space, it follows that linear combinations of 
such products could also represent possible pure quantum states for the system. 
Such quantum superpositions which cannot be expressed as a single product of 
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sub-system states are known as entangled (or non-separable) states. 

The concept of entanglement can then be extended to mixed states, where 
quantum states for the system and the sub-systems are specified by density op¬ 
erators. The detailed definition of entangled states is set out in Section [51 This 
definition is based on first carefully defining the non-entangled (or separable) 
states. The set of non-entangled states must allow all possible quantum states 
for the given sub-system, but in addition these states must be preparable via pro¬ 
cesses involving separate operations on each sub-system after which correlated 
sub-system quantum states are combined in accordance with classical probabil¬ 
ities. Thus, although the sub-system states retain their quantum natures the 
combination resulting in the overall system state is formed classically rather 
than quantum mechanically. This overall process then involves local operations 
on the distinct sub-systems and classical communication [LOCC) to prepare a 
general non-entangled state. The entangled states are then just the quantum 
states which are not non-entangled states. The general idea that in all com¬ 
posite systems the non-entangled states all involve LOCC preparation processes 
was first suggested by Werner [^. The notion of quantum states, the nature 
of the systems and sub-systems involved and the specific features required in 
the definition of non-entangled states when identical particles are involved is 
discussed in detail in Section [3] Entangled states underlie a number of effects 
that cannot be interpreted in terms of classical physics, including spin squeezing, 
non-local measurement correlations - such as for the Einstein-Podolski-Rosen 
(EPR) paradox and violations of Bell Inequalities ([1], [2], [3], [6], [7]). The 
quantum theory of measurement i, i, m, [H] invokes entangled states of 
the system and measuring apparatus as key concepts in the theory. More re¬ 
cently, entangled states have been recognised as a resource that can be used in 
various quantum technologies for applications such as teleportation, quantum 
cryptography, quantum computing and so on. Recent expositions on the effects 
of entanglement and its role in quantum information science include 

US], m, [la, [Sl¬ 
it would be pointless to characterise states as entangled unless such states 
have some important properties. The key requirement is that entangled states 
exhibit a novel quantum feature that is only found in composite systems. As will 
be seen in SubSection l2.3l separable states are such that the joint probability for 
measurements of all physical quantities associated with the sub-systems can be 
found from separate measurement probabilities obtained from the sub-system 
density operators and the overall classical probability for creating particular 
products of sub-system states. In general, entangled states do not exhibit this 
feature of separable probabilities, and it is this key non-separability feature 
that led Schrodinger to call these states ’’entangled”. Where the sub-systems 
are spatially separated, one can define spacelike separated local measurement 
events on each. This was historically the reason why the sub-systems and their 
measurements are often referred to as local. The EPR paper [T] suggested the 
possibility that although the predictions of quantum theory were correct, the 
theory was incomplete and there was an underlying reality in the form of classical 
hidden variables. Averaging over the unknown values of the hidden variables 
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would be required to produce the same measurement probabilities as quantum 
theory. Local hidden variable theories (LHV) are discussed in SubSection 12.51 
and it will be seen that the joint probabilities for measurements of sub-system 
physical quantities are of the same form as for separable states. As will be seen 
in Subsections 12.61 and O EPR or Bell inequality violations do not occur for 
states described either by LHV theories or as quantum separable states. Hence 
there is a direct link between EPR and Bell violations and both the failure of 
LHV theories and the presence of entanglement. The fact that certain entangled 
states do not exhibit the feature of separable probabilities shown in classical 
LHV theories highlights entanglement being a non-classical feature found in 
composite systems. 

Note that although an EPR or Bell inequality violation requires the quantum 
state to be entangled, there are examples of mixed entangled states that do 
not violate a Bell inequality. For pure states of qubits Gisin HZ] showed that 
entangled states always violated Bell inequality, but for mixed states Werner 
and others m, m. m have shown there are entangled states ( Werner states 
ca) for which a hidden variable theory can be constructed that gives the same 
joint probability function for measurement outcomes as quantum theory. Of 
course the quantum measurement outcomes must be given before the hidden 
variable model is constructed - there is no known way to determine the LHV 
theory distribution functions independently. These specific entangled states will 
therefore satisfy Bell inequalities. 

The mixed entangled states considered by Werner for which a hidden 
variable theory could be constructed were of a special form. Two distinguishable 
sub-systems each with d basis states \ur) were considered, for which the states 
could be transformed by unitary operators U, and the combined density operator 
'p was required to be invariant under all unitary transformations of the form 
U so that p = {U ®Uyp{U^ ®W). Werner considered the following 
unitary operators: (a) U-r such that U-r \ur) = — \ur), U-r |us) = -|- |ms) for 
s ^ r (b) P{r —>■ pr), which permute basis states P{r —>• fir) \ur) = (c) 

Urotin-jm), which transform basis states |u„), \um) into linear combinations of 
each other Urot (^, rn) | — bJnn \ '^n) T Unm |^m), bJrot (^, rn) I ^m) — bjYnn 

Umm \um)- As a consequcncc of these invariances Werner showed that the 
density operator could be expressed in terms of a single parameter $ in the 
form p = [d^ — d)~^ |(d— $)1 -I- (d$ — l)ld|, where 1 is the unit operator 

and V is the flip operator. These have matrix elements ( 1) = dm dsm 

\ / rs, nm 

and (v) = Srm dsn ■ From this form of the density operator Werner [5] 

V / rs, nm 

showed that the probability function for joint measurement outcomes on the two 
sub-systems could be expressed in the same form as applied in hidden variable 
theory. So although the mixed entangled state Werner considered were of a 
restricted type, the work demonstrated that entanglement did not preclude all 
hidden variable theory interpretations of the joint measurements. The fact that 
some entangled states do not violate a Bell inequality is another consequence of 
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Werner’s result. 

The issue of how best to treat the quantum aspects of correlations in mea¬ 
surement outcomes in composite quantum systems is still an active area of 
research and is beyond the scope of these two papers. Quantum entanglement 
is clearly relevant to the discussion, but concepts such as quantum discord |21j , 
and EPR steering [231, being used to describe quantum 

correlations. The link between these concepts is discussed in [26]. In these re¬ 
cent discussions of quantum correlation, it turns out that some separable states 
are regarded as exhibiting quantum correlations. 

It is now generally recognised that entanglement is a relative concept 
pS] . P5]i. [T3] . [30]) [SI] and not only depends on the quantum state under 
discussion but also on which sub-systems are being considered as entangled (or 
non-entangled). A quantum state may be entangled for one choice of the sub¬ 
systems but may be non-entangled if another choice of sub-systems is made. 
A simple example often cited is that for the hydrogen atom [20], a system 
made up of two distinguishable particles, a proton and an electron. Here the 
energy eigenstates are non-entangled if the sub-systems refer to the centre of 
mass of the entire atom and the relative position of the electron and the pro¬ 
ton, but which would be entangled if the sub-systems were the positions of 
the electron and proton. It could be argued that the centre of mass and the 
relative position are not really independent sub-systems - one always accompa¬ 
nies the other - but as unrelated centre of mass and relative position quantum 
states can be prepared, they can be regarded as distinct sub-systems. The in¬ 
dividual positions of the electron and the proton are also distinct sub-systems, 
and the ground state of the hydrogen atom is indeed entangled - the electron 
position is tightly correlated with the proton position. Another example in¬ 
volves two different choices of single particle states in a two mode Bose-Einstein 
condensate (BEC) - a system with a large number of identical particles. The 
issue of defining sub-systems will be dealt with below, but taking the origi¬ 
nal two sub-systems to be bosonic modes (or single particle states) denoted 
1^^) and a well known N boson entangled state of these two modes 

A and B (with mode annihilation operators d and h) is the binomial state 
given by |$) = ((cos0 exp(zx/2) a-I-sin0 exp(—ix/2) |0) (see [32] 

and Paper II, Section 3.7) which is a quantum superposition of Fock states 
^{N — k)l |0) with k = 0, ... ,N. Introducing new modes 
via c = {cos 9 exp(ix/2) a -1- sin0 exp(—ix/2) b) and d = {— sin0 exp(ix/2) a -b 
cosO exp(—*x/2) b) we see that we can also write |$) = (c^)^/\/M |0), so that 
the same quantum state is a separable state if the sub-systems are chosen to 
be the new modes C and D. Another example is the ground state of the sin¬ 
gle mode non-interacting BEC trapped in a harmonic oscillator (HO) potential. 
This is a separable state, with all bosons in the lowest energy mode if the sub¬ 
systems are chosen as the HO modes. However, if single particle position states 
spatially localised in two different regions are chosen as two sub-systems, then 
the same ground state for the identical particle system is spatially entangled, as 
pointed out by Goold et al [55] . 
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1.2 Measures and Tests for Entanglement 

Various measures of entanglement have been defined for certain types of quan¬ 
tum state - see m. M, ISO], m, EH. m. Ea, for details of these, and are 
aimed at quantifying entanglement to determine which states are more entangled 
than others. This is important since entanglement is considered as a resource 
needed in various quantum technologies. Calculations based on such measures 
of entanglement confirm that for some choices of sub-systems the quantum state 
is entangled, for others it is non-entangled. For two mode pure states the en¬ 
tanglement entropy - being the difference between the entropy for the pure state 
(zero) and that associated with the reduced density operator for either of the 
two sub-systems - is a useful entanglement measure. As entropy and information 
changes are directly linked m. |14j . this measure is of importance to quantum 
information science. Measurements of entanglement based on Renyi entropy 
and purity are discussed in ES], m and |38] . Another entanglement measure 
is particle entanglement, defined by Wiseman et al [33] j [10] j E3] identical 
particle systems and based on projecting the quantum state onto states with 
definite particle numbers. One of the problems with entanglement measures is 
that there is often no simple way to measure the quantities required. 

In the case of bipartite entanglement in qubit systems [H] and [33] obtained 
a sufficient condition for a quantum state to be entangled {PPT condition) (see 
|3T] . [30] for details). Suppose the density operator p is changed into p by 
mapping the matrix elements associated with one of the sub-systems into their 
transpose. Then provided the new operator p^ is a valid density operator (with 
real, non-negative eigenvalues that add to unity), the original density operator 
represents a separable state. Thus, if it is shown that some of the eigenvalues of 
p are negative, then the state p is entangled. However, it is often not practical 
to use this as an entanglement test for systems with large numbers of basis 
states, as it requires being able to measure all the density matrix elements. It 
was also later realised [43] that in general, the PPT condition was not a necessary 
condition for entanglement, apart from cases of 2 x 2 and 2x3 subsystems - 
that is, showing that p has only positive eigenvalues will not guarantee that p 
is separable, as counter-examples for 2 x 4 and 3x3 subsystems showed. 

Although not directly related to the various quantitative measures of en¬ 
tanglement, the results for certain measurements can play the role of being 
signatures or witnesses or tests of entanglement EQ], [3T], [34]. These are in the 
form of inequalities for variances and mean values for certain physical quanti¬ 
ties, which are dependent on the inequalities applying for non-entangled quan¬ 
tum states. If such inequalities are violated then it can be concluded that the 
state is entangled for the relevant sub-systems. In the case of entanglement 
witnesses, Jdie idea is to find a hermitian operator W such that for separable 
states Tr{Wp) > 0, so that if Tr{Wf)) < 0 the state must be entangled. Here 
we note that the density operator occurs linearly when evaluating the quan¬ 
tities involved. Some of the correlation tests discussed in paper H are cases 
involving entanglement witnesses. However, in more general tests for entangle¬ 
ment the density operator appears non-linearly. For example, a spin squeezing 
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test for entanglement may require showing that the variance for a spin oper¬ 
ator is less than a multiple of the magnitude of the mean value of another 
spin operator - thus for example < | ^5*2^ |/2. This could be written 

as Tr{{Sl ± S^/2)'p) — {Tr{Sx)'p)Y < Oj which is of a more general form than 
for an entanglement witness. Non-linear tests are discussed in Ref. |31j . One 
of the advantages of entanglement tests is that the quantities involved can be 
measured. It cannot be emphasised enough that these tests provide sufficiency 
conditions for establishing that a state is entangled. So if the test is satisfied 
we can conclude that the state is not separable. The failure of a test does not 
mean that the state is not entangled - sufficiency does not imply necessity. The 
violation of a Bell inequality is an example of such a signature of entanglement, 
and the demonstration of spin squeezing is regarded as another. However, the 
absence of spin squeezing (for example) does not guarantee non-entanglement, 
as the case of the NOON state in SubSection 3.6 of the accompanying paper 
II shows. A significant number of such inequalities have now been proposed 
and such signatures of entanglement are the primary focus of the accompanying 
paper, which is aimed at identifying which of these inequalities really do iden¬ 
tify entangled states, especially in the context of two mode systems of identical 
bosons. 

At present there is no clear linkage between quantitative measures of entan¬ 
glement (such as entanglement entropy) and the quantities used in conjunction 
with the various entanglement tests (such as the relative spin fluctuation in 
spin squeezing experiments). Results for experiments demonstrating such non- 
classical effects cannot yet be used to say much more than the state is entangled, 
whereas ideally these experiments should determine how entangled the state is. 
Again we emphasise that neither the entanglement tests nor the entanglement 
measures are being used to define entanglement. Entanglement is defined first 
as being the quantum states that are non-separable, the tests for and measures 
of entanglement are consequential on this definition. 

1.3 Particle versus Mode Entanglement 

These two papers deal with identical particles - bosons or fermions. In the sec¬ 
ond quantisation approach used here the system is regarded as a set of quantum 
fields, each of which may be considered as a collection of single particle states 
or modes. We now take into account the situation where systems of identi¬ 
cal particles are involved. This requires us to give special consideration to the 
requirement that quantum states in such cases must conform to the symmetri- 
sation principle. What sub-systems are possible must take into account that 
entanglement requires the specification of sub-systems that are distinguishable 
from each other and on which measurements can be made. In addition, the 
sub-systems must be able to exist as separate systems which can be prepared in 
quantum states for that sub-system alone. These key requirements that the sub¬ 
systems must be distinguishable, susceptible to measurements and can exist in 
separate quantum states are necessary for the concept of entanglement to make 
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physical sense, and have important consequences for the choice of sub-systems 
when identical particles are involved. These three key logical requirement for 
sub-systems rule out considering labelled identical particles as sub-systems and 
lead to the conclusion that sub-systems must be modes. Thus both the system 
and sub-systems will be specified via the modes that are involved, so here the 
sub-systems in terms of which non-entangled (and hence entangled) states are 
defined are modes or sets of modes, not particles m. Es, [29]), m, m, m- 
In this approach, particles are associated with the occupancies of the various 
modes, so that situations with differing numbers of particles will be treated as 
differing quantum states of the same system, not as different systems - as in 
the first quantisation approach. Note that the choice of modes is not unique - 
original sets of orthogonal one particle states (modes) may be replaced by other 
orthogonal sets. An example is given in Section 2 of accompanying paper II. 
Modes can often be categorised as localised modes, where the corresponding 
single particle wavefunction is confined to a restricted spatial region, or may 
be categorised as delocalised modes, where the opposite applies. Single parti¬ 
cle harmonic oscillator states are an example of localised modes, momentum 
states are an example of delocalised modes. This distinction is significant when 
phenomena such as EPR violations and teleportation are considered. 

However, even if the system consists entirely of distinguishable particles we 
can still regard the sub-systems as collections of modes. Each distinguishable 
particle is still associated with a set of single particle states or modes (momen¬ 
tum eigenstates, harmonic oscillator eigenstates, etc.) that can be occupied. 
More general states associated with a single particle may be quantum superposi¬ 
tion states of those with a single particle occupancy of the modes. If the overall 
system consists of a number of distinguishable particles each of which is con¬ 
sidered as a sub-system, then each such sub-system can equally be regarded as 
the set of modes associated with the particular distinguishable particle. Overall 
system states involving just one particle of each type would be simultaneous 
eigenstates of the number operators for each of the distinguishable particles, 
with an eigenvalue of unity corresponding to there being only one particle of 
each type. The second quantisation approach can still be used, but is somewhat 
superfluous when the modes for each particle are only occupied once. 

Although multi-mode systems are also considered, in this paper we mainly 
focus on two mode systems of identical bosonic atoms, where the atoms at most 
occupy only two single particle states or modes. For bosonic atoms this situa¬ 
tion applies in two mode interferometry, where if a single hyperfine component 
is involved the modes concerned may be two distinct spatial modes, such as 
in a double well magnetic or optical trap, or if two hyperfine components are 
involved in a single well trap each component has its own spatial mode. Large 
numbers of bosons may be involved since there is no restriction on the number 
of bosons that can occupy a bosonic mode. For fermionic atoms each hyper¬ 
fine component again has its own spatial mode. However, if large numbers 
of fermionic atoms are involved then as the Pauli exclusion principle only al¬ 
lows each mode to accommodate one fermion, it follows that a large number of 
modes must considered and two mode systems would be restricted to at most 
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two fermions. Consideration of multi-mode entanglement for large numbers of 
fermions is outside the scope of the present paper (see [4^ for a treatment of 
this), and unless otherwise indicated the focus will be on bosonic modes. The 
paper focuses on identical bosonic atoms - whether the paper also applies to 
photons is less clear and will be discussed below. 

1.4 Symmetrization and Super-Selection Rules 

The work presented here begins with the fundamental issue of how an entangled 
state should be defined in the context of systems involving identical particles. 
To reiterate - in the commonly used mathematical approach for defining entan¬ 
gled states, this requires first defining a general non-entangled state, all other 
states therefore being entangled. We adhere to the original definition of Werner 
in which the separable states are those that can be prepared by local op¬ 
erations and classical communication [LOCC). This approach is adopted by 
other authors, see for example [47], [24], [48]. However, in other papers - see for 
example |49] . [50] so-called separable non-local states are introduced in which 
LOCC is not required (see SubSection ]3.4.3l for an example). By contrast (and 
consistent with Werner’s approach), in the present paper it is contended that 
the density operators both for the quantum states of the overall system and 
those for the non-entangled (local) sub-systems in the context of non-entangled 
states must be compatible with certain principles and rules that have been 
found to be both necessary and sufficient for understanding physical experi¬ 
ments in non-relativistic many-body systems. In some other work (discussed 
below) this has not been the case. A key feature required of all quantum states 
for systems involving identical particles., entangled or not is that they satisfy 
the symmetrization principle ED- This places restrictions both on the form of 
the overall density operator and also on what can be validly considered to be 
a sub-system. In particular this rules out individual identical particles being 
treated as sub-systems, as is done in some papers (see below). If the system 
consists entirely of distinguishable particles then the symmetrization principle 
is not relevant. In addition, super-selection rules (SSR) [52] only allow den¬ 
sity operators which have zero coherences between states with differing total 
numbers of particles to represent valid quantum states, and this will be taken 
into account for all quantum states of the overall system, entangled or not. 
This is referred to as the global particle number super-selection rule In non- 
entangled or separable states the density operator is a sum over products of 
sub-system density operators, each product being weighted by its probability 
of ocurring (see below for details). For the non-entangled or separable quan¬ 
tum states, a so-called local particle number super-selection rule will also be 
applied to the density operators describing each of the sub-systems. These sub¬ 
system density operators must then have have zero coherences between states 
with differing numbers of sub-system particles. This additional restriction ex¬ 
cludes density operators as defining non-entangled states when the sub-system 
density operators do not conform to the local particle number super-selection 
rule. Consequently, density operators where the local particle number SSR does 
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not apply would be regarded as entangled states. This viewpoint is discussed 
in papers by Bartlett et al [47], [53] as one of several approaches for defining 
entangled states. However, other authors such as [H], [SD] state on the contrary 
that states where the sub-system density operators do not conform to the local 
particle number super-selection rule are still separable, others such as |54] . [55] 
do so by implication - the latter papers applied to atomic as well as photon 
modes. So in these two papers we are advocating a different definition to some 
other definitions of entanglement in identical particle systems, the consequence 
being that the set of entangled states is now much larger. This is a key idea in 
this paper - not only should super-selection rules on particle numbers be applied 
to the the overall quantum state, entangled or not, but it also should be applied 
to the density operators that describe states of the modal sub-systems involved 
in the general definition of non-entangled states. 

Note that for systems entirely consisting of N distinguishable particles the 
super-selection rules are still true, but are now superfluous. Each sub-system is 
the set of modes or one particle states of the specific distinguishable particle and 
the overall state is an N particle state in which the sub-systems only contain one 
particle. Consequently there are no sub-system or system coherences between 
states with differing particle number. 

The detailed reasons for adopting the viewpoint that the entanglement cri¬ 
teria be compliant with the requirement of the local particle number super- 
selecrtion rule (SSR) for the sub-system are set out below. As will be seen, 
the local particle number super-selection rule restriction firstly depends on the 
fundamental requirement that for all composite systems - whether identical 
particles are involved or not - non-entangled states are only those that can 
be prepared via processes that involve only LOCC. The requirement that the 
sub-system density operators in identical particle cases satisfy the local particle 
number SSR is consequential on the sub-system states being possible sub-system 
quantum states. As mentioned before, the general definition of non-entangled 
states based on LOCC preparation processes was first suggested by Werner [ 3 . 
Apart from the papers by Bartlett et al we are not aware that this 

LOCC/SSR based criteria for non-entangled states has been invoked previously 
for identical particle systems, indeed the opposite approach has been proposed 
[49] . [50] . However, the idea of considering whether sub-system states should 
satisfy the local particle number SSR has been presented in several papers - [49] , 
m. m, [53], [56], [57], [58] . mainly in the context of pure states for bosonic 
systems, though in these papers the focus is on issues other than the definition 
of entanglement - such as quantum communication protocols [49] , multicopy dis¬ 
tillation [ 43 , mechanical work and accessible entanglement [56] . m and Bell 
inequality violation [58] . The consequences for entanglement of applying this 
super-selection rule requirement to the sub-system density operators are quite 
significant, and in the accompanying paper H important new entanglement tests 
are determined. Not only can it immediately be established that spin squeezing 
requires entangled states, but though several of the other inequalities (accom¬ 
panying paper H) that have been used as signatures of entanglement are still 
valid, additional tests can be obtained which only apply to entangled states that 
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are defined to conform to the symmetrisation principle and the super-selection 
rules. 

It is worth emphasising that requiring the sub-system density operators sat¬ 
isfy the local particle number SSR means that there are less states than other¬ 
wise would be the case which are classed as non-entangled, and more states will 
be regarded as entangled. It is therefore not surprising that additional tests for 
entanglement will result. If further restrictions are placed on the sub-system 
density operator - such as requiring them to correspond to a fixed number of 
bosons again there will be more states regarded as entangled, and even more 
entanglement tests will apply. A particular example is given in SubSection 4.3 
of paper II, where the sub-systems are restricted to one boson states. 

The symmetrisation requirement for systems involving identical particles is 
well established since the work of Dirac. There are two types of justification 
for applying the super-selection rules for systems of identical particles (both 
massive and otherwise). The first approach for invoking the superselection rule 
to exclude quantum superposition states with differing numbers of identical 
particles is based on simple considerations and may be summarised as: 

1. No way is known for creating such SSR non-compliant states. 

2. No way is known for measuring the properties of such states. 

3. Coherence and interference effects can be understood in terms of SSR 
compliant states. 

The second approach is more sophisticated and involves linking the absence 
or presence of SSR to whether or not there is a suitable reference frame in 
terms of which the quantum state is described [59], [60] . m, m, [ 50 ], [ 62 ], 
[63] . [64] . [53], [56], [^, [34]. This approach will be described in SubSection l3.2l 
and Appendix [ini the key idea being that SSR are a consequence of considering 
the description of a quantum state by a real observer (Charlie) whose phase 
reference frame has an unknown phase difference from that of a hypothetical 
observer ((Alice), both studying the same system. Alice is assumed to possess 
a phase reference frame such that her description of the quantum state of the 
system violates the SSR. Charlie, on the other hand is an actual observer with 
no such phase reference frame. Thus, whilst Alice’s description of the system 
involves a quantum state may violate the SSR, the description of the same 
system by Charlie will involve a quantum state that is SSR compliant. In the 
main part of this paper the density operator p used to describe the various 
quantum states will be that of the external observer (Charlie). Note that if 
well-defined phase references do exist and the relationship between them is 
known, then the SSR can be challenged (see SubSection 13.31 and Appendix [TU]) . 
but this situation does not apply in the case of massive bosons (or fermions). 

It should be noted that both of these justifications for applying the SSR 
are dependent on what is practical in terms of measurements in non-relativistic 
quantum physics. Here the situation is clearer for systems of massive particles 
such as atoms than for massless particles such as photons. Applying SSR for 
photons is discussed in SubSubSection 13.2.91 

However, to allow for quantum states that as far as we know cannot be made 
or measured, and for which there are no known physical effects that require 
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their presence is an unnecessary feature to add to the non-relativistic quantum 
physics of many body systems or to quantum optics. Considerations based on 
the general principle of simplicity (Occam’s razor) would suggest not doing so 
until a clear physical justification for including them is found. The quantum 
state is intended to specify what is known about a quantum system and how it 
was prepared. It is used to determine the probabilities for possible measurements 
on the system. Clearly there is no point in including non-SSR compliant terms 
in the density operator for the quantum state. Such terms would neither allow 
for possible preparation processes, or contribute to measurement probabilities 
associated with physical efects. Furthermore, experiments can be carried out on 
each of the mode sub-systems considered as a separate system, and essentially 
the same reasons that justify applying the super-selection rule to the overall 
system also apply to the separate mode sub-systems in the context of defining 
non-entangled states. Hence, unless it can be justified to ignore the super¬ 
selection rule for the overall system it would be inconsistent not to apply it 
to the sub-system as well. As we will see, for separable states the requirement 
that the overall state is SSR compliant generally implies that the sub-system 
states are SSR compliant - though in some special cases this is not the case (see 
Subsection 4.3.3 of paper II). The onus is on those who wish to ignore the super¬ 
selection rule for the separate sub-systems to justify why it is being applied to 
the overall system. In addition, joint measurements on all the sub-systems can 
be carried out, and the interpretation of the measurement probabilities requires 
the density operators for the sub-system states to be physically based. The 
general application of super-selection rules has however been challenged (see 
Subsection on the basis that super-selection rules are not a fundamental 
requirement of quantum theory, but are restrictions that could be lifted if there 
is a suitable system that acts as a reference for the coherences involved. In 
Section [3] and in Appendix [TT] an analysis of these objections to the super¬ 
selection rule is presented, and in Appendix 1101 we see that the approach based 
on phase reference frames does indeed justify the application of the SSR both 
to the general quantum states for multi-mode systems of identical particles and 
to the sub-system states for non-entangled states of these systems. 

The sceptic who wishes to ignore the super-selection rules in the definition of 
entanglement - and consequently only consider as valid tests for entanglement 
where SSR compliance is not used in their derivation - needs to carry out a 
research program analogous to that which resulted in parity non-conservation 
becoming a basic feature of weak interaction theory. The successful incorpora¬ 
tion of parity non-conservation involved first proposing (on symmetry grounds) 
possible interactions in which parity was not conserved, second working out pos¬ 
sible experiments that could confirm parity non conservation and third carrying 
out key experiments that did confirm this. At this stage no such work in regard 
to SSR violation in non-relativistic many body physics has been carried out or 
is likely to be in the near future (except possibly for photons). As we will see 
in paper II, none of the experimental methods for entanglement tests that we 
examine can detect SSR non-compliance - none involve a suitable phase refer¬ 
ence. To ignore SSR in non-relativistic entanglement theory and experiment on 
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the grounds of scepticism would be analogous to ignoring parity conservation 
in quantum chemistry or atomic physics - areas which are well-understood in 
terms of parity being conserved (apart from the well-known parity violating ef¬ 
fects of external electric fields). When and if SSR violation in non-relativistic 
many body physics is found would then be the time to revise the definition 
of quantum entanglement. In these two papers we will utilise the definition 
of entanglement and derive tests based on SSR compliance, though of course 
recognising that there are also tests that do not require SSR compliance which 
are also valid for SSR compliant states. Although other definitions of entangle¬ 
ment will be considered for comparison, to avoid confusion the SSR compliant 
definition will be the one which we mean when we refer to entanglement. 

A further sound scientific argument can be presented in favour of studying 
SSR compliant entanglement tests (as is our aim in these papers (I) and (II)). 
This involves a consideration of what can be concluded from such tests by 
suppporters or sceptics of SSR. For example, one such test (see paper II) involves 
spin squeezing in two mode systems. If the state is separable and the sub-system 
states comply with local SSR then there is no spin squeezing. However, if 
experimental tests do demonstrate spin squeezing, then what can we conclude? 
The supporters of SSR compliance being required for the sub-systems would 
conclude that the state was not separable and hence entanglement is present 
between the subsystems. On the other hand, the sceptic who does not believe 
local SSR compliance is required would have no option but to conclude either 
that entanglement is present between the subsystems or (if they argue there is no 
entanglement) the state is separable but one or both of the quantum subsystems 
violates the SSR. The sceptic may favour the second conclusion, but that would 
then imply an actual experimental circumstance where superselection rules did 
not apply to the sub-system states. In that case, the issues raised in the last 
four paragraphs regarding laek of phase references or SSR violating preparation 
processes etc. must be addressed directly. Either way, the study of such SSR 
based experiments is clearly important. Put another way, suppose the sceptic 
were to derive a different test using the separability requirement alone, for 
which an experimental outcome shows that the two subsystems were indeed 
not entangled. This would seem to require a test for entanglement which is 
necessary as well as being sufficient - the latter alone being usually the case for 
entanglement tests. Such criteria and measurements are a challenge, but not 
impossible even though we have not met this challenge in these two papers. If 
the conclusion from the earlier SSR based experiment was either entanglement 
or separability with non SSR compliance, then if the result from the different 
test based only on separability ruled entanglement out, it follows that the system 
must be in a separable state in which the sub-system states violate the SSR. 
Conversely, the latter test may confirm the entanglement possibility found in 
the earlier test. Thus, in principle there could be a pair of experiments that 
give evidence of entanglement, or failure of the Super Selection Rule. For such 
investigations to be possible, the use of entanglement criteria that do invoke the 
local super-selection rules is also required. 
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1.5 Entanglement Tests and Experiments - Paper II 

The main focus of the accompanying paper II is to derive the SSR compliant 
criteria and to consider the experimental implementation. This leads to impor¬ 
tant links between spin squeezing and entanglement. The link with quantum 
correlation functions (as proposed in Refs. [51]. ISSj ! is also treated. Heisen¬ 
berg Uncertainty Principle inequalities involving spin operators |65j and the 
consequent property of spin squeezing have been well-known in quantum op¬ 
tics for many years. The importance of spin squeezing in quantum metrology 
is discussed in the paper by Kitagawa et al [66] for general spin systems. It 
was suggested in this paper that correlations between the individual spins was 
needed to produce spin squeezing, though no quantitative proof was presented 
and the more precise concept of entanglement was not mentioned. For the case 
of two mode systems the earliest paper linking spin squeezing to entanglement is 
that of Sorensen et al m, which considers a system of identical bosonic atoms, 
each of which can occupy one of two internal states. This paper states that spin 
squeezing requires the quantum state to be entangled, with a proof given in the 
Appendix. A consideration of how such spin squeezing may be generated via 
collisional interactions is also presented. The paper by Sorensen et al is often 
referred to as establishing the link between spin squeezing and entanglement - 
see for example Micheli et al (5^, Toth et al [53] j Hyllus et al m- However, 
the paper by Sorensen et al m is based on a definition of non-entangled states 
in which the sub-systems are the identical particles, and this is inconsistent 
with the symmetrization principle. However, the accompanying paper H estab¬ 
lishes the link between spin squeezing and entanglement based on a dehnition 
of entanglement consistent with the system and sub-system density operators 
representing quantum states. 

It is also important to consider which components of the spin operator vec¬ 
tor are squeezed, and this issue is also considered in the accompanying paper. 
In the context of the present second quantisation approach to identical particle 
systems the three spin operator components for two mode systems are expressed 
in terms of the annihilation, creation operators for the two chosen modes. Spin 
squeezing can be defined (see Section 2 in the accompanying paper) in terms 
of the variances of these spin operators, however the covariance matrix for the 
three spin operators will in general have off-diagonal elements, and spin squeez¬ 
ing is also defined in terms of rotated spin operators referred to as principal spin 
operators for which the covariance matrix is diagonal. The principal spin oper¬ 
ators are related to new mode annihilation, creation operators in the same form 
as for the original spin operators, where the new modes are two orthogonal lin¬ 
ear combinations of the originally chosen modes. In discussing the relationship 
between spin squeezing and entanglement, the modes which may be entangled 
are generally those associated with the definition of the spin operators. 

A further focus of the accompanying paper is on the relationship between 
entanglement and certain correlation properties of sub-system operators. Tests 
for entanglement based on such correlations have also been published - see for 
example [54], [55]. These tests were based on ignoring the super-selection rules. 
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so in the accompanying paper we present revised correlation tests for entangle¬ 
ment when the super-selection rules are definitely complied with. We also show 
the link between correlation tests and tests involving spin operators. 

The accompanying paper also deals with the important question of what 
measurement systems are suitable for making spin and correlation tests for en¬ 
tanglement. We first consider a simple two mode interferometer which involves 
coupling the two modes employing a resonant classical field pulse which is asso¬ 
ciated with a variable pulse area for its amplitude and has an adjustable phase. 
It is shown that measurement of the mean value and variance of the popula¬ 
tion difference between the two modes after the interferometer pulse enables 
measurements of the mean value and covariance matrix elements of the spin 
operators for the quantum state that existed before the pulse was applied. The 
mean values and variances of certain spin operators are relevant for correlation 
and spin squeezing entanglement tests. 

Paper II is focused on two mode bosonic systems. These are of particular 
interest because cold atomic gases cooled well below the Bose-Einstein conden¬ 
sation (BEC) transition temperature can be prepared where essentially only two 
modes are occupied ([ZI],I32]). This can be achieved for cases involving a single 
hyperfine components using a double well trap potential or for two hyperfine 
components using a single well. At higher temperatures more than two modes 
may be occupied, so multi-mode systems are also of importance and thus are 
considered in paper II. 

1.6 Outlines of Papers I and II 

The plan of the present paper is as follows. In Section!^ the key definitions of 
entangled states are covered, and a detailed discussion on why the symmetri- 
sation principle and the super-selection rule is invoked in discussed in Section 
[3l Challenges to the necessity of the super-selection rule are outlined, with ar- 
guements against such challenges dealt with in Appendices (TU] and [TT] Two key 
mathematical inequalities are derived in Appendix [T] Details for the spin EPR 
paradox are given in Appendix [8] The final Section |4] summarises and discusses 
the key features about entanglement treated in this paper. 

In the accompanying paper II, Section 2 sets out the definitions of spin 
squeezing and in the following Section 3 it is shown that spin squeezing is a 
signature of entanglement, both for the original spin operators with entangle¬ 
ment of the original modes and the principle spin operators with entanglement 
of the two new modes, and also for multi-mode cases. Details are in Appen¬ 
dices A and B. A number of other tests for entanglement proposed by other 
authors are considered in Sections 4, 5 and 6, with details of these treatments 
set out in Appendices B, C and D. In Section 7 it is shown that a simple two 
mode interferometer can be used to measure the mean values and covariance 
matrix for the spin operators involved in entanglement tests. The treatment is 
then generalised to situations involving measurements on multi-mode systems. 
Details are covered in Appendices G and H. Actual experiments aimed at 
detecting entanglement via spin squeezing tests are examined in Section 8. The 
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final Section 9 summarises and discusses the key results regarding entanglement 
tests. Appendices E, F and I provide details regarding certain important states 
whose features are discussed in the paper - the ’’separable but non-local ” states 
and the relative phase eigenstate. 
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2 


Entanglement - General Features 

2.1 Quantum States 

The standard Copenhagen quantum theory notions of physical systems that 
can exist in various states and have associated properties on which measure¬ 
ments can be made are presumed in this paper. The measuring system may 
be also treated via quantum theory, but there is always some component that 
behaves classically, so that quantum fluctuations in the quantity recorded by 
the observer are small. The term quantum state (or ’’physical quantum state” 
or just ’’state” for short) refers to a state that can either be prepared via a 
process consistent with the laws of quantum physics and on which measure¬ 
ments can be then performed and the probabilistic results predicted from this 
state (prediction), or a state whose existence can be inferred from later quantum 
measurements (retrodiction). We may also refer to such states as allowed quan¬ 
tum states, and our approach is intended to be physically based. In quantum 
theory, quantum states are represented mathematically by density operators for 
mixed states or state vectors for pure states. For identical particle systems these 
representations must satisfy symmetrisation and other basic requirements in ac¬ 
cordance with the laws of quantum theory. The probabilities of measurement 
outcomes and the probabilities associated with retrodiction can be interpreted 
as Bayesian probabilities m m, and the quantum state is observer depen¬ 
dent. The quantum state, the system it is associated with and the quantities 
that can be measured are considered here as entities that are viewed as be¬ 
ing both ontological and epistimological. Different observers may have different 
information about how the quantum state was prepared, hence the quantum 
state is in part epistimological, and would be described differently by different 
observers. Hence the observer is important, but as there is actually something 
out there to be studied, quantum states also have an ontological aspect. We 
will avoid the unqualified term ’’physical state” because this term is generally 
invoked in discussions about the pre-Copenhagen notion of reality and refers to 
some as yet unknown but more fundamental description of the system which 
underlies the quantum state m- Hidden variable theories attempt to describe 
this more fundamental physical state that is assumed to exist - attempts that 
so far have been unsuccessful if locality is also invoked (see below). In addition 
to those associated with physical quantum states, other density operators and 
state vectors may be introduced for mathematical convenience. For physical 
quantum states, the density operator is determined from either the preparation 
process or inferred from the measurement process - quantum tomography - and 
in general it is a statistical mixture of density operators for possible prepara¬ 
tion processes. Measurement itself constitutes a possible preparation process. 
Following preparation, further experimental processes may change the quantum 
state and dynamical equations give the time evolution of the density opera¬ 
tor between preparation and measurement, the simplest situation being where 
measurement takes place immediately after preparation. A full discussion of the 
predictive and retrodictive aspects of the density operator is given in papers by 


21 


Pegg et al m, m- Whilst there are often different mathematical forms for the 
density operator that lead to the same predictive results for subsequent mea¬ 
surements, the results of the measurements can also be used to retrodictively 
determine the preferred form of the density operator that is consistent with the 
available preparation and measurement operators. An example is given in m- 


2.2 Entangled and Non-Entangled States 

2.2.1 General Considerations 

Here the commonly applied physically-based approach to mathematically defin¬ 
ing entangled states will be described M- The definition involves vectors and 
density operators that represent states than can be prepared in real experi¬ 
ments, so the mathematical approach is to be physically based. The concept 
of quantum entanglement involves composite systems made up of component 
sub-systems each of which are distinguishable from the other sub-systems, and 
where each could constitute a stand-alone quantum system. This means the 
each sub-system will have its own set of physically realisable quantum states - 
mixed or pure - which could be prepared independently of the quantum states 
of the other sub-systems. As will be seen, the requirement that sub-systems 
be distinguishable and their states be physically preparable will have important 
consequences, especially in the context of identical particle systems. The formal 
definition of what is meant by an entangled state starts with the pure states, 
described via a vector in a Hilbert space. The formalism of quantum theory 
allows for pure states for composite systems made up of two or more distinct 
sub-systems via tensor products of sub-system states 

1$) = 1$^) (g) |$b) 0 |$C7) - (1) 


Such products are called non-entangled or separable states. However, since these 
product states exist in a Hilbert space, it follows that linear combinations of 
such products of the form 


a/S-y.. 



( 2 ) 


could also represent possible pure quantum states for the system. Such quantum 
superpositions which cannot be expressed as a single product of sub-system 
states are known as entangled (or non-separable) states. 

The concept of entanglement can be extended to mixed states, which are 
described via density operators in the Hilbert space. If A, B, ... are the sub- 
systems with p^, being density operators the sub-systems A, B, .then a 
general non-entangled or separable state is one where the overall density opera¬ 
tor p can be written as the weighted sum of tensor products of these sub-system 
density operators in the form [S] 


? = ^Pr'Pr®?'r®?%® ■■■ 

R 


(3) 
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with J2rPr — 1 Pr > 0 giving the probability that the specific product 
state Pfi = Pr® Pr® Pr® occurs. It is assumed that at least in princi¬ 
ple such separable states can be prepared [5]. This implies the possibility of 
turning off the interactions between the different sub-systems, a task that may 
be difficult in practice except for well-separated sub-systems. Entangled states 
(or non-separable states) are those that cannot be written in this form, so in 
this approach knowing what the term entangled state refers to is based on first 
knowing what the general form is for a non-entangled state. The density opera¬ 
tor p = 1$) ($1 for the pure state in ([2]) is not of the form ([31), as there are cross 


terms of the form 0 ( 






) ® ... involved. 


The concepts of separability and entanglement based on the Eqs. o and m 
for non-entangled states do not however just rest on the mathematical forms 
alone. Implicitly there is the assumption that separable quantum states de¬ 
scribed by the two expressions can actually be created in physical processes. 
The sub-systems involved must therefore be distinguishable quantum systems 
in their own right, and the sub-system states I'I'a) > I'f’s)) ••• or p^, pf, ... must 
also be possible quantum states for the sub-systems. We will return to these 
requirements later. The issue of the physical preparation of non-entangled (sep¬ 
arable) states starting from some uncorrelated fiducial state for the separate 
sub-systems was introduced by Werner , and discussed further by Bartlett et 
al (see im, Section IIB). This involves the ideas of loeal operations and elassieal 
communieation (LOCC) dealt with in the next SubSection. 

The key requirement is that entangled states exhibit a novel quantum feature 
that is only found in composite systems. Separable states are such that the joint 
probability for measurements of all physical quantities associated with the sub¬ 
systems can be found from separate measurement probabilities obtained from 
the sub-system density operators pf, etc and the overall classical proba¬ 
bility Pr (see SubSectioi l2.3p . This feature of separable probabilities is absent 
in certain entangled states, and because of this key non-separability feature 
Schrodinger called these states ’’entangled”. The separability feature for the 
joint probabilities is essentially a classical feature and applies in hidden variable 
theories (HVT) (see SubSectioi l2.5l) applied to quantum systems - as well as to 
quantum separable states. The fact that entangled states are quantum states 
that can exhibit the failure of this separability feature for classical LHV theories 
highlights entanglement being a non-classical feature for composite systems. 

An alternative operational approach to defining entangled states focuses on 
whether or not they exhibit certain non-classical features such as Bell Inequal¬ 
ity violation or whether they satisfy certain mathematical tests such as having 
a non-negative partial transpose |41| . |30) . and a utilititarian approach focuses 
or whether entangled states have technological applications such as in various 
quantum information protocols. As will be seen in SubSection l3.41 the particular 
definition of entangled states based on their non-creatability via LOCC essen¬ 
tially coincides with the approach used in the present paper. It has been realised 
for some time that different types of entangled states occur, for example states 
in which a Bell inequality is violated or states demonstrating an EPR paradox 
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[75]. Wiseman et al [73], |73j, [TS] and Reid et al [77], [T5], [T5], [75] discuss the 
concept of a heirarchy of entangled states, with states exhibiting Bell nonlocality 
being a subset of states for which there is EPR steering, which in turn is a sub¬ 
set of all the entangled states, the latter being defined as states whose density 
operators cannot be written as in Eq. (jS]) though without further consideration 
if additional properties are required for the sub-system density operators. The 
operational approach could lead into a quagmire of differing interpretations of 
entanglement dependng on which non-classical feature is highlighted, and the 
utilitarian approach implies that all entangled states have a technological use - 
which is by no means the case. For these reasons, the present physical approach 
based on the quantities involved representing allowed sub-system states is gen¬ 
erally favoured [14] . It is also compatible with later classifying entangled states 
in a heirarchy. 

2.2.2 Local Systems and Operations 

As pointed out by Vedral m, one reason for calling states such as in Eqs. (HD and 
m separable is associated with the idea of performing operations on the separate 
sub-systems that do not affect the other sub-systems. Such operations on such 
local systems are referred to as local operations and include unitary operations 
Ua, Ub, that change the states via —>■ UaPbU^ , —>■ UbPrU^ , etc as in 
a time evolution, and could include processes by which the states p(^, pf, are 
separately prepared from suitable initial states. 

We note that performing local operations on a separable state only produces 
another separable state, not an entangled state. Such local operations are obvi¬ 
ously faciltated in experiments if the sub-systems are essentially non-interacting 
- such as when they are spatially well-separated, though this does not have to 
be the case. The local systems and operations could involve sub-systems whose 
quantum states and operators are just in different parts of Hilbert space, such as 
for cold atoms in different hyperfine states even when located in the same spatial 
region. Note the distinction between local and localised. As described by Werner 
[5] , if one observer (Alice) is associated with preparing separate sub-system A in 
an allowed quantum state p^ via local operations with a probability Pr, a sec¬ 
ond observer (Bob) could be then advised via a classical communication channel 
to prepare sub-system B in state p^ via local operations. After repeating this 
process for different choices R of the correlated pairs of sub-system states, the 
overall quantum state prepared by both observers via this local operation and 
classical communication protocol (LOCC) would then be the bipartite non- 
entangled state p = PrPr ® p^. Multipartite non-entangled states of the 
form (jSD can also be prepared via LOCC protocols involving further observers. 
As will be seen, the separable or non-entangled states are just those that can 
be prepared by LOCC protocols. 
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2.2.3 Constraints on Sub-System Density Operators 

A key issue however is whether density operators p and in Eq. ([3]) always 

represent possible quantum states, even if the operators p and etc satisfy 

all the standard mathematical requirements for density operators - Hermitiancy, 
positiveness, trace equal to unity, trace of density operator squared being not 
greater than unity. In this paper it will be argued that for systems of identical 
massive particles there are further requirements not only on the overall density 
operator, but also (for separable states) on those for the individual sub-systems 
that are imposed by symmetrisation and super-selection rules. 

2.2.4 Classical Entanglement 

In addition to quantum entanglement there is a body of work (see [7B], [73]1. 
|80j dealing with so-called classical entanglement. Here the states of classical 
systems - such as a classical EM field - are represented via a formalism involving 
linear vector spaces and classical entanglement is defined mathematically. A 
discussion of classical entanglement is beyond the scope of this paper. Although 
there are some formal similarities with quantum entanglement - and even Bell 
type inequalites which can be violated, there are key features that is not analo¬ 
gous to that for composite quantum systems - quantum non-locality being one 
m- In the end, rather than just focusing on similarities in the mathemati¬ 
cal formalisms, classical and quantum entanglement are seen as fundamentally 
different when the physics of the two different types of system - one classical 
and deterministic, the other quantum and probabilistic are taken into account. 
In particular, the key feature of quantum entanglement relating to joint mea¬ 
surement probailities is quite different to the corresponding one for classical 
entanglement. 

2.3 Separate and Joint Measnrements, Reduced Density 
Operator 

In this Subsection we consider separate and joint measurements on systems 
involving several sub-systems and introduce results for probabilities, mean val¬ 
ues for measurements on one of the sub-systems which are conditional on the 
results for measurements on another of the sub-systems. This will require con¬ 
sideration of quantum theoretical conditional probabilities. The measurements 
involved will be assumed for simplicity to be von Neumann projective mea¬ 
surements for physical quantities represented by Hermitian operators fl, which 
project the quantum state into subspaces for the eigenvalue Xi that is measured, 
the subspaces being associated with Hermitian, idempotent projectors H^ whose 
sum over all eigenvalues is unity. These concepts are treated in several quantum 
theory textbooks, for example [n], m- Eor completeness, an account setting 
out the key results is presented in Appendix |6l 
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2.3.1 Joint Measurements on Sub-Systems 

For situations involving distinct sub-systems measurements can be carried out 
on all the sub-systems and the results expressed in terms of the joint probability 
for various outcomes. If is a physical quantity associated with sub-system 
A, with eigenvalues Xf and with Wf the projector onto the subspace with 
eigenvalue Af, VLb is a physical quantity associated with sub-system B, with 
eigenvalues Xf and with the projector onto the subspace with eigenvalue 
Xf etc., then the joint probability PAB..{i,j, ■••) that measurement of VIa leads 
to result Aj , measurement of fls leads to result Xj ,etc is given by 

PAB.Xi,j, ...) = rr(nfnf ...p) (4) 

This joint probability depends on the full density operator p representing the 
allowed quantum state as well as on the quantities being measured. Here the 
projectors (strictly Ilf 0 0 ..., <8> Hf 0 ..., etc) commute, so the order 

of measurements is immaterial. An alternative notation in which the physical 
quantities are also specified is Pab.X^aA^XIb, j; ■■■)■ 

2.3.2 Single Measurements on Sub-Systems and Reduced Density 
Operator 

The reduced density operator p^ for sub-system A given by 


'Pa=Ttb,c,..Xp) (5) 

and enables the results for measurements on sub-system A to be determined for 
the situation where the results for all joint measurements involving the other 
sub-systems are discarded. The probability PaX) that measurement of Ha leads 
to result xf irrespective of the results for meaurements on the other sub-systems 
is given by 


PA{i) = E PAB..{i,j, ■•■) 

= TriUfp) ( 6 ) 

= TrAiUfpA) (7) 

using E nf = 1, etc. Hence the reduced density operator p^ plays the role of 

3 

specifying the quantum state for mode A considered as a separate sub-system, 
even if the original state p is entangled. An alternative notation in which the 
physical quantity is also specified is Pa{XIaX)- 
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2.3.3 Mean Value and Variance 

The mean value for measuring a physical quantity Qa will be given by 

(17^) = 

= TvAin^PA) ( 8 ) 

where we have used ■ 

The variance of measurements of the physical quantity Qa will be given by 

Af 

= TrAi{n^ - {hA)y Pa) (9) 

so both the mean and variance only depend on the reduced density operator 
On the other hand the mean value of a product of sub-system operators 
Q.A O where Q,a, ^b, are Hermitian operators representing 

physical quantities for the separate sub-systems, is given by 

(^VIa ® ^B ® O . ) = EE- .XfXf...PAB..ihJ, ■•■) 

A<^ A® 

1 J 

= Tr (^riA® ^B ® Xlc ® ^P (10) 

which involves the overall system density operator, as expected. 

2.3.4 Conditional Probabilities 

Treating the case of two sub-systems for simplicity we can use Bayes theorem 
(see Appendix [51 Eg. (1155^ 1 to obtain expressions for conditional probabilities 
M- The conditional probability that if measurement of fls associated with 
sub-system B leads to eigenvalue Xf then measurement of flA associated with 
sub-system A leads to eigenvalue is given by 

PABm = Tr{flfflfp)/Tr{flfp) (11) 

In general, the overall density operator is required to determine the conditional 
probability. An alternative notation in which the physical quantities are also 
specified is Pab(17a, j)- 

As shown in Appendix [6] the conditional probability is given by 

PABm = Trifimondm, xf)) m 

where 

AonA,xf) = nf pfif/Tr(nf p) (13) 
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is the so-called conditioned density operator, corresponding the quantum state 
produced following the measurement oi fls that obtained the result Xf. The 
conditional probability result is the same as 

PABm = rr(nf Af)) (14) 


which is the same as the expression (jS]) with p replaced by 7>condi^B, Xf). This 
is what would be expected for a conditioned measurement probability. 

Also, if the measurement results for fig are not recorded the conditioned 
density operator now becomes 


Pcondi^B) 'y PBjj^Pcondi^Bj Xj ) 

= Enfpfl? 


(15) 


This is still different to the original density operator p because a measurement of 
Qb has occured, even if we dont know the outcome. However, the measurement 
probability for is now 

PAB{i\Anyj) = rr(nf 

= TriUfp) (16) 

= PA{i) (17) 

where we have used the cyclic properties of the trace, and 

nf = 1- The results in Eqs. ([T6|) and (fT7|) are the same as the mea¬ 
surement probability for if no measurement for Hb had taken place at all. 
This is perhaps not surprising, since the record of the latter measurements was 
discarded. Another way of showing this result is that Bayes Theorem tells us 
that PAB{i\j)PBij) = Xj PAB^iJ) = PA{i), since Xj PAB{i,j) is the prob¬ 
ability that measurement of 14^4 will lead to Xf and measurement of Hb will 
lead to any of the Xf . This result is called the no-signalling theorem [14]. 

Also, as PAB{i\Any j) = Tr{Ilfp^g^^{nB)) we see from 0 that 

PA = TrBiPcond(.^B)) (18) 


showing that the trace over B of the conditioned density operator for the state 
obtained by measuring any observable XIa and then discarding the results just 
gives the reduced density operator for sub-system A. 


2.3.5 Conditional Mean and Variance 

As explained in i^pendix|6l to determine the conditioned mean value of A after 
measurement of XI has led to the eigenvalue Xi we use 'p^g^^{Xl,i) rather than 
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p in the mean formula = Tr{Kp) and the result is given in terms of the 

conditional probability P{Aj\ni). Here we refer to two commuting observables 
and include the operators in the notation to avoid any misinterpretation. Hence 

( a ) = Tr(Ap^ond(Mp)) 

= (19) 

3 


For the conditioned variance of A after measurement of H has led to the 
eigenvalue \i we use 'Pcondi^^ *) rather than p and the conditioned mean ( A )_ rather 

than ^A^ in the variance formula ^AA^^ = Tr((A — ^A^)^p). Hence 

(aA^) = Tr((A-(A) )2p_,(H,*)) 

\ / U,,t \ / 12,2 

= ( 20 ) 

3 


If we weighted the conditioned mean by the probability P{n,i) that mea¬ 
suring n has led to the eigenvalue Xi and summed over the possible outcomes 
Xi for the fl measurement, then we obtain the mean for measurements of 
A after un-recorded measurements of H have occured. From Bayes theorem 


P{A, j\it,i)P{rt,i) = P{A,j) so this gives the unrecorded mean 

I 



as 


= Y^^3P(^^d) 

3 

= (a) (21) 


which is the usual mean value for measurements of A when no measurements 
of H have occured. Note that no such similar result occurs for the unrecorded 
variance ^AA^^_ 


(aa^)^ = E(aa"). P{n,i) 

^ (aa^) 


( 22 ) 
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2.4 Non-Entangled States 

In this Subsection we will set out the key results for measurements on non- 
entangled states. 


2.4.1 Non-Entangled States - Joint Measurements on Sub-Systems 


In the case of the general non-entangled state we hnd that the joint probability 
PAB.ihJ, (23) 

R 

where 

P«(z) = Tr(nfp^) P«(j)=Tr(nf p|) .. (24) 

are the probabilities for measurement results for Ha, Hb, ... on the separate sub- 
systems with density operators etc and the overall joint probability is 

given by the products of the probabilities E 4 (*), (j), •• for the measurement 

results Xf, Xf, ... for physical quantities Ha, Hb, ... if the sub-systems are in 
the states p^, pf, etc. Note that here E 4 (*), (j), are given by quantum the¬ 

ory formulae for the sub-systenm states. For simplicity only quantized measured 
values will be considered - the extension to continuous values is straightforward. 
Thus the results for the probabilities of joint measurements when the system is 
in a separable quantum state are determined by the measurement probabilities 
in possible quantum states for the sub-systems, combined with a classical prob¬ 
ability for creating the particular pair of sub-system quantum states. Note the 
emphasis on ’’possible” - some of the separable states described in [13] are not 
possible. 

Furthermore, if we consider measurements of the physical quantity Ha 0 Hb 
then for a separable state the mean value for measurement of this quantity is 
given by 

A B 

^Ha®^b) =Fr (Ha®52bp) (Ha) (Hb) (25) 

R R R 


where ^Ha^ = Tva ^HaPbJ and = Tvb are the mean values 

of Ha and Hb for the sub-system states p^ and pf respectively. If ( Ha < 8 ) Hb ) = 


^Ha^ then the state is said to be uncorrelated. Separable states are corre- 

lated except for the case where p^^^ = p 0 p , but the correlation is essentially 
non-quantum and attributable to the classical probabilities Pr. However, for 
separable states the inequality | ^Ha 0 H^^ p < ^hJ^Ha 0 H^Hb^ applies, so 

that if I ^Ha 0 p > ^H^Ha 0 H^Ha^ then the state is entangled. 

In the simple non-entangled pure state situation in Eq. (HD the joint proba- 
bilty only involves a single product of sub-system probabilities 


PAB..{i,j, •■•) = PA{i)PB{j)--- 


( 26 ) 


30 


where 


(27) 


PA(^) = ($a| nf \<^a) PbU) = ($b| nf |$b) 

just give the probabilities for measurements in the separate sub-systems. 

This key result (|23l) showing that the joint measurement probability for a 
separable state only depends on separate measurement probabilities for the sub¬ 
systems, together with the classical probability for preparing correlated product 
states of the sub-systems, does not necessarily apply for entangled states [5]. 
However the key quantum feature for composite systems of non-separability for 
joint measurement probabilites applies only to entangled states. This strange 
quantum feature of entangled states has been regarded as particularly unusual 
when the sub-systems are spatially well-separated (or non-local) because then 
measurement events can become space-like separated. This is relevant to quan¬ 
tum paradoxes such as Einstein-Poldolsky-Rosen (EPR) and Bell’s theorem 
which aim to show there could be no causal classical theory explaining quan¬ 
tum mechanics n, m- Measurements on sub-system A of physical quantity 
Qa affect the results of measurements of Hs at the same time on a distant 
sub-system R, even if the choice of measured quantity is unknown to the 
experimenter measuring ^Ia- As will be shown below, a similar result to (I23p 
also occurs in hidden variable theory - a classical theory - so non-separability for 
joint measurements resulting from entanglement is a trnly non-classical feature 
of composite systems. 

2.4.2 Non-Entangled States - Single Sub-System Measurements 

Eor the general non-entangled state, the reduced density operator for sub-system 
A is given by 

(28) 

R 

A key feature of a non-entangled state is that the results of a measurement on 
any one of the sub-systems is independent of the states for the other subsystems. 
Prom Eqs. © and ((^ the probability PA{i) that measurement of ft a leads to 
result Xf is given by 

Pa(*) = E^«^a(*) (29) 

R 

where the reduced density operator is given by Eq. (E51) for the non-entangled 
state in Eq. This result only depends on the reduced density operator 
which represents a state for sub-system A and which is a statistical mixture of 
the sub-system states p^, with a probability Pr that is the same for all sub¬ 
systems. The result for the measurement probability Pa(*) is just the statistical 
average of the results that would apply if sub-system A were in possible states 
Pfl. For all quantum states the final expression for the measurement probabil¬ 
ity PAii) only involves a trace of quantities Ilf, p^ that apply to sub-system 
A, but for a non-entangled state the reduced density operator p^ is given by 
an expression (l28l) that does not involve density operators for the other sub¬ 
systems. Thus for a non-entangled state, the probability PA{i) is independent 
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of the states pj^, associated with the other sub-systems. Analogous results 
apply for measurements on the other sub-systems. 

2.4.3 Non-Entangled States - Conditional Probability 

For a general non-entangled bipartite mixed state the conditional probability is 
given by 

PABm = E Pa{^Pb{3 )/E Pr PbU) (30) 

R R 

which in general depends on fls associated with sub-system B and the eigen¬ 
value Xf. This may seem surprising for the case where A and B are localised 
sub-systems which are well separated. Even for separable states a measurement 
result for sub-system B will give immediate information about a totally sepa¬ 
rated measurement on sub-system A - which is space-like separated. However it 
should be remembered that the general separable system can still be a correlated 
state, since each sub-system density operator for sub-system B is matched 
with a corresponding density operator for sub-system A. Results at A can 
be correlated with those at B, so the observer at A can potentially infer from 
a local measurement on the sub-system A the result of a local measurement on 
sub-system B. It is therefore not necessarily the case that measurement results 
for A are independent of those for B. However, as we will see below, such corre¬ 
lations (usually) hve a classical interpretation. Result (|5(I)) is not a case of the 
’’spooky action at a distance” that Einstein [T] referred to. 

However, for a non-entangled pure state where p = p ® p we do find that 

Pab{i\3) = PaH) (31) 

where PaII) = Trlllf'p^). For separable pure states the conditional probability 
is independent of associated with sub-system B and the eigenvalue X^. 
Also of course PAB{i\3)PB{3) = PA{i) is true for separable states since 

it applies to general bipartite states. Hence if the measurement results for 
are discarded then the probability distribution for measurements on will be 
determined from the conditioned density operator 'Pcondi^B) and just result in 
Pa(j) - as in shown in Ea. dTTI) for any quantum state. 

2.4.4 Non-Entangled States - Mean Values and Correlations 

For non-entangled states as in Eq. © the mean value for measuring a physical 
quantity Ha (E) Hb (E> Hq (E> ■■■, where Ha, Hb, Hq, ■■ are Hermitian operators 
representing physical quantities for the separate sub-systems can be obtained 
from Eqs.dSl) and (ITOll and is given by 

(Ha 0 Hb 0 He 0 . ( 32 ) 
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where 


K 


(33) 


Qk) =Tr{nKPR), {K = A,B,...) 


R 


is the mean value for measuring ^Ir in the K sub-system when its density 
operator is Since the overall mean value is not equal to the product of 
the separate mean values, the measurements on the sub-systems are said to be 
correlated. However, for the general non-entangled state as the mean value is 
just the products of mean values weighted by the probability of preparing the 
particular product state - which involves a LOCC protocal, as we have seen - 
the correlation is classical rather than quantum |14) . In the case of a single 


p 0 ... we have = 


product state where 'p = 'p 0 p" 

17^4^ "''^hich is just the product of mean values for the 

separate sub-systems, and in this case the measurements on the sub-systems 
are said to be uncorrelated. For entangled states however the last result for 
17.4 0 0 f7c 0 ■) does not apply, and the correlation is strictly quantum. 


2.5 Local Hidden Variable Theories 

In a general local hidden variable theory as envisaged by Einstein et al [I] and 
Bell [5], physical quantities associated with the sub-systems are denoted 17 a, 
I7b etc, which are real numbers not operators. Their values are assumed to 
be Aj , Xj etc - having the same ranges as in quantum theory, since HVT 
does not challenge the quantization feature. In the realist viewpoint of HVT 
all the physical quantities have definite values at any time, the probabilities 
for measuring these values being determined from a set of hidden variables 
which are themselves given by a probability function P(^) for each state 
preparation process. Measurement is not required for the values for physical 
quantities to be created, as in quantum theory, nor do the hidden variables 
change as a result of the act of measurement itself (though they may change 
as a result of local interactions of the system with the measurement apparatus 
[S2], [53] ■ As in classical physics, ideal measurement is assumed not to change 
the state of the system - the hidden variables would only change in accord 
with the (as yet unknown) dynamical equations that govern their evolution. 
The hidden variables are regarded as the elements of reality that constitute the 
fundamental way of describing the system [1]. There may be just a single hidden 
variable or a set, and the hidden variables could be discrete or continuous - these 
details do not matter in a general HVT. In the original treatment of Bell [3] the 
hidden variables uniquely determine the actual values that physical quantities 
would have when measured. However, in a so-called ’’fuzzy” hidden variable 
theory [82|, [53], [84], [85], [15], [6] (see also Section 7.1 of [IS]) the values 
for Qa, etc are determined prohahilistically from the hidden variables, the 
probability functions being classical and allow for the hidden variables not being 
known - just as in classical statistical mechanics, where the unknown (but real) 
positions and momenta of the classical particles are described via probabilities. 
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The probabilistic treatment of the hidden variables attempts to replicate the 
probabilistic nature of quantum theory. For our purposes we will consider only 
local hidden variable theories (LHV) - this is sufficient to demonstrate key results 
such as the Bell inequalities. For local hidden variable theories although the 
hidden variables ^ are global, they act locally even for spatially separable sub¬ 
systems. For particular hidden variables ^ the probability that VIa has value 

will be given by Pa(*,C) the probability that fls bas value Xf will be 
given by The LHVT joint probability for measurement outcome 

for Ha, XIb, etc will be given by 

PAB.Xh3,-)= ! d^P{0PA{i,S,)PB{3,0- (34) 

Here P{^)d^ is the probability that the hidden variables are in the range 
around the HV being assumed continuous - which is not a requirement [J. 
The probabilities satisfy the usual sum rules for all outcomes giving unity, thus 
^■Pa(*, 0 = 1, etc., f d^P(^) = 1. The sub-system probabilities Pa(*,C)) 
PbUjO etc only depend on the hidden variables Bell inequalities are con¬ 
straints derived on the basis of the assumption (IMl) . and if violated therefore 
falsify all LHV theories. 

The formal similarity between the hidden variable theory expression for the 
joint probability (|M1) and the quantum expression (1^ for a separable state is 

noticable. We could map f ^ R, P{f) -A Pr, J df ^ PA{i,£,) -A PaX) and 

R 

PB{j,f) -A PbU)- Tbe Werner preparation process [5] would then determine 
the setting for the hidden variables If a hidden variable theory underpinned 
quantum theory, it follows that the quantum probabilities PaX) and P^ij) 
would always be equivalent to hidden variable probabilites PbUjO 

for each of the sub-systems (it would not be consistent to only have this apply 
to one of the sub-systems and not the other). From the expression (1341) for the 
joint probability general hidden variable theory expressions for the mean value 
{Qa X XIb)RYR for the product of the measurement results for observables Ha 
and H b for subsystems A, B respectively (see (1551) below) can be obtained that 
are analogous to the quantum expression (15^ for a separable state. There is of 
course no independent fully developed classical HVT that can actually predict 
the Pa(L'C), PB(j,0-etc. 

However, as we will see both the HVT (see m for a proof) and the quantum 
separable state predictions are consistent with Bell Inequalities, and it therefore 
requires a quantum entangled state to violate Bell inequalities and to demon¬ 
strate failure of the LHV theory model (1341) . Naturally it follows that such 
quantum entangled states cannot be described via a LHV theory. Hence the 
experimental violation of Bell inequalities would also show that the particular 
quantum state must be entangled. 

2.5.1 LHV - Mean Values and Correlation 

The actual values that would be assigned to the physical quantities XIa, XIb etc 
will depend on the hidden variables but can be taken as the mean values of the 
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possible values Xf,Xf etc. We denote these mean values as (n^(^)), {^b{0) 
etc where 

{^Km = E>^k PKik,0 {K = A,B,...) (35) 

These expressions my be compared to Eg. (1551) for the mean values of physical 
quantities flyi, ^Ib etc in quantum separable states. 

We can then obtain an expression for the mean value in HVT of the physical 
quantity VLa x x ^Ic x ..., where 17^, XIb, etc. are physical quantities for the 
separate sub-systems. This is obtained from Eas. (|551) and (I55|) and is given by 

X X f7c X .)^HV =Id^P{0 i^AiO) {^b{0) {^€{0) (36) 

This may be compared to Eg. (1551) for the mean value of the physical quantity 
flA 0 fls 0 Xlc 0 .. in quantum separable states. 


2.5.2 LHV - GHZ State 


The GHZ state [SB] , [SI] is an entangled state of three sub-systems A, B and C, 
each of which is associated with two quantum states |-|-1) and | —1). Each sub¬ 
system has three physical quantities, which correspond to Pauli spin operators 
and The quantum states |-|-1) and | —1) are eigenstates of with 
eigenvalues -|-1 and —1 respectively. Note that the eigenvalues of the other two 
Pauli spin operators are also -|-1 and —1. The GHZ state is defined by 


\'^)ghz — (I+1)a I+1)b I+1)c + I“1)b I“1)c)/'^ 


(37) 


The GHZ state provides a clear example of an entangled quantum state 
which cannot be described via local hidden variable theory m, m- In a non- 
fuzzy version of the LHV model each of the nine physical quantities a^, 
a^, <7^, <7y, erf, erf, erf, erf will be associated with hidden variables that 
directly specify the values -|-1 and —1 that each one of these physical quantities 
may have. We denote these hidden variables as , where K = A,B,C and 
a = x,y,z and we have = -|-1 or —1. With this direct specification of the 
physical values Eg. (1551) just becomes (crf-(M^)) = and Eg. (1551) becomes 

= M¥ Mf .We can then derive a contradiction with 
LHV ^ ^ 

quantum theory regarding the LHV description of the GHZ state. 

Firstly, using the Pauli spin matrices for the |-|-1) and | —1) basis states 


(38) 


it is straightforward to show that the GHZ state satisfies three eigenvalue equa¬ 
tions 


(7^ X erf X erf. 


■ 0 

1 


1 

o 


■ 1 0 ■ 

1 

0 

[dy] = 

i 0 

[dx\ = 

0 -1 


—B'-^C i»Tf\ 

'^A' — B'^C i»Tf\ 
^y^x^y \^)gHZ 
'^A' — B'^C i»Tf\ 
^y^y^x \^)gHZ 


{-^)mGHZ 

{-^)mGHZ 

{-^)\'^)ghz 


(39) 
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Hence in LHV the three quantities a^a^ay, a^a^ay and must all 

have value —1 in the GHZ state, so that as the values for these quantities are 
just the products of the values for each of the factors we get three equations 

= -1 = -1 = -1 (40) 

Secondly, if we apply all three operators to the GHZ state we find 

another eigenvalue equation 

^x'^x^x \^)ghz ~ i~^^)\^)GHZ ("^1) 

which leads to 

Mf = +1 (42) 

However, if we multiply the three equations in Ea. (l40l) together and use = 

+ 1 we find that = — 1, in direct contradiction to the last equation. 

Thus the assignment of hidden variables for all the physical quantities fails 
to describe the GHZ state. As we will see in the next SubSection, there are tests 
involving the violation of Bell Inequalities that are satisfied by some entangled 
states which allow a demonstration of the failure of more general local LHV 
theories, even allowing for correlations that are less than ideal. 

The assumption of non-fuzzy LHV theories is not essential for the GHZ 
arguement in the case of the ideal GHZ state (1571) . This is because one may use 
the correlations of (1551) to establish a precise prediction of one of the spins at 
A, by measuring the spins at the other two locations. The assumption of local 
realism (on which the LHV theory is based) then establishes a precise value for 
the hidden variable [86], [87], [88]. In a more realistic scenario where the GHZ 
correlations are not perfect, the ’’elements of reality” established this way becme 
fuzzy, and in that case Mermin’s Bell inequality [55] can be used to establish a 
contradiction with LHV models. 

2.6 Paradoxes 

The EPR and Schrodinger Gat paradoxes figured prominently in early discus¬ 
sions about entanglement. Both paradoxes involve composite systems and the 
consideration of quantum states which are entangled Both these paradoxes 
reflect the conflict between quantum theory, in which the values for physical 
quantities only take on definite values when measurement occurs and classi¬ 
cal theory, in which the values for physical quantities always exist even when 
measurement is not involved. The latter viewpoint is referred to as realism. 
Quantum theory is also probabalistic, so although the possible outcomes for 
measuring a physical quantity can be determined prior to measurement, the 
actual outcome in a given quantum state for any measurement is only known in 
terms of a probability. However, from the realist viewpoint, quantum theory is 
incomplete and a future theory based around hidden variables would determine 
the actual values of the physical quantities, as well as the quantum probabilities 
that particular values will be found via measurement. 
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Whilst the EPR and Schrodinger Cat paradoxes are of historical interest 
and have provoked much debate, it was the formulation of the Bell inequalities 
(which are described in the next SubSection l2.7l) and the conditions under which 
they could be violated that provided the first clear case of where the predictions 
of quantum theory could differ from those of hidden variable theories. It then 
became possible to carry out actual experiments to distinguish these two funda¬ 
mentally different theories. The actual experimental evidence is consistent with 
quantum theory and (apart from a small number of remaining loopholes) rules 
out local hidden variable theories. 

2.6.1 EPR Paradox 

In the original version of the EPR paradox, Einstein et al [T] considered a two- 
particle system A, B in which the particles were associated with positions xa, 
Xb and momenta pA, Pb- They envisaged a quantum state in which the pairs of 
physical quantities xa, xb or pa, Pb had highly correlated values - measured or 
otherwise. To be specific, one may consider a simultaneous eigenstate |$) of the 
two commuting operators xa — xb and pa +Pb, where {xa — xb) |$) = 2a: 1$) 
and {pa + Pb) I'h) = 0 |$). This state is an example of an entangled state, as 
may be seen if it is expanded in terms of position eigenstates \xaxb)- If the 
system is in state |$) then from standard quantum theory if A had a mean 
momentum p then B would have a mean momentum —p. Alternatively, if 
A had a mean position x then B would have a mean position —x. Then if 
the eigenvalue 2x is very large so that the two particles will be well-separated 
(in quantum theory their spatial wave functions would be localised in separate 
spatial regions) it follows that if the position of B was measured then the position 
of A would be immediately known, even if the particles were light years apart. 
On the other hand, if the momentum of B was measured instead, then the 
momentum of A would immediately be known. From the realist point of view 
both A and B always have definite positions and momenta, even if these are not 
known, so all these measurements do is reveal these (hidden) values. It would 
seem then that measurements of position and momentum on particle B could 
lead to a knowledge of the position and momentum at a far distant particle 
A, perhaps with an accuracy that would violate the Heisenberg Uncertainty 
Principle (HUP). As we will see, this is not the case when quantum theory is 
applied correctly. However, what Einstein et al pointed out as being particularly 
strange was that the choice of whether the momentum or position of B was 
measured (and found to have a definite value) would instantly determine which 
of the position or momentum of A would then have a definite value - even if 
A and B were separated by such a large distance that no signal could have 
been passed from B to A regarding which quantity was measured. Einstein 
referred to this as ’’spooky action at a distance” to highlight the strangeness of 
what came to be referred to as entangled states. Thus a somewhat paradoxical 
situation would seem to arise. Einstein stated that this did not demonstrate 
that quantum theory was wrong, only that it was incomplete. 

The EPR argument assumes local realism, to justify that the posibility of 
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an exact prediction of the postion of the far-away particle A (based on the 
measurement of the position for the particle B) implies the realist viewpoint 
that the position of particle A was predetermined. The same argument applies 
to the momentum of particle A, and hence EPR conclude that both the position 
and momentum of particle A are precisely predetermined - in conflict with the 
Heisenberg Uncertainty Principle derived from quantum mechanics. Since the 
argument is based on the assumption of local realism, the modern interpretation 
of the EPR analysis is that it reveals (for the appropriate entangled state) the 
inconsistency of local realism with the completeness of quantum mechanics. 

Discussions of the EPR paradox [T] in terms of hidden variable theories has 
been given by numerous authors (see m. ca, M, [m for example). The 
papers and reviews by Reid et al [76], m, na, give a full account taking into 
consideration the ’’fuzzy” version of local HVT (LHV) and determining the pre¬ 
dictions for the conditional variances for xa and pA based both on separable 
quantum states and states described via local HVT. This treatment success¬ 
fully quantifies the somewhat qualitative considerations described in the pre¬ 
vious paragraph. If the position for particle B is measured and the result is 
X, then the original density operator 'p for the two particle system is changed 
into the conditional density operator 'p^^^^{xB,x) = 'pU!^/Tr{Il^^), where 
Hf = {\x) {x\)b is the projector onto the eigenvector \x)^ (the eigenvalues x are 
assumed for simplicity to form a quasi-continuum). Similarly, if the momentum 
for particle B is measured and the result is p, then the original density operator 
p for the two particle system is changed into the conditional density operator 
Pcoud{PB,p) = H^ pn ^/ rr ( n ® p ), where H® = {\p) {p\)b is the projector onto 
the eigenvector \p)b (the eigenvalues p are assumed for simplicity to form a 
quasi-continuum). Here we outline the discussion based on quantum separable 
states. Conditional variances for position and momentum for sub-system A are 
considered based on measurements for sub-system B of position. It is shown 
that for these conditional variances the Heisenberg uncertainty principle still 
applies. The same conclusion is obtained if the measurements on sub-system B 
had been the momentum. As the experimenter on sub-system A could not know 
whether the measurement on sub-system B was on position or momentum, the 
action at a distance feature of quantum entanglement is confirmed. 

The question is whether the conditional variances for measuring 

XA for sub-system A having measured xb for sub-system B, and for 

measuring pA for sub-system A having measured pb for sub-system B violate 
the Heisenberg Uncertainty Principle m 

{ASbs. (Apbf. < (43) 

where the measurements on sub-system B are left unrecorded. If this inequality 
holds we have an EPR violation. 

For separable states the conditional probability that measurement of xa on 
sub-system A leads to eigenvalue xa given that measurement oi xb on sub- 
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system B leads to eigenvalue xs is obtained from Eg. (15(11) as 


P{xa-,xa\xb,xb) = ^PrPa{xa,xa)P§{xb,xb)/^PrPb{xb,xb) (44) 

R R 


where 


P^{xa,xa) = TrA{Btj>i) P^{xb.xb) = TrB(nf^pi) (45) 

are the probabilities for position measurements in the separate sub-systems. 
The probability that measurement of xr on sub-system B leads to eigenvalue 
Xb is 

P{xb,xb) ='^PrPb{xb,xb) (46) 

R 

The mean result for measurement of xa for this conditional measurement is 
from Ea. (fT51) 


= '^xaP{xa,xa\xb,xb) 

XA 

= ^Pr {xa)rPE{xb,xb)/P{xb,xb) (47) 

R 

where 

{xa)r = ^ xaPa{xa, xa) (48) 

XA 

is the mean result for measurement of ir^ when the sub-system is in state 
The conditional variance for measurement of ir^ for the conditional mea¬ 
surement oixB on sub-system B which led to eigenvalue is from Ea. (|20l) 

= Y^{xa-{xa)^^^^^)"^ P{xa.xa\xb.xb) 

XA 

= '^Pr {^^A)^s,XB^BiXB,XB)/P{XB,XB) (49) 

R 

where 

(^^a)^b,xb = PAixA,XA) 

XA 

is a variance for measurement of ir^ for when the sub-system is in state but 
now with the fluctuation about the mean (xa)'- for measurements condi- 
tional on measuring xr- 

However, for each sub-system state R the quantity is minimised 

if {xa)-Sjj xb replaced by the unconditioned mean {xa)r just determined from 
Pfl. Thus we have an inequality 

(^^a)Z,xb ^ ( 50 ) 
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where 


(51) 


(Ax^) = y^jxA - {xA)f Pa{xa,xa) 

XA 

is the normal variance for measurement of xa for when the sub-system is in 
state 'p"^. 

Now if the measurements oixB are unrecorded - as would be the case from 
the point of view of the experimenter on spatially well-separated sub-system A 
when measurements on this sub-system take place at the same time - then the 
conditioned variance is 

XB 

= J:J2Pr{Ax\)^^^^pB{xb,xb) (52) 

XB ^ 

which in view of inequality (15(11) satisfies 

P^Cxb.xb) 

XB ^ 

= Y^PniAxlf (53) 

R 

using '^P§{xb,xb) = 1- Thus the variance for measurement of position xa 

XB 

conditioned on unrecorded measurements for position xb satisfies an inequality 
that only depends on the variances for measurements of xa in the possible 
sub-system A states 

Now exactly the same treatment can be carried out for the variance of mo¬ 
mentum PA also conditioned on unrecorded measurements of measurements for 
momentum xb ■ Details are given in Appendix [6] We have with 





PB 


(Ap\)^ 

= '^(PA - {pa)p^^pJ'^ P{pa,Pa\pb,Pb) 



PA 


{P^)pB,PB 

= '^PaP{.Pa,Pa\pb,Pb) 



PA 


the inequality 

R 

(54) 

with 


{^Pa) = - {PA)f Pa{va,Pa) 

(55) 


PA 


is the normal variance 

for measurement of pA for when the sub-system 

is in 

state 'p^. 
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We now multiply the two conditional variances, which it is important to 
note were associated with two different conditioned states based on two different 
measurements - position and momentum - carried out on sub-system B. 

^ ^ Ps {Aplf (56) 

R S 

However, from the general inequality in Eq.([T77|) 

Pr^rJ2 Pr Pr ^ Pr VCrP^ (57) 

we then have 

- (5 

= (^E Pr \/ X \J(Ap^}^'j (58) 

But we know from the HUP that for any given state 'p'^ that (( Ap^)^ > 
so for the conditioned variances associated with a separable state 

(59) 

showing that for a separable state the conditioned variances (Ax^)^ and 
still satisfy the HUP. It is important to note that these variances were 
associated with two different conditioned states based on two different measure¬ 
ments - position and momentum - carried out on sub-system B, the results of 
which the observer for sub-system A would be unaware of. Thus if the EPR 
violations as defined in Ea. (H5)) are to occur then the state must be entangled. 
Progress towards experimental confirmation of EPR violations is reviewed in 
Refs. [TS], 0. 

In [85) an analogous treatment based on local hidden variable theory (LHV) 
also shows that the HUP is satisfied for the conditioned variances. The details 
of this treatment will not be given here, but the formal similarity of expressions 
for conditional probabilities in LHV theories and for separable states indicates 
the steps involved. 

The EPR paradox is not confined to position and momentum measurements 
on two sub-systems. A related paradox [89] occurs in the case of measurements 
on spin components Sai and Sa 2 - with a = x,y, z - associated with two sub¬ 
systems 1 and 2. The spin operators also satisfy non-zero commutation rules 
(see paper H for details) 

[•Sal, 5/3 i] = iS^l [Sa2, -5/32] = iSj2 (60) 
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where a, /3 ,7 are x,y, z in cyclic order. Different spin components for each sub¬ 
system do not have simultaneous precise measurements leading to Heisenberg 
Uncertainty principle relations involving the variances and mean values 

(a 5 L) (A5I2) > Jl (§72) P ( 61 ) 

As in the case of position and momentum a special state of the combined system 
has interesting features. For the case where the spin quantum number of each 
sub-system is 1/2 the measured values for any spin component of either system 
is either - 1 - 1/2 or In terms of eigenstates for S^i and 8^2 we consider the 

state 

(k, +)i ® k, -)2 - \x, -)i ® |a;, -f) 2 ) (62) 

This is actually one of the Bell states. In this form it shows that measurements 
of the X components of the spins are perfectly correlated, so that for example 
if the measurement of 8 x 2 for sub-system 2 results in the value — 1 / 2 , then a 
subsequent measurement of 8 xi for sub-system 1 must result in the value - 1 - 1 / 2 . 
However, the same state can be expressed in terms of eigenstates for and 
8 y 2 as 

1^") = (l2/>+)i ® \y^-)2 - \y^-)i ® \y^+)2) (63) 

and analogous statements regarding measurement correlations apply if the mea¬ 
surements were for 8 y 2 on sub-system 2 with a subsequent measurement of Syl 
on sub-system 1 . If the two sub-systems were well-separated it might be ex¬ 
pected that first measuring 8 x 2 for sub-system 2 would determine the result 
of measuring for sub-system 1 , and then measuring 8 y 2 for sub-system 2 
would determine the result of measuring 8 yi for sub-system 1 - and as the sec¬ 
ond { 8 y 2 ) measurement on far distant sub-system 2 should not affect the former 
measurement on sub-system 1 this would appear to result in precise measured 
values for 8 xi and 8 yi on sub-system 1, which conflicts with the Heisenberg 

Uncertainty principle requirement that (^A 8 xi^ P- 

However we can consider the variances for 8 xi and which are conditional 
on measurements for 8 x 2 and 8 y 2 for sub-system 2 and show that for separable 
states of the two sub-systems we have 




'^Sy2 

IV 


(64) 

Thus if we find that 







N A 
<73 

<1 

'^Sy2 

<il 


(65) 


then we have an example of a spin EPR violation. Such a violation requires that 
the quantum state is entangled. The derivation of the result (IMl) for separable 
states is set out in Appendix [51 
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An effect related to the EPR paradox is EPR Steering. As we have seen, the 
measurement of the position for particle B changes the density operator and con¬ 
sequently the probability distributions for measurements on particle A will now 
be determined from the conditional probabilities, such as PABixA,XA\xB,XB) 
or PABiPA,PA\xB,XB).AiiXLS measurements on B are said to steer the results 
for measurements on A. Steering will of course only apply if the measurement 
results for xb are recorded, and not discarded. A discussion of EPR Steering 
(see [in], Hm) is beyond the scope of this article. 

2.6.2 Schrodinger Cat Paradox 

The Schrodinger Cat Paradox 0, mni relates to composite systems where one 
sub-system (the cat) is macroscopic and the other sub-system is microscopic 
(the radioactive atom). The paradox is a clear consequence of quantum theory 
allowing the existence of entangled states. Schrodinger envisaged a state in 
which an alive cat and an undecayed atom existed at an initial time, and because 
the decayed atom would be associated with a dead cat, the system after a 
time corresponding to the half-life for radioactive decay would be described in 
quantum theory via the entangled state 

1^) = ■^(|e)Aiom \Alive)cat + Id)Atom \Dead)cat) (66) 

in an obvious notation. The quantum state defined by (l66l) represents the 
knowledge that an observer outside the box would have about the combined 
atom-cat system one hour after the live cat was placed in the box along with 
an undecayed atom. The combined system is in an enclosed box, and opening 
the box and observing what is inside constitutes a measurement on the system. 
According to quantum theory if the box was opened at this time there would be 
a probability of 1/2 of finding the atom undecayed and the cat alive, with the 
same probability for finding a decayed atom and a dead cat. From the realist 
viewpoint the cat should be either dead or it should be alive irrespective of 
whether the box is opened or not, and it was regarded as a paradox that in 
the quantum theory description of the state prior to measurement the cat is in 
some sense both dead and alive. This paradox is made worse because the cat 
is a macroscopic system - how could a cat be either dead or alive at the same 
time, it must be one or the other? From the quantum point of view in which 
the actual values of physical quantities only appear when measurement occurs, 
the Schrodinger cat presents no paradox. The two possible values signifying 
the health of the cat are ” alive” and ” dead”, and these values are found with a 
probability of 1/2 when measurement takes place on opening the box, and this 
would entirely explain the results if such an experiment were to be performed. 
There is of course no paradox if the quantum state is only considered to represent 
the observer’s information about what is inside the box. If the box is closed 
then at one half life after the cat was put into the box, the state vector 
enables the outside observer to correctly assess the probability that the cat will 
be alive is 1/2. If the box is then opened and the cat is found to be dead, then 
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the observer’s information changes and the state vector for the cat-atom system 
is now 

= \9)Atom Cat ( 67 ) 


In this interpretation of quantum states, the notion of there being some sort of 
underlying reality that exists prior to measurement is rejected. It is only this 
notion that such a reality must exist - perhaps described via hidden variables 
- that leads to the paradox. EPR paradoxes can also be constructed from the 
entangled state (1551) . as outlined in Refs. m. m- 

In recent times, experiments based on a Rydberg atom in a microwave cavity 
[93] involving states such as (I5S1) have been performed showing that entangle¬ 
ment can occur between macroscopic and microscopic systems, and it is even 
possible to prepare states analogous to Alive)f^^^ + lDead) in the macro¬ 

scopic system itself. In such experiments the different macroscopic states are 
large amplitude coherent states of the cavity mode. Coherent states are pos¬ 
sible for microwave photons as they are created from classical currents wiith 
well-defined phases. A coherent superposition of an alive and dead cat within 
the cat sub-system itself can be created by measurement. The entangled state 
in (l66l) can also be written as 


|vb) 


i\^)Atom + \9)Atom) Cat + \Dead)cat) 

+ ^ (|e)Atom - \9)Atom) Cat “ \Dead)cat)} (68) 


so that measurement on the atom for an observable in which the superposition 
states {\^)Atom ^ \9)Atom) eigenstates for this observable would result 

in the cat then being in the corresponding.macroscopic superposition states 
Alive) ± \Dead)(jg^.f) of an alive and dead cat. 


2.7 Bell Inequalities 

Violations of Bell’s Inequalities represent situations where neither hidden vari¬ 
able theory nor quantum theory based on separable states can account for the 
result, and therefore provide a clear case where an entangled quantum state is 
involved. 


2.7.1 Local Hidden Variable Theory Result 

A key feature of entangled states is that they are associated with violations of 
Bell inequalities [5] and hence can exhibit this particular non-classical feature. 
The Bell inequalities arise in attempts to restore a classical interpretation of 
quantum thory via hidden variable treatments, where actual values are assigned 
to all measureable quantities - including those which in quantum theory are 
associated with non-commuting Hermitian operators. In this case we consider 
two different physical quantities 11^ for sub-system A, which are listed Ai, 
A 2 , etc , and two fls for sub-system B, which are listed Bi, B 2 , etc. The 
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corresponding quantum Hermitian operators ^Ia, ^b, etc are Ai, A 2 and, Bi, 
i? 2 - The Bell inequalities involve the mean value {Ai x Bj) of the product 
of observables Ai and Bj for subsystems A, B respectively, for which there are 
two possible measured values, +1 and —1. For simplicity we consider a local 
HVT. In a local hidden variable theory.(LHV) we see using (IMl) that the mean 
values {Ai x Bj)j^^y are given by 


{A, X B,)^^y = I d^P{0 {Am {Bm) (69) 

where and {Bj{^)) (as in Ea. (l551) l are the values are assigned to Ai and 

Bj when the hidden variables are and P(^) is the hidden variable probability 
distribution function. If the corresponding quantum Hermitian operators are 
such that their eigenvalues are +1 and —1 - as in the case of Pauli spin operators 

- then the only possible values for {Ai{^)) and {Bj{^)) are betweem +1 and 

— 1, since HVT does not conflict with quantum theory regarding allowed values 
for physical quantities. However, local hidden variable theory predicts certain 
inequalities for the mean values of products of physical quantities for the two 
sub-systems. 

The form given by Clauser et al [82] for Bell’s inequality is 

l^l < 2 (70) 


where 

S = {Al X Hi) j^jjy -\- {Al X B 2 )]^uy + {A 2 X Hi )^^y - {A 2 X H 2 ) j^jjy (71) 

The minus sign can actually be attached to any one of the four terms. 

Following the proof of the Bell inequalities in m we have 

{A 2 X Bi)^^y - {A 2 X B2)^hv = J ((^2(0) {Bl{0) - (712(0) (^2(0)) 

= jdmiOUMO) (i?i(0)(i±(7ii(0) (^2(0)) 

- [dmionMO) (i?2(0)(i±(7ii(0) {Blim 

(72) 

Now all the quantities {Ai{^)), {Bj{^)) are bounded by -1-1 or —1, so the ex¬ 
pressions (1 ± (Hi(0) (^ 2 ( 0 )) and (1 ± (^i(O) (^i(O)) are never negative. 
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Taking the modulus of the left side leads to an equality 


1(^2 X - {A 2 X B2)i^fjy\ 

< Id^p{0i\{A2m\ (^ 2 ( 0 )) 

+ 1 d^p{oi\{A2m\ mm a ± {Aim {sm) 

< j dmiO{i±{Am) {B2m)+ J dmiOi^^iAm {Siim 
= 2±{Jd^p{o {Aim {B2m+Jdmio {a^)) {Bi{0)) 

= 2±{{Ai X B 2 )[^uy + {Ai X Bi)j^jjy) (73) 

where we have used the results that |(A 2 (^))| , |(7?i(C))l ^.nd |(-B 2 ( 0 )I Isss 

than unity and that J d^iO = 1 - Hence since |(Hi x B 2 ) j^jjy + {Ax x Bi)^^y\ = 

+ ((Hi X B 2 )]^uy + {Ai X Bi)^^y) or -((Hi x B 2 )j^^y + (Hi x Bi)^jjy) we 
have 

|(H2 X - (H 2 X 52 ) 1 ^ 1 ^ 1.1 ±|(Hi X B2)lhv + (^1 X - 81 ) 1 ^ 111.1 ^ 2 (74) 

But since jX —F| < |X| + |F| we see that from the + version of the last 
inequality that 

|(H 2 X Bi)^fjy - (H 2 X B2)i^hv + (^1 X B 2 ) HVT + (Hix 

Bi) hvt I < 2 (75) 

This is a Bell inequality. Interchanging H 2 O Hi and repeating the derivation 
gives |(Hi X Bi) jjyj, — (Hi x B 2 ) j^fjy + (H 2 x B 2 ) j^^jy + (H 2 x Bi) j^^y\ < 2, 
which is another Bell inequality. Interchanging Bi B 2 and repeating the 

derivation gives |(H 2 x B 2 )~ ("^2 x Bi)^j^y + (Hi x Bi)^^y + (Hi x B 2 )^hv\ ^ 
2 , and interchanging H 2 O Hi and Bi ^ B 2 and repeating the derivation gives 
|(Hi X B 2 )[^uy - (Hi X Bi)^jjy + (H 2 X Bi)^^y + (H 2 X B 2 )i^uy\ < 2. Thus 
the minus sign can be attached to any one of the four terms. An example of an 
entangled state that violates the Bell inequality is given in SubSection 12.7.31 

2.7.2 Non-Entangled State Result 

It can be shown that the Bell inequalities also always occur for non-entangled 
states (see Section 7.3 of the book by Vedral [13]). For Bell’s inequalities we 
consider Hermitian operators Ai and Bj for subsystems H, B respectively, for 
which there are two eigenvalues +1 and — 1 , where examples of the operators are 
given by the components H^ = • a a and Bj = bj - as of Pauli spin operators 

tT .1 and ds along directions with unit vectors and bj. The corresponding 
quantum theory quantity for the Bell inequality is 

S = 5(Hi ® Bi) + E{Ai 0 B 2 ) + E{A2 0 Bi) — E{A2 0 B 2 ) (f 6 ) 
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where in quantum theory the mean value is given by E{Ai®Bj) = l^Ai ® Bjj = 

Tr(p Ai ® Bj). For the general bipartite non-entangled state given by [3] it is 
easy to show that 


5=E ((^.); + ml («■ - m,) (”) 

R 


where = Tr{Aip^) and = Tr(Bjp^) are the expectation values 

of Ai and Bj for the sub-systems A, B in states 'p^ and respectively. Now 
Ai) and ( Bi ) must lie in the range —1 to -bl, so that ( Bi ± ^2 > must 


R \ / R 

each lie in the range —2 to +2- Hence 


R 


Bn 


< 

'^Pr I 

[\{m 


R 


< 

T.Pr I 

(|(Bi 


R 


< 

2 


i obtain | (Hi 

+ B 2 



\s' 

B 2 

>d = 1 



' R ' 

- \ -B 

/ ^ 

; Bi = 

-(B2 


/ R 

\ 

1 = 

0. 



Bi + Bn 


R 


Hi — Bn 


R 


B\ — Bn 


I = 2 requires (b^ = (b 2 

^ ^ \B ^ 

B 2 ) I = 0 , or to obtain | 
' R 

B ! ^ ^ 

= ±1 and then \IBi+B 2 
R \ 


(78) 


= ±1 and then 


B^ — Bn 


= 2 


= \{B 


R 


2.7.3 Bell Inequality Violation and Entanglement 

It follows that for a general two mode non-entangled state 151 cannot violate the 
Bell inequality upper bound of 2. Thus, the violation of Bell inequalities proves 
that the quantum state must be entangled for the sub-systems involved, so Bell 
inequality violations are a test of entanglement. For entangled states such as 
the Bell state |'I>_) (see [14], Section 2.5) written in terms of eigenstates of 
and 

|4-_) = ^(I+1)a ® |-l)s - |-1 )a ® I+1)b) (79) 

we find that 

£’(_^ • (g) (80) 

The Bell inequality (1751) can be violated for the choice where 61 and 62 are or¬ 
thogonal and oi, 02 are parallel to &i + 62 , — i >2 respectively (see [H, Section 

5.1). Furthermore, such a quantum state cannot be described via a hidden 
variable theory, since Bell inequalities are always satisfied using a hidden vari¬ 
able theory. Experiments have been carried out in optical systems providing 
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strong evidence for the existence of quantum states that violate Bell inequali¬ 
ties with only a few loopholes remaining(see [30], [15] and [6] for references to 
experiments). Such violation of Bell inequalities is clearly a non-classical fea¬ 
ture, since the experiments rule out all local hidden variable theories. As Bell 
inequalities do not occur for separable states, the experimental observation of 
a Bell inequality indicates the presence of an entangled state. These violations 
are not without applications, since such Bell entangled states can be useful in 
device-independent quantum key distribution m, m, ca. 

2.8 Non-local Correlations 

Another feature of entangled states is that they are associated with strong cor¬ 
relations for observables associated with localised sub-systems that are well- 
separated, a particular example being EPR correlations between non-commuting 
observables. Entangled states can exhibit this particular non-classical feature, 
which again cannot be accounted for via a hidden variable theory. 

2.8.1 Local Hidden Variable Theory 

Consider two operators id a and Ob associated with sub-systems A and B. These 
would be Hermitian if observables are involved, but for generality this is not 
required. In a local hidden variable theory these would be associated with func¬ 
tions Oc(0 (C = A,B) of the local hidden variables with the Hermitean 
adjoints O^ being associated with the complex conjugates Op(^). In local hid¬ 
den variable theory correlation functions are given by the following mean values 

(0*^0^ X O^Ob)bbV = J d^P{O^AmAif)il*BmBiO (81) 

which then can be shown to satisfy the following correlation inequality 

I (O^ X Hb)bhv I — {^a^a X O^Ob) bhv (82) 

This result is based on the inequality 

I dfPiOCiO > {jd^PiOV^y (83) 

for real, positive functions C(^),P(^) and where J d^P{f) = 1, and which is 

proved in Appendix[7j In the present case we have C)^) = r2J^(^)nyi(f) ilBiO^B{f), 
which is real, positive. A violation of the inequality in Eq. (1821) is an indication 
of strong correlation between sub-systems A and B. 
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2.8.2 Non-Entangled State Result 

It can be shown that the correlation inequalities are always satisfied for non- 
entangled states. In quantum theory the correlation functions are given by 
® = Tr{pQ^^ ® Qb) and (8) (g) 

For a non-entangled state of sub-systems A and B we have 

R 

A B 

= ^^Pr (84) 


Now 

\(n\<g)nB)\<Y.PR\{n^y\\{nB)l\ (85) 

R 

since the modulus of a sum is always less than the sum of the moduli. Using 

we obtain the Schwarz 

inequality - which is true for all states - ~ ^ ~ 

I ( UM p, and hence 


I 0 Ub ) I < Pr 

R 



B 


B 

R 


( 86 ) 


Next we use the inequality 



(87) 


for real, positive functions Cb,Pr and where J2 rPr ~ This inequality, 
which was used in the paper by Hillery et al [S3] , is proved in Appendix |7| In 

the present case we have Cr = so that 

I 0 1^ < Pr ^U^Ur^ =^U^Ua0U^Ur^ (88) 

R R R 

Thus for a non-entangled state we obtain the correlation inequality 

I 0 Ur^ P = I (Qa 0 P < (n'^A^A 0 ^b^b) (89) 

where the general result ^r2^0UR^ = ^U^0U^^ has been used. Thus 
non-entangled states have correlation functions that are consistent with hidden 
variable theory. 
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2.8.3 Correlation Violation and Entanglement 

Hence if it is found that the correlation inequality is violated | P = 

I (tIa 0 ® then the state must be entangled, so the 

correlation inequality violation is also a sufficiency test for entanglement. 
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3 Identical Particles and Entanglement 

We now take into account the situation where systems of identical particles are 
involved. This requires us to give special consideration to the requirement that 
quantum states in such cases must conform to the symmetrisation principle m- 
Further, entanglement is defined as a property that involves systems with two (or 
more) sub-systems, and the definition requires the specification of sub-systems 
that are distinguishable from each other and on which measurements can be 
made. In addition, the sub-systems must be able to exist as separate systems 
which can in principle be prepared in quantum states for that sub-system alone. 
This feature is vital to the definition of separable (or non-entangled) states on 
which the defintion of entangled states is based. These key requirements that 
the sub-systems must be distinguishable, susceptible to measurements and can 
exist in separate quantum states are necessary for the concept of entanglement 
to make physical sense, and will have important consequences for the choice 
of sub-systems when identical particles are involved. These three key logical 
requirement for sub-systems rule out considering labelled identical particles as 
sub-systems and lead to the conclusion that sub-systems must be modes or sets 
of modes. 


3.1 Symmetrisation Principle 

Whether entangled or not the quantum states for systems of identical parti¬ 
cles must conform to the symmetrisation principle, whereby for mixed states 
the overall density operator has to be invariant under permutation operators, 
or if pure states are involved, the state vector is either unchanged (bosons) or 
changes sign (fermions) if the permutation operator is odd. Either a first quan¬ 
tisation approach in which the basis states are written as symmetrised products 
of single particle states occupied by labeled identical particles can be used, or a 
second quantisation approach where the basis states are products of Fock states 
for all single particle states (modes), each Fock state specifying the number of 
identical particles occupying the particular mode. In first quantisation the sym¬ 
metrisation process removes any distinction between identical particles, whereas 
in second quantisation only mode creation operators are involved, and these do 
not involve labeled particles. Symmetrization is built into the definition of the 
Fock states. The two approaches are equivalent, but as we will see the second 
quantisation approach is more suited to identifying sub-systems and defining 
entanglement in systems of identical particles. 

It is useful to clarify some of the issues involved by considering a simple 
example. Since density operators can always be expressed in a diagonal form 
involving their orthonormal eigenstates 1$) with real, positive eigenvalues P{^) 
as p = ^ P(<I>) |<I>) ($1 and each |<I>) can always be written as a linear combina- 

<j> 

tion of basis vectors |'I'), we will focus on these basis vectors and their forms in 
both first and second quantisation. We consider a system with N = 2 particles, 
which may be identical and are labeled 1 and 2, or they may be distinguishable 
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and labeled a and /3. In each case a particle has a choice of two modes which 
it may occupy. Thus there are two distinct single particle states (modes) des¬ 
ignated as |A) and \B) in the identical particle case, and four distinct single 
particle states (modes) designated as \Aa), \Ba) and \Ap), jS/g) in the distin¬ 
guishable particle case for particles a and P respectively. The notation in first 
quantisation is that \C{i)) refers to a vector in which particle i is in mode \C). 
The notation in second quantisation is that |n)(^ refers to a vector where there 
are n particles in mode \C). 

For the case of the identical particles we consider basis states for two bosons 
or for two fermions, which are written in terms of first quantization as 

= ■^{\A{l))®\B{2)) + \B{l))®\Am (90) 

= ■^{\A{l))®\B{2))-\B{l))®\Am (91) 

and clearly satisfy the symmetrization principle. In second quantization the 
basis state in both the fermion and boson cases is 

\^)boson, fermion = \^)A \^)b (^ 2 ) 

In both first and second quantisation this basis state involves one identical 
particle in mode \A) and the other in mode \B). 

These examples highlight two possibilities for specifying sub-systems for sys¬ 
tems of identical particles. The two possibilities have differing consequences in 
terms of whether specific pure states are regarded as separable or entangled 
in terms of the general form in Eq. o for separable pure states, depending on 
whether the first or second quantisation approach is used. The first option is 
to regard the labeled identical particles as sub-systems - in which case using 
first quantisation the boson or fermion basis states in Eas. (l90l) and (IM]) would 
be regarded as entangled states of the two sub-systems consisting of particle 
1 and particle 2 m, m, m- This is a more mathematical approach, and 
suffers from the feature that the sub-systems are not distinguishable and mea¬ 
surements cannot be made on specifically labelled identical particles. In the 
case of identical particles the option of regarding labeled identical particles as 
the sub-systems leads to the concept of entanglement due to symmetrisation. 
In the textbook by Peres m, see ppl26-128) it is stated that ’’two particles of 
the same type are always entangled”. Peres obviously considers such entangle¬ 
ment is a result of symmetrization. The second option would be to regard the 
modes or single particle states as sub-systems [29) - in which case using second 
quantisation the basis state for both fermions or bosons in Ea. dMl) would be 
regarded as a separable state of two sub-systems consisting of modes |^) and 
\B). This is a more physically based approach, and has the advantage that 
the sub-systems are distinguishable and measurements can be made on specific 
modes. Noting that in the example the same quantum state is involved with 
one identical particle in mode |A) and the other in mode \B), the different cat¬ 
egorisation is disconcerting. It indicates that a choice must be made in regard 


52 


to defining sub-systems when identical particles are involved (see SubSection 

imi). 

Now consider the case where the particles are distinguishable. Each distin¬ 
guishable particle a, j5 has its own unique set of modes Ba, Ajs, Bp . 
There are two cases in which one particle a occupies mode \Aa) or \Ba) and 
the other particle /3 occupies mode \Ap) or \Bp). Basis states analogous to the 
previous ones are given in first quantization as 

\'^)di8t = \^c.{a)) (^\Bp{f)) or \'i’)dist = \Bc{a)) (^\Ap{l3)) (93) 

The somewhat surplus particle labels (a) and (/3) have been added for com¬ 
parison with o and m- The states (US]) are not required to satisfy the 
symmetrization principle since the particles are not identical. Each may be 
either a boson or a fermion. In second quantisation the basis states are 

= (| 1 )ac ® | 0 )bc) ® ® 11 ) 5 ^) 

or 

= (|0) Ac ® ® ® 1^)53) ( 94 ) 

In both first and second quantisation, the first case corresponds to particle a 
being in mode | A„) and particle /3 being in mode \Bp) with the other two modes 
empty, and the second case corresponds to particle a being in mode \Ba) and 
particle /3 being in mode \Ap) with the other two modes empty. 

These examples also highlight two possibilities for specifying sub-systems 
for systems of distinguishable particles. In this case the two possibilities have 
similar consequences in terms of whether specific pure states are regarded as 
separable or entangled, based on the general form in Eq. m for separable pure 
states, irrespective of whether the first or second quantisation approach is used. 
Here the first option is to regard the labeled distinguishable particles as sub¬ 
systems - in which case using first quantisation the boson or fermion basis 
states in Eas. (|M)l would be regarded as separable states of the two sub-systems 
consisting of particle a and particle /3. The second option would be to regard 
the modes or single particle states as sub-systems - in which case using second 
quantisation the basis state for both fermions or bosons in Ea. dMl) would be 
regarded as a separable state of four sub-systems consisting of modes |Hq), \Ba,) 
and \Ap)^ \Bp). Both expressions refer to the same quantum state, and the same 
result regarding separability is obtained in both first and second quantisation, 
even though the number of sub-systems differ. It indicates that either option 
may be chosen in regard to defining sub-systems when distinguishable particles 
are involved. However, it is simpler if the same option - particles or modes as 
sub-systems - is made for treating either identical or distinguishable particle 
systems and we will adopt this approach. 

To highlight the distinction between the identical and distinguishable par¬ 
ticles situation, we note that for the two distinguishable particle case treated 
previously we can also form entangled states from the basis states dMj) or dMj) 

IT) = ^(|H„(a)) 0 \Bpm ± |i?„(a)) 0 \Apm) (95) 
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which are similar in mathematical form to and m when written in first 
quantisation, and which are given by 

~ ® ® (|0)y4p ® ^ ® ® ® |0)_B^)) 

(96) 

when written in second quantisation. However, in this case both the first and 
second quantisation forms are clearly cases of entangled states. Whether they 
are regarded as entangled states of two sub-systems consisting of particle a and 
particle /3 (first option) or entangled states of the four sub-systems consisting of 
modes |Aq), \Ba) and \Ap), \Bp) (second option) depends on whether particle 
or modes are chosen as sub-systems. 

Note however that not all basis states result in separable/entangled dis¬ 
tinctions even in the case of identical particles. For the same two mode, two 
particle case as considered previously for bosons the basis vectors |H(l))(g)|H(2)) 
or \B{1)) 0 \B{2)) (first quantisation) or equivalently |2)^ 0 lO)^ or |0)^ 0 \2)g 
(second quantisation) would be regarded as separable states irrespective of 
whether particle or modes were chosen as the sub-systems. Entangled states 
such as (1^(1)) 0 1^(2)) ± |i?(l)) 0 \B(2)))/y/2 (first quantisation) or equiva- 
lently (|2)^ 0 |0)g ± |0)^ 0 \2)g)/\/2 (second quantisation) can also be formed 
from the two doubly occupied basis states. There are no analogous states for 
fermions due to the Pauli principle. 

It is worth noting that these examples illustrate the general point that just 
the mathematical form of the state vector or the density operator alone is not 
enough to determine whether a separable or an entangled state is involved. The 
meaning of the factors involved also has to be taken into account. Failure to 
realise this may lead to states being regarded as separable when they are not 
(see Subsection 13.41 for further examples). 

In the above discussion the symmetrisation principle was complied with both 
in the first and second quantisation treatments. It should be noted however 
that some authors disregard the symmetrisation principle. In describing Bose- 
Einstein condensates (162], [94]) consider states of the form 

p = ^PrP^0P^0P^0 ... (97) 

R 

as defining non-entangled states, where is a density operator for particle i. 
However such a state would not in general be allowed, since the symmetrisation 
principle would be violated unless the 'p'^ were related. For example, consider 
the state for two identical bosonic atoms given by 

p = 0 ^ + Pgr, 0 0 (98) 

and apply the permutation P = P(1 •«->■ 2). The invariance of p in general re- 

^ ^ ^ ^ ^ ^1 ^2 ''1 ''2 
quires a = f and 0 = ij, giving p = P„a 0ct + PeO 00. This is a statistical 

mixture of two states, one with both atoms in state a, the other with both atoms 
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in state 9. Thus only special cases of (IMl) are compatible with the symmetrisa- 
tion principle. Of course if the atoms were all different (atom 1 a Rb®^ atom, 
atom 2 a Na^® atom, ..) then the expression (155)) would be a valid non-entangled 
state, but there the atomic sub-systems are distinguishable and symmetrisation 
is not required. Such authors are really ignoring the symmetrisation principle, 
and in addition are treating the individual identical particles in the BEC as sep¬ 
arate sub-systems - a viewpoint we have described previously and will discuss 
further in the next SubSubSection. For the present we just point out that valid 
quantum states must comply with the symmetrisation principle. 


3.1.1 Sub-Systems - Particles or Modes ? 

As highlighted in the previous SubSection 13.11 when the quantum system in¬ 
volves identical particles the very definition of entanglement itself requires spe¬ 
cial care in regard to identifying legitimate sub-systems. There is a long¬ 
standing debate on the issue, with at present two schools of thought - see reviews 
such as [30] or [M] . As explained in the previous SubSection, the first approach 
is to identify mathematically labelled individual identical particles as the sub¬ 
systems m. m. m, izni. Sub-systems may of course also be sets of such 
individually labeled particles. This approach leads to the conclusion that sym¬ 
metrisation creates entanglement of identical particles. The second approach 
is to identify single particle states or modes that the identical particles may 
occupy as the sub-systems [29]. The sub-systems may of course also be sets 
of distinguishable modes. This approach leads to the conclusion that it is m- 
teraetion processes between modes that creates entanglement of distinguishable 
modes. 

The approach based on particle entanglement is still being used |70j . As 
explained in SubSection 13.II this is not the same as mode entanglement so tests 
and measures for particle entanglement will differ from those for mode entan¬ 
glement. A further discussion about the distinction is given in |35] . In a recent 
paper Killoran et al [5S] considered original states such as (|aO(l)) ® |al(2)) ± 
|a0(2)) 0 |al(l)))/-\/2 involving two modes aO and al - which were considered 
(based on first quantisation) as an entangled state for two sub-systems consist¬ 
ing of particles 1 and 2, but would be considered (in second quantisation) as 
a separable state |1 )qq 0 |l)ai for two sub-systems consisting of modes aO and 
al. In addition there were two modes 50 and 51 which are intially unoccupied. 
The particles may be bosons or fermions. They envisaged converting such an 
input state using interferometer processes which couple A modes aO and al 
to previously unoccupied B modes 50 and 51, into an output state - which is 
different. Projective measurements would then be made on the output state, 
based on having known numbers of particles in each of the A mode pairs aO and 
al and in the B mode pairs 50 and 51. The projected state with one particle in 
the A modes and one particle in the B modes would be of the form (in second 
quantisation) (|l)_^o 0 10)^^ (g) |0)jo (g) |l)f,i ± |0)^o (g) |l)^i (g) |l)j,o (g) |0)jJ/\/2 , 
which is a bipartite entangled state for the two pairs of modes A and B and is 
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mathematically of the same form as the first quantisation form for the original 
A modes state considered as an example of particle entanglement if the corre¬ 
spondences |aO(l)) ^ |l)^o (g) |0)^i, |al(2)) ^ |0)^o 0 |a0(2)) ^ |l)^o 10)^^ 

and |al(l)) —>■ |0)^g |l)ai ^^e made. Even the minus sign is obtained in the the 
fermion case. Details are given in Appendix [^Killoran et al stated that this 
represented a way of extracting the original symmetrization generated entan¬ 
glement. However, another point of view is that the two mode interferometer 
process created an entangled state from a non-entangled state, and as the final 
measurements are still based on entanglement of modes it is hard to justify the 
claim that entanglement due to symmetrization exists as a directly observable 
basic feature in composite quantum systems - though the mapping identified 
in |95j is.mathematically correct. Furthermore, all quantum states for identical 
particles are required to be symmetrized, so if symmetrization causes entan¬ 
glement it differs from the numerous other controllable processes that produce 
entanglement by coupling the sub-systems. Since the idea of extracting entan¬ 
glement due to symmetrisation is of current interest, a fuller discussion of the 
approach by Killoran et al [SS] is set out in Appendix O 

However, it is generally recognised that sub-systems consisting of individ¬ 
ually labeled identical particles are not amenable to measurements. What is 
distinguishable for systems of identical bosons or fermions is not the individual 
particles themselves - which do not carry labels, boson 1, boson 2, etc. - but 
the single particle states or modes that the bosons may occupy. For bosonic or 
fermionic atoms with several hyperfine components, each component will have 
its own set of modes. For photons the modes may be specihed via wave vectors 
and polarisations. Although the quantum pure states can be specified via sym¬ 
metrized products of single particle states occupied by specific particles using 
a first guantization approach, it is more suitable to use second guantization. 
Here, a basis set for the quantum states of such sub-systems are the Fock states 
\n)A (n = 0,1,2, ...) etc, which specify the number of identical particles occu¬ 
pying the mode A, etc., so in this approach the mode is the sub-system and the 
Fock states give different quantum states for this sub-system. Symmetrization 
is built into the definition of the Fock states, so the symmetrisation principle is 
automatically adhered to. If the atoms were fermions rather than bosons the 
Pauli exclusion principle would of course restrict n = 0, 1 only. In this second 
quantization approach situations with differing numbers of identical particles 
are recognised as being different states of a system consisting of a set of modes, 
not different systems as would be the case in hrst quantisation. The overall 
system will be associated with quantum states represented in the theory by 
density operators and state vectors in Fock space., which includes states with 
total numbers of identical particles ranging from zero in the vacuum state right 
up to infinity. Finally, the artificial concept of entanglement due to symmetrisa¬ 
tion is replaced by the physically realistic concept of entanglement due to mode 
coupling. 

The point of view in which the possible sub-systems A, B, etc are modes (or 
sets of modes) rather than particles has been adopted by several authors (|27j. 
EH], [29]), [13], [44] , [45] and will be the approach used here - as in |4] . To empha- 
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sise - what are or are not entangled in the present treatment involving systems 
of identical particles are distinguishable modes not labelled-indistinguishable- 
particles. Overall, the system is a collection of modes, not particles. Particles 
are associated with mode occupancies, and therefore related to specifying the 
quantum states of the system, rather than the system itself. In terms of this ap¬ 
proach, for non-interacting identical particles at zero temperature, the ground 
states for Bose-Einstein condensates and Fermi gases trapped in a harmonic 
potential provide examples of non-entangled states for bosonic and fermionic 
atoms respectively, when the sub-systems are chosen as the harmonic oscillator 
(HO) modes. In the bosonic case all the bosons occupy the lowest energy HO 
state, in the fermionic case one fermion occupies each HO state from the lowest 
up to a high energy state (the Fermi energy) until all the fermions are accom¬ 
modated. On the other hand, if one particle position states spatially localised 
in two different regions are chosen as two sub-systems, then the same zero tem¬ 
perature state for the identical particle system is spatially entangled, as pointed 
out by Goold et al [33] . Note that in this approach states where there is only a 
single atom may still be entangled states - for example with two spatial modes 
A, B the states which are a quantum superposition of the atom in each of these 
modes, such as the Bell state (|1)^ |0)^ -I- |0)^ \l)g)/\/2 are entangled states. 
For entangled states associated with the EPR paradox or for quantum telepor¬ 
tation, the mode functions may be localised in well-separated spatial regions - 
spooky action at a distance - but spatially overlapping mode functions apply in 
other situations. This distinction is important in discussions of quantum non¬ 
locality. Atoms in states with overall spin zero only have one internal state, 
but two mode systems can be created for their spatial motion using double-well 
trap potentials. If the wells are separated then two spatially separated modes 
can be created for studies of quantum non-locality. On the other hand atoms 
with spin 1/2 have two internal states, which constitute a two mode system. 
However these two modes may be associated with the same or overlapping spa¬ 
tial wave functions, in which case studies of quantum non-locality are precuded. 
These latter situation can however still lead to what is referred to as intrasys¬ 
tem entanglement [50] • Furthermore, as well as being distinguishable the modes 
can act as separate systems, with other modes being ignored. For interacting 
bosonic atoms this is much harder to accomplish experimentally than for the 
case of photons, where the relatively slow processes in which photons are de¬ 
stroyed in one EM field mode and created in another may require the presence 
of atoms as intermediaries. Two bosonic atoms in one mode may collide and 
rapidly disappear into other modes. However, atomic boson interactions can be 
made very small via Feshbach resonance methods. Near absolute zero the basic 
physics of a BEG in a single trap potential is describable via a one mode theory. 
Hence with A, B, .. signifying distinct modes, the general non-entangled state 
is given in Eq. (|3|) though the present paper mainly involves only two modes. 

As pointed out in SubSection l3.ll in the case of systems consisting entirely of 
distinguishable particles the sub-systems may still be regarded as sets of modes, 
namely those single particle states associated with the particular distinguish¬ 
able particle. In this case the particle descriptor (He atom, Na atom, ..) is 
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synonomous with its collection of modes. Here all the sub-system states are one 
particle states. 


3.1.2 Multi-Mode Sub-Systems 

As well as the simple case where the sub-systems are all individual modes, the 
concept of entanglement may be extended to situations where the sub-systems 
are sets of modes, rather than individual modes, In this case entanglement 
or non-entanglement will be of these distinct sets of modes. Such a case in 
considered in SubSection 4.3 of paper II, where pairs of modes associated with 
distinct lattice sites are considered as the sub-systems. Another example is 
treated in He et al [30], which involves a double well potential with each well 
associated with two bosonic modes, these pairs of modes being the two sub¬ 
systems. Entanglement criteria for the mode pairs based on local spin operators 
associated with each potential well are considered (sse SubSections 4.2 and 5.3 
of paper H). A further example is treated by Heaney et al [37] , again involving 
four modes associated with a double well potential. As in the previous example, 
each mode pair is associated with the same well in the potential, but here a Bell 
entanglement test was obtained for pairs of modes in the different wells. The 
concept of entanglement of sets of modes is a straightforward extension of the 
basic concept of entanglement of individual modes. 

3.2 Super-Selection Rule 

As well as the symmetrisation principle there is a further requirement that quan¬ 
tum states of systems of identical particles must satisfy - these are known as 
super-selection rules. These rules restrict the allowed quantum states of such 
systems to those in which the coherences between states with differing numbers 
of particles are zero. This applies at the global level for the overall quantum 
state, but also - as will be discussed in a later sub-section - to the sub-system 
states involved in the definition of separable or non-entangled states. The just- 
hcation of the SSR at both the global and local level will be considered both in 
terms of simple physics arguements and in terms of reference frames. Examples 
of SSR and non-SSR compliant states will be given, both for overall states and 
for separable states. The validity of the SSR for the case of massive bosons 
or fermions is generally accepted, but in the case of photons there is doubt re¬ 
garding their applicability -as will be discussed below. As pointed out in the 
Introduction, in the case of systems consisting entirely of single distinguish¬ 
able particles the sub-systems may still be regarded as sets of modes, namely 
those single particle states associated with the particular distinguishable parti¬ 
cle. Here all the sub-system states are one particle states and the overall system 
is an N particle state, so the local and global particle number super-selection 
rules, though true are irrelevant. 
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3.2.1 Global Particle Number SSR 


The question of what quantum states - entangled or not - are possible in the 
non-relativistic quantum physics of a system of identical bosonic particles - such 
as bosonic atoms or photons - has been the subject of much discussion. Whether 
entangled or not it is generally accepted that there is a super-selection rule that 
prohibits quantum superposition states of the form 

OO OO OO OO 

|d>) = J:Cn\N) p= J: \N) (iV|+ E E il-SN,M)CN C*M |A^) (M| 

N=0 N=0 N=0 M=0 

(99) 

being quantum states when they involve Fock states \N) with differing total 
numbers N of particles. The density operator for such a state would involve 
coherences between states with differing N. Although such superpositions - 
such as the Glauber coherent state |a), where Cn = exp(—|ap/2) /'/N\ - do 

have a useful mathematical role, they do not represent actual quantum states 
according to the super-selection rule. The papers by Sanders et al m and 
Cable et al [35] are examples of applying the SSR for optical fields, but also 
using the mathematical features of coherent states to treat phenomena such as 
interference between independent lasers. The super-selection rule indicates that 
the most general quantum state for a system of identical bosonic particles can 
only be of the form 


P = E ’^P'S'N (I'&Tv) (‘fArl) 

N=o $ 

\^n) = ( 100 ) 

i 

where \^n) is a quantum superposition of states \Ni) each of which involves 
exactly N particles, and where different states with the same N are designated 
as \Ni). This state p is a statistical mixture of states, each of which contains 
a specific number of particles. Such a SSR is referred to as a global SSR, as it 
applies to the system as a whole. Mathematically, the global particle number 
SSR can be expressed as 

[iV,)3]=0 (101) 

where N is the total number operator. 

3.2.2 Examples of Global Particle Number SSR Gompliant States 

Examples of a state vector |<i)vr) for an entangled pure state [55] and a density 
operator 'p for a non-entangled mixed [33] state for a two mode bosonic system, 
both of which are possible quantum states are 

N 

1$^) = j:CiN,k)\k)^(i^\N-k)s (102) 

N 

P = E P(fc) |fc)^(fc|^ (g) |N-fc)^(fV-(103) 

/c=0 


59 


The entangled pure state is a superposition of product states with k bosons 
in mode A and the remaining N — k bosons in mode B. Every term in the 
superposition is associated with the same total boson number N. The non- 
entangled mixed state is a statistical mixture of product states also with k 
bosons in mode A and the remaining N — k bosons in mode B. Every term in 
the statistical mixture is associated with the same total boson number TV. For 
the case of a two mode fermionic system the Pauli exclusion principle restricts 
the number of possible fermions to two, with at most one fermion in each mode. 
Expressions for a state with exactly N = 2 fermions are 

|$2) = |1)a®|1)b (104) 

P = |1)a (l|y4 ® |1)b (1|b (105) 

Neither state is entangled and both are the same pure state since p = |<i) 2 ) (‘1*21- 
Although the super-selection rules and symmetrisation principle also applies 
to fermions, as indicated in the Introduction this paper is focused on bosonic 
systems, and it will be assumed that the modes are bosonic unless indicated 
otherwise. 

Bell states m, [la for = 1 bosons provide important examples of entan¬ 
gled two mode pure quantum states that are compliant with the global particle 
number SSR. The modes are designated A, B and the Fock states are in general 
IbajBb). These Bell states may be written 



h^AB) 

= -^(IOa, 1b) - |1a,0b)) 



I^-ab) 

= -^(IOa, 1b) + |1a,0b)) 

(106) 

Neither of these states is separable. There are also two other two mode Bell 
states given by 



= ;^(|0a,0b) - |1a, 1b)) 



|4>1b> 

= ■^(|0a,0b) + |1a, 1b)) 

(107) 


These however are not compliant with the global particle number SSR. Linear 
combinations = IO^jOb) and (- = 

I 1a, 1b) are global particle number SSR compliant and also separable, corre¬ 
sponding to states with N = 0 and N = 2 bosons respectively. 

3.2.3 Super-Selection Rules and Conservation Laws 

It is important to realise that such super-selection rules [52] are different con¬ 
straints to those imposed by conservation laws, as emphasised by Bartlett et al 
|60j . For example, the conservation law on total particle number only leads to 
the requirement on the superposition state |$) that the \Cn\'^ are time indepen¬ 
dent, it does not require only one Cn being non-zero. They are however related. 
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as is discussed in Section IB.d.ll and Appendix [TUI where the super-selection rules 
based on particle number are related to invariances of the density operator un¬ 
der changes of phase reference frames. This involves considering groups of phase 
changing operators T[9a) = exp(iNaOa) when considering local particle number 
SSR for single modes in the context of separable states, or T{9) = exp(*iV0) 
when considering global particle number SSR in the context of multimode entan¬ 
gled states. Super-selection rules are broad in their scope, forbidding quantum 
superpositions of states of systems with differing charge, differing baryon num¬ 
ber and differing statistics. Thus a combined system of a hydrogen atom and 
a helium ion does not exist in quantum states that are linear combinations of 
hydrogen atom states and helium ion states - the super-selection rules on both 
charge and baryon number preclude such states. The basis quantum states for 
such a combined system would involve symmetrised tensor products of hydrogen 
atom and helium ion states, not linear combinations - symmetrisation being re¬ 
quired because the system contains two identical electrons. On the other hand, 
super-selection rules do not prohibit quantum superpositions of states of sys¬ 
tems with differing energy, angular or linear momenta - other physical quantities 
that may also be conserved. Thus in a hydrogen atom quantum superpositions 
of states with differing energy and angular momentum quantum numbers are 
allowed quantum states. 

However, conservation laws on total particle number (such as apply in the 
case of massive bosons) are relevant to showing that multi-mode states generated 
via total particle number conserving processes from an initial separable state 
will be global SSR compliant if the sub-systems in the initial state are local 
particle number SSR compliant, and will not be if the initial state involves a 
sub-system state that is not local particle number SSR compliant. For simplicity 
we consider two sub-systems A and B with the initial state 

= (108) 
R 

If U{t) is the evolution operator where )o(t) = U{t)'p{Q)U and the processes 
are number conserving then [N, U{t)] = 0. We then have 

[iV,p(t)] = U{t)[N,m]U{t)^ 

= U{t)Y,PR{[NA,pi](^PR + pi^[NB,Pn])uity (109) 

R 

Hence if and are local particle number SSR compliant, then [Na, 'pr] 
and [Nb, p^] are zero, showing that [N,'p{t)] = 0 so the state is global par¬ 
ticle number SSR compliant. On the other hand if [fV,p(t)] = 0 we see that 
[iV,p(0)] = J2rPr Pr]‘^Pr + Pr‘^ [Nb, p|]) = 0. By taking Tta and 

TrB of this result gives Y^rPr [Na, Pk] ^ J2rPr [Nb, pi] = 0. This shows 
that both of the reduced density operators PrT^r and Pr7>r nrust be lo¬ 
cal particle number SSR compliant, which amounts to requiring the sub-system 
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density operators to be local particle number SSR compliant. This situation 
applies even when there is coupling between modes, provided the interaction is 
number conserving - such as a coupling given by R = Xab^ + HC. Analogous 
results apply for systems of massive bosons if there are more than two modes 
involved, where again global SSR compliance involves the total particle number 
since even with interactions there is total number conservation. For example 
with three modes in coupled BECs, interactions of the form V = X{c)^a^b^ + HC 
in which two bosons are annihilated in mode C and one boson is created in each 
of modes A and B are consistent with total particle number conservation and 
lead to global SSR involving the total particle number. 

Although outside the focus of this paper, it is worth pointing out that some¬ 
what different considerations apply to photons. Single non-interacting modes, 
such as are discussed in the context of separable states do have a conservation 
law for the photon number in that mode. The applicability (or otherwise) of the 
local particle number SSR for the sub-system density operators in a separable 
state is discussed in Section 13.2.91 In the case of interacting photonic modes 
there may be no conservation law associated with total photon number and it 
may be thought that no global SSR would apply. However, other global SSR 
involving combinations of the mode photon numbers may still apply. As an 
example, we consider a three mode situation in a non-degenerate parametric 
amplifier, where the basic generation process involves one pump photon of fre¬ 
quency ujc = WA +wb being destroyed and one photon created in each of modes 
A and B. The interaction term is H = XaAb^ + HC. It is straight-forward to 
show the a total quanta number operator Ntot = Na + Nb + 2Nc commutes 
with the Hamiltonian. The situation is analogous to the atom-molecule system 
treated in Appendix[TTJ Thus Ntot is conserved and we can then consider a group 
of phase changing operators T{9) = exp{iNtotd) and show that there could be a 
global SSR for the three mode system, but now involving the total quanta num¬ 
ber Na + Nb + ‘2Nc. The pure state which is often used in a quantum treatment 
of the non-degenerate parametric amplifier |'h) ='^Cn\n)^(§i\n) ^ (Si\N — n)^ 


n 


is global SSR compliant in terms of the modified Ntot, since in every term 
Na + Nb + 2Nc = 2N and there are no coherences between terms with dif¬ 
ferent Ntot- For the non-degenerate parametric amplifier case an analogous 
treatment to that for number conserving processes shows that if 


^( 0 ) = '^PrPr<^Pr^Pr 


( 110 ) 


R 


then using [Ntot,U{t)] = 0 we have 
[Ntot,m] = u{t)[Ntot,pmu{t)^ 



( 111 ) 
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Hence if and p^ are local particle number SSR compliant, then [iVtoi, p(t)] = 

0 so the state is SSR compliant, but with global total quanta number Ntot- On 
the other hand if [7Vtot,p(t)] = 0 we find that p;r] =J2rPr[Nb,Pr] = 

J2rPr Wc, Pfl] = 0. This shows that each of the reduced density operators 
^rPrPr , ^rPr7>r and ^rPr'Pr must be local particle number SSR com¬ 
pliant, which amounts to requiring the sub-system density operators to be local 
particle number SSR compliant. 

3.2.4 Global SSR Compliant States and Quantum Correlation Func¬ 
tions 

We now prove a theorem concerning quantum correlation functions for bosonic 
systems with two modes A and B. 

Theorem. If a state is global particle number SSR compliant then all quan¬ 
tum correlation functions for which n + I ^ m + k must 

be zero. 

Proof: If the state is global particle number SSR compliant then if we choose 
a complete orthonormal set of Fock states |A^, a) with a = 1,2, ... ,dN listing 
states which are eigenstates of the total number operator N with eigenvalue N 
we can write the density operator in the form 

N a,P 

where since Tr'p = 1 we must have 

i-EE^<(!« (113) 

N a 

Now (al)"(a)’”(&l)*(&)^ \N a) must be a linear combination of Fock states with 
N replaced byA^-|-n-|-Z — m — /cso we can write 

(al)”(a)™(Sl)'(S)'" \N,a) \{N + n + I - m - k),-f) 

7 

(114) 

Hence 

<^(al)”(a)’"(Sl)'(&)'=^ =Tr [ E E \{N + n + I - m - k),j) {N, /3| 

y N a,p 7 

(115) 

But Tr{\{N -\- n + I — m — k) {N, (3\) = 0 unless n-\-l — m — k = 0. Hence 
^(al)”(a)’”(6l)'(6)'=^ =0 ifn + l^m + k (116) 

which is the required theorem. 
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3.2.5 Testing the Super-Selection Rules 

The last result for the general two mode quantum correlation function 


is relevant to the various experimental measurements that are discussed in the 
accompanying paper II. For example, as we will see (Sx\ is a combination 


of with n = 1, m = 0, / = 0, fc = 1 and n = 0,m = 


IJ = I,fc = 0, and (ASl) = 

p 


Si) - 


Sx ) would involve terms such as 

p 


with n = 2,m = 0,1 = 0,k = 2 and n = 0,m = 2,1 = 

All of these have 


2,k = 0, and n = 1, m = 1,1 = 1, fc = 1 from {^S-^ 
n + I = m + k, so they can be non-zero for globally SSR compliant states. The 
question then arises - what sort of quantity of the form 
could be used to see if the quantum state was not globally SSR compliant? The 
answer is seen in terms of two corollaries to the last theorem. 

Corollary 1. If we find that any of the quantum correlation functions ^(ai)"(a)’"(6i)*(6)* 
are non-zero when n + I ^ m + k then the state is not global particle number 
SSR compliant. 

This result indicates what type of measurement is needed to see if SSR 
non compliant states exist. Quantities of the type are 

measured for which n + l m + k. If we find any that are non-zero we can then 
conclude that we have found a state which is not global particle number SSR 
compliant. 

Corollary . Measurements of the QCF ^(a'i)”(a)’”(6'i)*(&)^^ when n + l = 
m + k cannot determine whether or not a state is global particle number SSR 
non-compliant. 

If the state is non-compliant then its density operator must contain a contri¬ 
bution which allows for non-zero coherences between Fock states with different 
N. We can therefore write the density operator as 


N a,P 


N^M a,P 


Tp^f \N, a) {M, P\ 


(117) 


where the second term is the SSR non-compliant contribution. 

A similar calculation to before for the situation when n -\-1 = m -\- k gives 


(at)"(a)-( 5 t)'( 6 )'=) = Tr | ^ ^ ^ m, I, fc) |IV , 7 ) (A^,/3| 

N OL,j3 7 


-l-Tr 

yN^M a,/3 7 

(118) 
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But Tr{\N , 7) (M , / 3 |) = 0 for iV ^ M, so the non-compliant contribution gives 
zero and as Tr{\N, 7) {N, (3\) = we end up with 


{a^r{ar{Vf(bf)= (119) 


N a,13 


which is entirely dependent on the contribution to the density operator that is 
globally SSR compliant. Measurements of this type with n + I = m + k would 
not respond to the presence of contribution to the density operator that is not 
globally SSR compliant. The BS measurements discussed in this paper are all 
of this type, so will not test the super-selection rule. 

Hence the conclusion is that a quantum correlation function of the form 
must be measured for cases where n + I ^ m + k and 
a non-zero measurement result must be found. If it is, then we would have 
demonstrated that the state is not globally SSR compliant. The simplest case 
would be to find a non-zero result for (a)^ oi (b^ . 

Similar considerations apply to local SSR compliance in the sub-system 
states. For sub-system a a QCF of the form must be measured 

for cases where n ^ m and a non-zero measurement result must be found. If it 
is, then we would have demonstrated that the state is not locally SSR compliant. 
The simplest case would be to find a non-zero result for (d)^. 
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3.2.6 SSR Justification and No Suitable Phase Reference 

There are two types of justification for applying the super-selection rules for sys¬ 
tems of identical particles. The first approach is based on simple considerations 
and will be outlined below in this subsection. The second approach [53], [50] . 
ini], [SD, [53], [51], [53], [55], [57], [31] is more sophisticated and involves linking 
the absence or presence of SSR to whether or not there is a suitable reference 
frame in terms of which the quantum state is described, and is outlined in the 
next subsection and Appendix HU] The key idea is that SSR are a consequence of 
considering the description of a quantum state by an external observer (Charlie) 
whose phase reference frame has an unknown phase difference from that of an 
observer ((Alice) more closely linked to the system being studied. Thus, whilst 
Alice’s description of the quantum state may violate the SSR, the description 
of the same quantum state by Charlie will not. In the main part of this paper 
the density operator p used to describe the various quantum states will be that 
of the external observer (Charlie). 

3.2.7 SSR Justication and Physics Considerations 

A number of straightforward reasons have been given in the Introduction for 
why it is appropriate to apply the superselection rule to exclude quantum su¬ 
perposition states of the form (133)) as quantum states for systems of identical 
particles, and these will now be considered in more detail. 

Firstly, no way is known for creating such states. The Hamiltonian for such a 
system commutes with the total boson number operator, resulting in the |CAr|^ 
remaining constant, so the quantum superposition state would need to have 
existed initially. In the simplest case of non-interacting bosonic atoms, the Fock 
states are also energy eigenstates, such Fock states involve total energies that 
differ by energies of order the rest mass energy mc^, so a coherent superposition 
of states with such widely differing energies would at least seem unlikely in 
a non-relativistic theory, though for massless photons this would not be an 
issue as the energy differences are of order the photon energy hw. The more 
important question is: Is there a non-relativistic quantum process could lead to 
the creation of such a state? Processes such as the dissociation of M diatomic 
molecules into up to 2M bosonic atoms under Hamiltonian evolution involve 
entangled atom-molecule states of the form 

M 

1$) = S c™ |M - (120) 

m—0 

but the reduced density operator for the bosonic atoms is 

M 

Patoms= E \Cm\^ (l^m) {2m\) (121) 

m—0 

which is a statistical mixture of states with differing atom numbers with no co¬ 
herence terms between such states. Such statistical mixtures are valid quantum 
states, corresponding to a lack of a priori knowledge of how many atoms have 
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been produced. To obtain a quantum superposition state for the atoms alone, 
the atom-molecule state vector would need to evolve at some time into the form 

M M 

1$) = J: Bm\M- ® E I2n),,_ (122) 

m—0 n—0 

where the separate atomic system is in the required quantum superposition 
state. However if such a state existed there would be terms with at least one 
non-zero product of coefficient involving product states \M — 0 

\2n)afom n ^ TO if the state |<i)) is not just in the entangled form p20D . 
However, the presence of such a term would mean that the conservation law 
involving the number of molecules plus two times the number of atoms was 
violated. This is impossible, so such an evolution is not allowed. 

Secondly, no way is known for measuring all the properties of such states, 
even if they existed. If a state such as dMI) did exist then the amplitudes Cn 
would oscillate with frequencies that differ by frequencies of order mc^/H (the 
Compton frequency, which is > 10^® Hz for massive bosons) even if boson-boson 
interactions were included To distinguish the phases of the Cat in order to verify 
the existence of the state, measurement operators would need to include terms 
that also oscillate at relativistic frequencies, and no such measurement operators 
are known. 

Thirdly, there is no need to invoke the existence of such states in order to un¬ 
derstand coherence and interference effects..It is sometimes thought that states 
involving quantum superpositions of number states are needed for discussing 
coherence and interference properties of BECs, and some papers describe the 
state via the Glauber coherent states. However, as Leggett m has pointed out 
(see also Bach et al [99], Dalton and Ghanbari [32]), a highly occupied number 
state for a single mode with N bosons has coherence properties of high order n, 
as long as n <C iV. The introduction of a Glauber coherent state is not required 
to account for coherence effects. Even the well-known presence of spatial inter¬ 
ference patterns produced when two independent BEGs are overlapped can be 
accounted for via treating the BEGs as Fock states. The interference pattern is 
built UD as a result of successive boson oosition measurements [TOO]. m. [98]. 

3.2.8 SSR Justification and Galilean Frames ? 

Finally, in addition to the previous reasons there is an arguement that has 
been proposed based on the requirement that the dynamical equations for such 
non-relativistic quantum systems should be invariant under a Galilean trans¬ 
formation which has been proposed [101] as a proof of the super-selection rule 
for atom number. This approach is linked to the reference frame based justifi¬ 
cation of SSR (see Appendix fTUl). However, whilst the paper shows that under 
a Galilean transformation - corresponding to describing the system from the 
point of view of an observer moving with a constant velocity v with respect 
to the original observer, and where the two observers have identical clocks - 
the terms in a superposition state with different numbers N of massive bosons 
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would oscillate like expi /fi, and may be expected if the same quan¬ 

tum state is described by a moving observer. This feature alone does not seem 
to require the super-selection rule, since here the moving observer’s reference 
frame has a well-defined velocity with respect to that attached to the system. 
However, the moving observer’s reference frame may actually have an unknown 
relative velocity, in which case a twirling operation resulting in the elimination 
of number state coherences could be involved (see Appendix 1101) . This will be 
not be considered further at this stage. 

On the other hand, an approach of this kind involving rotation symmetry 
would seem to rule out such states as quantum superpositions of a boson (spin 
0) and a fermion.(spin 1/2). Let such a state be prepared in the form {\F) + 
\B))/^/2. Consider an observer whose cartesian reference frame is X, Y, Z. This 
is a classical system that can be rotated in space. If the observer rotates with 
his frame through 27r about any axis they are then back in the same position, 
but the observer now sees the state as (— \F) + \B))/\f2. This state is apparently 
orthogonal to the one observed before the rotation, and this is paradoxical since 
the observer would be in the same position. Thus there is a super-selection rule 
excluding states such as (|F) -|- \B))l\f2. A similar argument based on the time 
reversal anti-unitary operator was given by Wick et al [52]. 

3.2.9 SSR and Photons 

Though this paper is focused on massive bosonic atoms the question is whether 
similar SSR also apply to the optical quantum EM field, which involve massless 
bosons - photons. Here the situation is not so clear and we therefore merely 
present the differing viewpoints in the current literature. Some of the same gen¬ 
eral reasons for applying the super-selection rule to systems of identical massive 
bosons also apply here, though the details differ, but others do not. The situa¬ 
tion depends also on whether optical or microwave photons are involved. The 
issue is whether for individual photon modes, states can be prepared that are 
not local particle number SSR compliant, and if so can the effects of the non- 
SSR compliant terms be observed and furthermore do we need to invoke the 
existence of non SSR compliant states to understand interference and coherence 
effects. As we will see, some SSR non compliant feature needs to be present in 
order to prepared allegedly non SSR compliant states. Another way of looking 
at the issue is to ask whether phase reference systems exist for photon modes. 
In addition, there is the issue for multi-mode situations whether states can be 
prepared, observed or are needed that are not global particle number SSR com¬ 
pliant. As we will see, SSR may now involve a modified total particle number 
involving combinations of the mode numbers diferent to the total particle num¬ 
ber - because it is these combinations that are conserved in mode interaction 
processes. In the approach adopted in this paper, the states in question are those 
described by so-called external observers - not hypothetical observers that are 
somehow attached to phase reference systems internal to the experiment (see 
below. Section IS3|). The SSR issue for systems involving massless photons is 
particularly important in regard to describing entanglement. As explained be- 
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low in Section [311 if separate sub-system states in photonic systems can be 
prepared with density operators that violate the local particle number SSR, 
then these so-called ’’separable but non-local” states [H] would be classihed as 
separable rather than as entangled states. Some of the tests for entanglement 
described in Paper II for systems of massive boson (such as spin squeezing in 
any spin component) would then no longer apply for photonic systems, though 
others (such as the Hillery spin variance test) which do not depend on SSR 
would still apply. 

We first consider the requirement of showing how a non SSR compliant 
states can be prepared for single modes. In the case of photons, Molmer m 
has argued that the quantum state for a single mode optical laser field oper¬ 
ating well above threshold is not a Glauber coherent state, and the density 
operator would be a statistical mixture of the form (llOOl) . with I^at) = \N) 
and P$Ar = exp(—TV) N /N\. Here the density operator is a statistical mixture 
of photon number states with a Poisson distribution, or equivalently a statisti¬ 
cal mixture of equal amplitude coherent states |a) with a = ''/n exp{i4>) and 
all phases (j) having equal probability. In either form, the quantum state is 
SSR compliant. In terms of possible processes for preparing states for single 
mode optical laser fields, this feature is confirmed in theories for single mode 
lasers involving atomic gain media energised via incoherent pumping processes 

- there is no well defined optical phase that is imposed on the process. The 
Scully-Lamb theory (see Mandel and Wolf |103) . p935) gives the above thresh¬ 
old steady state density operator for the laser mode can be written in the form 
of a statistical mixture of number states (somewhat broader than for a Poisson 
distribution), which again is an SSR compliant state with no well-defined optical 
phase. Further detailed discussion of laser light generation processes by Pegg 
and Jeffers m confirms this. An alternative approach is presented by Wiseman 
et al [104] , [105) , in which the optical laser is treated via a master equation, but 
where monitoring of the laser environment (difficult!) is required to determine 
whether certain pure state ensembles - such as those involving coherent states 

- are physically realisable. The conclusion reached is that for finite self energy 
the coherent state ensemble is not physically realisable, the closest ensemble 
being that involving squeezed states, though for zero self energy coherent state 
ensembles are obtained. On the other hand, microwave photons in single mode 
high Q cavities can be generated by oscillating electric currents having a well- 
defined phase. In this case, as shown in experiments on the Jaynes-Cummings 
model bv Rempe et al |106| and Brune et al m demonstrating collapses and 
revivals of Rydberg atom population differences, it is possible to create Glauber 
coherent states in microwave cavity modes, and the presence of these states are 
necessary to explain the collapse and revival effects. 

We next consider cases where interacting photon modes are involved. The 
two mode squeezed states generated for example in a non-degenerate paramet¬ 
ric amplifier are often written in the form ^ Cn \n)j^ \n )corresponding to the 

n 

basic generation process in which a pump photon of frequency ujq = <^A + <^b 
is destroyed and one photon is created in each of modes A and B. Such a state 
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is not even global SSR compliant, but is used in describing various quantum in¬ 
formation processes as well as describing two mode squeezing. However, whilst 
mathematically convenient for treating such applications this do not demon¬ 
strate that this two mode pure state has actually been created. This state vector 
is in fact based on a very simplified version of the process, in which the pump 
mode is treated classically. If it is treated quantum mechanically and there were 
N photons initially in mode C, the interaction term V = Xca^b^ + HC would 
result in a global SSR compliant state vector like ^ |iV — |n)^ |n)^, 

n 

but now involving a total quanta number Ntot = Na + Nb + 2Nc (see Sec¬ 
tion [3^2^ for details). The state describing modes A and B alone would be 
'Pab — S l^'nP 1?^)^ ® In)g{n\g. This state is global SSR compliant in 

n 

terms of Na,b = Na + Nb but is not the same as the pure state '^Cn\n)^\n)g 

n 

- which is not SSR compliant. Even more elaborate quantum treatments al¬ 
lowing for irreversible damping processes for all three modes (see for example 
McNeil et al |108) 1 that result in non SSR compliant steady state solutions, 
include assumptions such as the pump mode being coupled to a laser mode 
that is treated classically - and thus begging the question of whether non SSR 
compliant states were prepared, since the classical treatment of the laser mode 
is itself SSR non compliant. We are unaware of any situation where non SSR 
compliant states are claimed to have been created for optical photons, where the 
theoretical treatment of the preparation process has not assumed the presence 
of non SSR compliant states for some key sub-system involved - usually in an 
input pump mode. In Section 13.2.31 we considered a preparation process for 
the non-degenerate parametric amplifier involving conservation of Ntot starting 
from an initial separable state in which all the sub-system density operators are 
local particle number SSR compliant and show that this results in a quantum 
state that is global SSR compliant in terms of total quanta number Ntot- There 
must therefore be some non SSR compliant feature in the initial state (which 
could include pump modes) to produce global non SSR compliant states, so 
then the issue shifts back to how these non SSR compliant states are prepared 
in the first place. 

Second, there is the requirement of being able to measure the non SSR 
compliant terms. For the free quantum EM field there is a conservation law 
for the photon number in each mode, so in a pure state such as in Eq. (IMl) 
the ICatP would be time independent. However, for photons the Cn would 
oscillate with frequencies that only differ by non-relativistic frequencies of order 
u! rather than the Compton frequency that applies for massive bosons, so the 
arguement against being able to detect coherent states based on this frequency 
being so large that the oscillations cannot be followed do not necessarily apply. 
Clocks that can follow microwave oscillations are common-place, and even at 
optical frequencies the development of atomic clocks based on optical atomic 
transitions that may enable optical frequency oscilations to be observed. So 
this consideration does not rule out non SSR compliant states, even for optical 
photons. 
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Finally, we consider the requirement of non SSR compliant states being 
needed to explain interference, coherence effects etc. We need to distinguish 
the sitation where it is mathematically convenient to invoke SSR non compliant 
states to explain these effects from the situation where it is essential to do 
so. Thus it may be convenient to explain the presence of interference patterns 
in position measurements for bosons from two independent BECs by choosing 
Glauber coherent states to represent their states, but as pointed out in Section 
13.2.71 such interference patterns are accounted in terms of Fock states, together 
with quantum interference of probability amplitudes associated with bosons 
being taken from the two different sources (see Refs. [TUO] . [SI], ISH])- In 
fact, the more detailed feature that although the separation of the peaks is 
well defined the actual position of the peaks are random, is inconsistent with 
the Glauber coherent state description. We also point out below (see Section 
13.3.21 and Appendix ITT]) that the interpretation of Ramsey fringes in a proposed 
experiment to detect a coherent superposition of an atom and a molecule does 
not show that such a state was created or that the BEG involved had to be 
described by a Glauber coherent state. Many experiments in which coherence, 
interference effects are observed do not depend on SSR non compliant states 
being created. Optical interference and coherence effects can also be explained 
without invoking Glauber coherent states, as shown by Molmer [102] and in 
other papers such as m- In regard to two mode squeezing in the non-degenerate 
parametric amplifier described above, the observation of squeezing effects is 
often discussed in terms of SSR non compliant states of the form \n)j^ W)^ 

_ n 

(see Ref. |15j for example). However, as may be seen from the experimental 
paper of Ou et al [109] . the way in which two mode squeezing is observed in 
the non-degenerate parametric amplifier involves generating the pump field by 
frequency doubling from a lower frequency laser. That lower frequency laser is 
also used to provide the local oscillator fields for the homodyne measurements on 
modes A and B used to detect squeezing - these modes are coupled to the local 
oscillator fields using beam splitters. The original lower frequency laser acts as 
an internal phase reference for the overall experiment, as may be seen in Eig 2 of 
Ref [109] which involves the relative phase between the local oscillator and the 
squeezed input field. But as there is no external phase reference system involved, 
only the relative phases of the A, B and local oscillator modes are well defined, 
and not the overall phase as would be required for preparing non SSR compliant 
states. Again the convenient use of SSR non compliant states to understand 
experiments in which only internal phase references are involved does not show 
that SSR non compliant states are necessary to interpret the experiments. A 
similar arrangement occurs for the degenerate parametric amplifier experiment 
of Wu et al m, where Eig 2 clearly shows how the local oscillator field derives 
from the original Nd-YAG laser. 

Perhaps the best way to approach the question of whether SSR compliance 
is required for optical photons for example (see next sub-section and SubSection 
110.41 in Appendix 1101) involves the consideration of phase reference frames. The 
quantum state of a single mode laser may be described as a Glauber coherent 
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state by an observer (Alice) with one reference frame, but would be described 
as a statistical mixture of photon number states by another observer (Charlie) 
with a different reference frame whose phase reference is completely unrelated 
to the previous one. This arguement against the presence of coherent state in 
Charlie’s viewpoint is only overcome if inter-related phase reference frames at 
the relevant photon frequencies actually exist. 

3.3 Reference Frames and Violations of Snperselection Rules 

Challenges to the requirement for quantum states to be consistent with super¬ 
selection rules have occured since the 1960’s when Aharonov and Susskind [53] 
suggested that coherent superpositions of different charge eigenstates could be 
created. It is argued that super-selection rules are not a fundamental require¬ 
ment of quantum theory, but the restrictions involved could be lifted if there is 
a suitable system that acts as a reference for the coherences involved - [59] , [60], 
|6T], [62], [53], [64], [53], [56], [57], [34] provide discussions regarding reference 
systems and SSR. 

3.3.1 Linking SSR and Reference Frames 

The discussion of the super-selection rule issue in terms of reference systems is 
quite complex and too lengthy to be covered in the body of this paper. However, 
in view of the wide use of the reference frame approach a full outline is presented 
in Appendix (TU] The key idea is that there are two observers - Alice and 
Charlie - who are describing the same prepared system in terms of their own 
reference frames and hence their descriptions involve two different quantum 
states. The reference systems are macroscopic systems in states where the 
behaviour is essentially classical, such as large magnets that can be used to 
define cartesian axes or BEC in Glauber coherent states that are introduced to 
define a phase reference. The relationship between the two reference systems 
is represented by a group of unitary transformation operators listed as T{g), 
where the particular transformation (translation or rotation of cartesian axes, 
phase change of phase references, ..) that changes Alice’s reference system into 
Charlie’s is denoted by g. Alice describes the quantum state via her density 
operator, whereas Charlie is the external observer whose specification of the 
same quantum state via his density operator is of most interest. There are 
two cases of importance. Situation A - where the relationship between Alice’s 
and Charlie’s reference frame is is known and specified by a single parameter g, 
and Situation B - where on the other hand the relationship between frames is 
completely unknown, all possible transformations g must be given equal weight. 
Situation A is not associated with SSR, whereas Situation B leads to SSR. The 
relationship between Alice’s and Charlie’s density operators is given in terms 
of the transformation operators (see Eq. (12251) for Situation A and Eq. (I226|) 
for Situation B). In Situation B there is often a qualitative change between 
Alice’s and Charlie’s description of the same quantum state, with pure states 
as described by Alice becoming mixed states when described by Charlie. It is 
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Situation B with the U(l) transformation group - for which number operators 
are the generators - that is of interest for the single or multi-mode systems 
involving identical bosons on which the present paper focuses. An example of the 
qualitative change of behaviour for the single mode case is that if it is assumed 
that Alice could prepare the system in a Glauber coherent pure state - which 
involves SSR breaking coherences between differing number states - then Charlie 
would describe the same state as a Poisson statistical mixture of number states - 
which is consistent with the operation of the SSR. Thus the SSR applies in terms 
of external observer Charlie’s description of the state. This is how the dispute 
on whether the state for single mode laser is a coherent state or a statistical 
mixture is resolved - the two descriptions apply to different observers - Alice 
and Charlie. On the other hand there are quantum states such as Fock states 
and Bell states which are described the same way by both Alice and Charlie, 
even in Situation B. The general justification of the SSR for Charlie’s density 
operator description of the quantum state in Situation B is derived in terms 
of the irreducible representations of the transformation group, there being no 
coherences between states associated with differing irreducible representations 
(see Eq. (12501) '). For the particular case of the U(l) transformation group the 
irreducible representations are associated with the total boson number for the 
system or sub-system, hence the SSR that prohibits coherences between states 
where this number differs. Finally, it is seen that if Alice describes a general non- 
entangled state of sub-systems - which being separable have their own reference 
frames - then Charlie will also describe the state as a non-entangled state and 
with the same probability for each product state (see Eqs. (12551) and (I256|) f. 
For systems involving identical bosons Charlie’s description of the sub-system 
density operators will only involve density operators that conform to the SSR. 
This is in accord with the key idea of the present paper. 

3.3.2 Can Coherent Superpositions of Atom and Molecules Occur ? 

Based around the reference frame approach Dowling et al |111] and Terra Cunha 
et al propose processes using a BEC as a reference system that would create 
a coherent superposition of an atom and a molecule, or a boson and a fermion 
m- Dunningham et al |112j consider a scheme for observing a superposi¬ 
tion of a one boson state and the vacuum state. Obviously if super-selection 
rules can be overcome in these instances, it might be possible to produce co¬ 
herent superpositions of Fock states with differing particle numbers such as 
Glauber coherent states, though states with N ^ 10® would presumably be 
difficult to produce. However, detailed considerations of such papers indicate 
that the states actually produced in terms of Charlie’s description are statistical 
mixtures consistent with the super-selection rules rather than coherent superpo¬ 
sitions, which are only present in Alice’s description of the state (see Appendix 
[ini). Also, although coherence and interference effects are demonstrated, these 
can also be accounted for without invoking the presence of coherent superpo¬ 
sitions that violate the super-selection rule. As the paper by Dowling et al 
m entitled "Observing a coherent superposition of an atom and a molecule.” 
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is a good example of where the super-selection rules are challenged, the key 
points are described in Appendix 1111 Essentially the process involves one atom 
A interacting with a BEC of different atoms B leading to the creation of one 
molecule AB, with the BEC being depleted by one B atom. There are three 
stages in the process, the first being with the interaction that turns separate 
atoms A and B into the molecule AB turned on at Feshbach resonance for a 
time t related to the interaction strength and the mean number of bosons in 
the BEC reference system, the second being free evolution at large Feshbach 
detuning A for a time r leading to a phase factor </> = A r, the third being 
again with the interaction turned on at Feshbach resonance for a further time 
t. However, it is pointed out in Appendix [TT] that Charlie’s description of the 
state produced for the atom plus molecule system is merely a statistical mixture 
of a state with one atom and no molecules and a state with no atom and one 
molecule, the mixture coefficients depending on the phase (j) imparted during 
the process. However a coherent superposition is seen in Alice’s description of 
the final state, though this is not surprising since a SSR violating initial state 
was assumed. The feature that in Charlie’s description of the final state no 
coherent superposition of an atom and a molecule is produced in the process is 
not really surprising, because of the averaging over phase differences in going 
from Alice’s reference frame to Charlie’s. It is the dependence on the phase (p 
imparted during the process that demonstates coherence (Ramsey interferome¬ 
try) effects, but it is shown in Appendix [TT] that exactly the same results can be 
obtained via a treatment in which states which are coherent superpositions of 
an atom and a molecules are never present, the initial BEC state being chosen 
as a Fock state. In terms of the description by an external observer (Charlie) 
the claim of violating the super-selection rule has not been demonstrated via 
this particular process. 

3.3.3 Detection of SSR Violating States 

Whether such super-selection rule violating states can be detected has also not 
been justified. For example, consider the state given by a superposition of a 
one boson state and the vacuum state (as discussed in [112) 1. We consider an 
interferometric process in which one mode A for a two mode BEC interferometer 
is initially in the state a |0) -I-11), and the other mode B is initially in the state 
|0) - thus |5'(i)) = (a |0)-|-/3 |l))yiC)|0)^ in the usual occupancy number notation, 
where |ap -I- |/3p = I. The modes are first coupled by a beam splitter, then a 
free evolution stage occurs for time r associated with a phase difference cp = Ar 
(where A = ws — uja is the mode frquency difference), the modes are then 
coupled again by the beam splitter and the probability of an atom being found 
in modes A, B finally being measured. The probabilities of finding one atom in 
modes A, B respectively are found to only depend on |/3p and (p. Details are 
given in Appendix 1111 There is no dependence on the relative phase between 
a and /3, as would be required if the superposition state a |0) -I- /3 |1) is to be 
specified. Exactly the same detection probabilities are obtained if the initial 
state is the mixed state p(i) = |Q;p(|0)^ (0|^ (g) |0)^ (Ol^) + |/3p(|l)^ (1|^ C) 
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|0)^ (0|g), in which the vacuum state for mode A occurs with a probability 
|ap and the one boson state for mode A occurs with a probability |/3p. In this 
example the proposed coherent superposition associated with the super-selection 
rule violating state would not be detected in this interferometric process, nor in 
the more elaborate scheme discussed in m- 

Of course, the claim that in isolated systems of massive particles it is not 
possible in non-relativistic quantum physics to create states that violate the 
particle number SSR - either for the sub-system states in a separable state or for 
any quantum state of the overall system - can be questioned. Ideally the claim 
should be tested by experiment, in particular when the number of particles is 
large in view of the interest in macroscopic entanglement since the Schrodinger 
cat was first described. The simplest situation would be to test whether states 
that violate the (local) particle number SSR could be created for a single mode 
system. Clearly, a specific proposal for an experiment in which the SSR could be 
violated is required, but to our knowledge no such proposal has been presented. 
Bose-Einstein condensates, in which all the bosons can occupy a single mode 
would seem an ideal candidate as a suitable bosonic system, and the Glauber 
coherent state is an example of a non SSR compliant state. For fermions, the 
Pauli exclusion principle would limit the number of fermions in a one mode 
system to be zero or one, but coherent superpositions of a zero and one fermion 
state are examples of non-SSR compliant pure states. As pointed out above, 
some authors such as [SU, [5S], [in], |5U] , base their definition of entanglement 
by allowing for the possible presence of non-SSR compliant sub-system states 
when defining separable states. The approach in Refs. [4] is based on the 
physical assumption that states that are non-compliant with particle number 
SSR - both local and global - do not come into the realm of non-relativistic 
quantum physics, in which the concept of entanglement is useful. Until clear 
evidence is presented that non-SSR- compliant states can be prepared, and in 
view of the theoretical reasons why they cannot^ it seems preferable to base the 
theory of entanglement on their absence when defining separable and entangled 
quantum states. 

3.4 Super-Selection Rule - Separate Sub-Systems 

In this sub-section the important case of SSR in separable states will be dealt 
with, since this is key to understanding what entangled states are allowed in 
systems involving identical particles. This forms the basis for the treatment of 
entanglement tests presented in the second part of this review (Paper II). 

3.4.1 Local Particle Number SSR 

We now consider the role of the super-selection rule for the case of non-entangled 
states. The global super-selection rule on total particle number has restricted 
the physical quantum state for a system of identical bosons to be of the form 
(llOOl) . Such states may or may not be entangled states of the modes involved. 
The question is - do similar restrictions involving the sub-system particle number 
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apply to the inodes, considered as separate sub-systems in the definition of non- 
entangled states ? The viewpoint in this paper is that this is so. Note that 
applying the SSR on the separate sub-system density operators ... is 

only in the context of non-entangled states. Such a SSR is referred to as a local 
SSR, as it applies to each of the separate sub-systems. Mathematically, the 
local particle number SSR can be expressed as 

[iVx,pf]=0 (123) 

where Nx is the number operator for sub-system X = A, B, ... .The SSR re- 
striction is based on the proposition that the density operators p^, p^, ... for 
the separate sub-systems A, B, ... should themselves represent possible quan¬ 
tum states for each of the sub-systems, considered as a separate system and 
thus be required to satisfy the super-selection rule that forbids quantum su¬ 
perpositions of Fock states with differing boson numbers. Note that if the local 
particle number SSR applies in each sub-system the global particle number SSR 
applies to any separable state. The proof is trivial and just requires showing 
that [N,Psep] = 0 

The justification of applying the local particle number super-selection rule 
to the density operators p^, p^, ... for the sub-system quantum states that 
occur in any separable state is simply that these are possible quantum states 
of the sub-systems when the latter are considered as separate quantum systems 
before being combined as in the Werner protocol [S] to form the separable state. 
Hence all the justifications based either on simple physical considerations or 
phase reference systems that were previously invoked for the density operator 
p of any general quantum states of the combined sub-systems to establish the 
global particle number super-selection rule apply equally well here. No more 
need be said. It is contended that expressions for the non-entangled quantum 
state p in which p^, p^, p^j... were not allowed quantum states for the sub¬ 
systems would only be of mathematical interest. 

Applying the local particle number SSR to the sub-system density operators 
for non-entangled states is discussed in papers by Bartlett et al HZ], m as 
one of several operational approaches for defining entangled states. As pointed 
out above other authors [33] , ISOj define separable (and hence entangled) states 
differently by specifically allowing sub-system density operators that are not 
consistent with the local particle number SSR, though the overall density op¬ 
erator is globally SSR consistent. The corresponding overall states are termed 
separable but non-local, and states that they would regard as separable would 
here be regarded as entangled. Examples of such states are given in Eqs. (|130l) 
and (I132|) . There are also other authors [54], [55] who define separable (and 
hence entangled) states via ([3]) but leave unspecified whether the sub-system 
density operators are consistent or inconsistent with the local particle number 
SSR. Note that any inequalities involving measured quantities that are found 
for separable states in which local SSR compliance is neglected must also ap¬ 
ply to separable states where it is required. The consequent implications for 
entanglement tests where local particle number SSR compliance is required is 
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discussed in the accompanying paper 11. Hence, in this paper we are advocating 
a revision to a widely held notion of entanglement in identical particle systems, 
the consequence being that the set of entangled states is now much larger. This 
is a key idea in this paper - not only should super-selection rules on particle 
numbers be applied to the the overall quantum state, entangled or not, but 
it also should be applied to the density operators that describe states of the 
modal sub-systems involved in the general definition of non-entangled states. 
The reasons for adopting this viewpoint have been discussed above - basically 
it is because in separable states the sub-system density operators must repre¬ 
sent possible quantum states for the sub-systems considered as isolate quantum 
systems, so the general reasons for applying the SSR will apply to these density 
operators also. Apart from the papers by Bartlett et al m , [53] we are not 
aware that this definition of non-entangled states has been invoked previously, 
indeed the opposite approach has been proposed m. [50]. However, the idea of 
considering whether sub-system states should satisfy the local particle number 
SSR has been presented in several papers - m], [5^, gZ], |S3], [SS], [57], [55], 
mainly in the context of pure states for bosonic systems, though in these papers 
the focus is on issues other than the definition of entanglement, such as quantum 
communication protocols [49] . multicopy distillation ET], mechanical work and 
accessible entanglement m, m and Bell inequality violation [55]. However, 
there are a number of papers that do not apply the SSR to the sub-system den¬ 
sity operators, and those that do have not studied the consequences for various 
entanglement tests. These tests are also discussed in the accompanying paper 
H. 

3.4.2 Criterion for Local and Global SSR in Separable States 

Theorem. A necessary and sufficient condition for all separable states for a 
given set of sub-system density operators to be global particle number 

SSR compliant is that all such sub-system states are local particle number SSR 
compliant. 

We first note that 

W^P] =^PRi{na,p‘}i]^PR +Pr ® (124) 

R 

Necessity: If the state p is globally SSR compliant then [N,'p] = 0. Taking 
the trace of both sides of (I124|) over sub-system space b, using Tr^ijo’^^ = 1 and 
Tr},{\nh,'p\f\) = 0 and then repeating the process for sub-system space a gives 

0 = Q = Y^P^{[n,X]) (125) 

R R 

which are operator equation in sub-system spaces a and b respectively. 

The Pr are not independent, satisfying Pr = 1. By choosing a particular 

R 
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(126) 


Ps we can write the last equation for sub-system a as 

0 = ^ + (1 - E PH){[na,rs]) 

R^S R^S 

where the remaining Pr are now independent. Differentiating the last equation 
with respect to Pr then gives 

0 = [na,PR] - [na,Ps] (127) 

for any two different R and S. Thus all the [na,Ps] must be the same. Using 
0 = II ?r]) again with equal [«„,?«] and ^ Pr = 1 we then see that 

R R 

all [na,'pR] must be zero. Similar considerations show that [nb,'pR] = 0. 

As these results apply for any choice of the Pr and of the (or we can 
then conclude that 

[«a,?R] = 0 [nb,PR] (128) 

which establishes that the sub-system states are local particle number SSR 
compliant. 

Note that the proof depended on the choice of the Pr being arbitrary apart 
from H-Pfl = 1- If the Pr are fixed then although we can show that 0 = 
R 

HT’fl([lia,?R]) = H-Pfl([lih,PR]), the steps leading to [uqjPr] = [«b,?R] = 0 
R R 

do not follow. The four sub-system states in Section [3.4.31 where all Pr = 1/4 
are not local particle number SSR compliant even though the overall state is 
global particle number SSR compliant. This would not be the case if any of the 
Pr differed from 1/4. 

Sufficiency: If the sub-system states are local particle number SSR compliant 
then [nQ,pR] = [nf,,pR] = 0. It then follows from p24|l that 

[fV,p]=0 (129) 

which establishes that the separable state is global particle number SSR com¬ 
pliant. This conclusion applies for arbitrary Pr- 

3.4.3 Global but not Local Particle Number SSR Compliant States 

However, it should be noted that some authors [49], [50] consider sub-system 
density operators in the context of two mode systems which comply with the 
global particle number SSR but not the local particle number SSR. Such a case 
involving four zero and one boson superpositions is presented by Verstraete et 
al [49], [50]. The overall density operator is a statistical mixture 

P = ^(IV'i)(V’il)AG|V'i)(V'il)B + ^(IV'i)(^il)AG|V'i)(V'*l)s 

+ ^(|V’-l) (V'-i|)a G \fj_i) (V'-i|)b + ® IV'-i) {V'-i|)s 

(130) 
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where \4’J) = (|0) + w\l))/\/2, with uj = The \'ip^) are superposi¬ 

tions of zero and one boson states and consequently the local particle number 
SSR is violated by each of the sub-system density operators IV'^) and 

\tp^) {'4’i^\)B- Although the expression in Eg. dldOl) is of the form in Eq.®, the 
subsystem density operators \ipa}) (V'cjDa and do not comply with 

the local particle number SSR, so in the present paper and in [3] the state 
would be regarded as entangled. However, Verstraete et al [49], [50] regard it as 
separable. They refer to such a state as separable but nonlocal. 

On the other hand, the global particle number SSR is obeyed since the 
density operator can also be written as 

p = i(|o)(o|Uo|o)(o|)B + i(|i)(i|Uo|i)(i|)B 

+ i(|vI/+)(vI/+|)^B (131) 

where 1^'+)^^ = (|0)^|1)5 -f |1)^ lO)^)/^^. This is a statistical mixture of 
N = 0,1,2 boson states. Note that Eg. (11311) indicates that the state could 
be prepared as a mixed state containing two terms that comply with the local 
particle number SSR in each of the sub-systems plus a term which is an entangled 
state of the two sub-systems. The presence of an entangled state in such an 
obvious preparation process challenges the description of the state as being 
separable. 

To further illustrate some of the points made about super-selection rules - 
local and global - it is useful to consider a second specific case also presented 
by Verstraete et al [49] . |50j . This mixture of two mode coherent states is 
represented by the two mode density operator 

P = I^\a,a}{a,a\ 

rif) 

= (I«)(«I).4®(I«)(«I)b (132) 

where \oi)Q is a one mode coherent state for mode C = A, B with a = |a| exp(—10), 
and modes A, B are associated with bosonic annihilation operators a, b. The 
magnitude |a| is fixed. 

This density operator appears to be that for a non-entangled state of modes 
A, B in the form 

'P = Y.Pr'Pr®'Pr (133) 

R 

with '^Pr ^ f ^ and —)> (|a) (q;|)^ and —>• (|a) (aDg. However al- 

R 

_ Q 

though this choice of p^ satisfy the Hermitiancy, unit trace, positivity fea¬ 
tures they do not conform to the requirement of satisfying the ( local) sub-system 
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boson number super-selection rule. From Eq. (11321) we have 

{n\{\a){a\)\m)^ = exp(-|an^i^ 

vn\ vm\ 

(p|(l«)(«l)k)B = (134) 

so clearly for each of the separate modes there are coherences between Fock 
states with differing boson occupation numbers. In the approach in the present 
paper the density operator in Eq. (11321) does not represent a non-entangled 
state. However, in the papers of Verstraete et al [H], [SD], Hillery et al [S3], 

[55] and others it would represent an allowable non-entangled (separable) state. 

Indeed, Verstraete et al [49] specifically state ”this state is obviously separable, 
though the states |a) are incompatible with the (local) super-selection rule.”. 

Verstraete et al [49] introduce the state defined in Eq. (11321) as an example of 
a state that is separable (in their terms) but which cannot be prepared locally, 
because it is incompatible with the local particle number super-selection rule. 

The mixture of two mode coherent states does of course satisfy the total or 
global boson number super-selection rule. The matrix elements between two 
mode Fock states are 

I in-i-m \fy\p+q ng 

{{n\A® {p\b)'p{\^)a®\<i)b) = exp(-2|ap) ^ exp(-ii(n -m+p- q)9) 

I In-l-m \fy\P+q 

= exp(-2|an ' ' ' ' (135) 

\/nWm\ VpWQ- 

These overall matrix elements are zero unless n-\-p = m + q, showing that there 
are no coherences between two mode Fock states where the total boson number 
differs. The mixture of two mode coherent states has the interesting feature of 
providing an example of a two mode state which satisfies the global but not the 
local super-selection rule. 

The reduced density operators for modes H, B are 

Pa = /^ (|a)(a|)A 'Pb = S^{W){a\)B 
and a straightforward calculation gives 

'Pa = exp(-|a|^)^^^^(|n)(n|)^ = exp(-|a|^) ^ (|p) (pD^ 

which are statistical mixtures of Fock states with the expected Poisson distribu¬ 
tion associated with coherent states. This shows that the reduced density oper¬ 
ators are consistent with the separate mode local super-selection rule, whereas 
the density operators 'p^ = (jo;) {oi\)a j = (|a) (aD^ are not . Later we will 
revisit this example in the context of entanglement tests. 

Note that if a twirling operation (see Eg. (12671) 1 were to be applied to mode A, 
the result would be equivalent to applying two independent twirling operations 
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to each mode. In this case the density operator for each mode is a Poisson 
statistical mixture of number states, so each mode has a density operator that 
complies with the local particle number SSR. 

3.4.4 Particle Entanglement Measure 

Wiseman et al have also treated entanglement for pure states [35] and mixed 
states m in identical particle systems, applying both the symmetrization 
principle and super-selection rules, invoking the arguement that without a phase 
reference the quantum state must be comply with the local (and global) particle 
number SSR. This is essentially the same approach as in El, El, 11 and in the 
present paper. For two mode systems the observable system density operator 
p is obtained from the density operator p that would apply if such a phase 
reference existed via the expression 

P= nn.4nfl = E (136) 

riAriB riAriB 


where Huatib = ^ projector onto sub-system states with ua, tib parti¬ 
cles in modes A, B respectively. Note that pn „^„3 is not nor¬ 

malised to unity. In fact the probability that there are tla, ns particles in modes 
A, B respectively is given by PnA-nB = = Tr(p(”^”®^), so 

Tr{'p) = PnAUB = 1- For separable states defined here as in Eq.(|3|), the 

riAriB _ 

expression in (11361) for the density operator is the same as that used here, since 
with p given by Eq.(|3|) and with (p^ (g) p|) Ubaub =Pr^Pr it is 

^ UAnB 

easy to show that p = Pj,gp. For general mixed states Wiseman et al introduce 
in Ref. m the idea of particle entanglement by defining its measure Epfp) 

Ep{T>)= E PriAriB = EpCp) (137) 

riAriB 

where i?M(p*'"^"®^) is a measure of the mode entanglement associated with the 
(unnormalised) state This might be taken as the entropy of mode 

entanglement Em{e) = —TrfaA InuA) for normalised density operators tr, 
where the reduced density operator for mode A is ua = FrB(u). Note that 
from Tln^nj^UmAmB = ^riAmA^riBniB ^uaub the particle entanglement measure 
Epifp) is the same for p, the observable density operator for the system. In the 
case of the separable state for modes Tl, B given in ([3|) it is straightforward to 
show that 


■AnAnB) 
r sep 


psep 

^riAriB 


Pr (htriA Pr ^riA ) ® (hlnB Pr ^ub ) 

R 

R 


(138) 

(139) 
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where n„^ and Xln^ are projectors onto sub-system states in modes A, B respec¬ 
tively with UA and ns particles in the respective modes = Bn a ® 

with PuAiPR) = TrAiBnAPn) and Pub^Pr) = being the proba¬ 

bilities of finding ua and ur particles in the respective modes when the corre¬ 
sponding sub-system states are and Since the state Psep is clearly a 
separable state of the form for the modes A, B, the corresponding measure 
of mode entanglement must be zero. It then follows from the general expression 
(|137l) that the particle entanglement measure is also zero for the separable state. 
This is as expected. 

Ep{p,,p)=0 (140) 

For the pure states considered in [39] we note that among them is the two 
boson state |1)^ 0 ll)^ which has one boson in each of the two modes A, 
B. The particle entanglement measure Epijj) is zero for this state (where 
p=(|l) (1|)^®(|1) (1 |)b), consistent with it being a separable rather than an 
entangled state. This indicates that Wiseman et al [35] do not consider that 
entanglement occurs due to symmetrization, as the first quantization form for 
the state might indicate. However, finding i?p(p) to be zero does not always 
shows that the state is separable, as the case of the relative phase state (defined 
in Appendix J of paper II, see also [32]) shows. As is shown there, Ep{'p) — 0 
for the relative phase state, yet the state is clearly an entangled one. Just as 
some entangled states have zero spin squeezing, some entangled states may be 
associated with a zero particle entanglement measure. Nevertheless a non-zero 
result for the particle entanglement measure i?p(p) shows that the state must 
be entangled - again we have a sufficiency test. However, as in the case of 
other entanglement measures the problem with using the particle entanglement 
measure to detect entangled states is that there is no obvious way to measure 
it experimentally. 

3.4.5 General Form of Non-Entangled States 

To summarise: basically the sub-systems are single modes that the identical 
bosons can occupy, the super-selection rule for identical bosons, massive or oth¬ 
erwise, prohibits states which are coherent superpositions of states with different 
numbers of bosons, and the only physically allowable Pr, Pr, ..for the separate 
mode sub-systems that are themselves compatiible with the local particle num¬ 
ber SSR are allowed. For single mode sub-systems these can be written as 
statistical mixtures of states with definite numbers of bosons in the form 

PR = J2^nAM {uaI PR=J2^nB\^B){nB\ ■■ (141) 

UA UB 

However, in cases where the sub-systems are pairs of modes the density 
operators pf, ..for the separate sub-systems are still required to conform 
to the symmetrisation principle and the super-selection rule. The forms for 
Pr, .. are now of course more complex, as entanglement within the pairs 
of modes Ai, A 2 associated with sub-system A, the pairs of modes Hi, B 2 
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associated with sub-system _B, etc is now possible within the definition for the 
general non-entangled state Eq. ([3]) for these pairs of modes. Within each pair 
of modes Ai, A 2 statistical mixtures of states with differing total numbers Na 
bosons in the two modes are possible and the sub-system density operators are 
based on states of the form given in Eq. (11021) . We have 






Pr 


Na 

Ca^{Na, k) 0 \Na — 

00 

EE P<i>NA 

Na=0 ■!> 


(142) 


with analogous expressions for the density operators etc for the other pairs 
of modes. Note that only involves quantum superpositions of states 

with the same total number of bosons Na. The expression (208) in Appendix 
B of paper 11 is of this form. 

3.5 Bipartite Systems 

We now consider the bipartite case where there are just two sub-systems in¬ 
volved. The simplest case is where each sub-system involves only a single mode, 
such as for two modes in a double well potential when only a single hyperfine 
state is involve. Another important case is where each sub-system contains two 
modes, such as in the double well case where modes with two different hyperfine 
states are involved. 

3.5.1 Two Single Modes - Coherence Terms 

The general non-entangled state for modes a and b is given by 

P = '^PrPr<S)Pr (143) 

R 


and as a consequence of the requirement that 'pr and 'pr are allowed quantum 
states for modes a and b satisying the super-selection rule, it follows that 

((«)")a - Tr{pi{ar) = 0 {(a^r)^=Tr{pi(a^r)=0 

^{br)^ = Trip^(br) = 0 ((6t)™)^ = rr(p|(&t)-) = 0 

(144) 

Thus coherence terms are zero. As we will see these results will limit spin 
squeezing to entangled states of modes a and b. Note that similar results also 
apply when non-entangled states for the original modes c and d are considered 
- ((c)”)c = etc.. 
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3.5.2 Two Pairs of Modes - Coherence Terms 

In this case the general non-entangled state where A and B are pairs of modes 
- Oi, 02 associated with sub-system A, and modes bi, 62 associated with sub¬ 
system B, the overall density operator is of the form (11431) . Consistent with 
the requirement that the sub-system density operators conform to the 

symmetrisation principle and the super-selection rule, these density operators 
will not in general represent separable states for their single mode sub-systems 
oi, 02 or &i, 62 - and may even be entangled states. As a result when considering 
non-entangled states for the sub-systems A and B we now have 

= T’r(p^(aJaj)”) 7^ 0 = 

= Tr{p^{blb,r)^0 i,j = l,2 (145) 

in general. In this case where the sub-systems are pairs of modes the spin 
squeezing entanglement tests as in Eqs.(50) - (52) in paper II for sub-systems 
consisting of single modes cannot be applied, as we will see. Nevertheless, there 
are still tests of bipartite entanglement involving spin operators. 
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4 Discussion and Summary of Key Results 

This paper is mainly concerned with two mode entanglement for systems of 
identical massive bosons, though multimode entanglement is also considered. 
These bosons may be atoms or molecules as in cold quantum gases. In the 
present paper we focus on the definition and general features of entanglement, 
whilst in the accompanying paper we consider spin squeezing and other tests 
for entanglement. 

The present paper starts with the general definition of entanglement for a 
system consisting of several sub-systems, and highlights the distinctive features 
of entangled states in regard to measurement probabilities for joint measure¬ 
ments on the sub-systems in contrast to the results for non-entangled or sepa¬ 
rable states. The relationship between entanglement and hidden variable theory 
is then explored followed by a discussion of key paradoxes such as EPR and 
violations of Bell inequalities. The notion of entanglement measures and entan¬ 
glement tests was briefly introduced, the latter being covered more fully in the 
accompanying Paper II. 

The paper then focuses on entanglement for systems of identical massive 
particles in the regime of non-relativistic quantum physics. A careful analysis 
is hrst given regarding the proper definition of a non-entangled state for sys¬ 
tems of identical particles, and hence by implication the proper definition of an 
entangled state. Noting that entanglement is meaningless until the subsystems 
being entangled are specihed, it is pointed out that whereas it is not possible 
to distinguish identical particles and hence the individual particles are not le¬ 
gitimate sub-systems, the same is not the case for the single particle states or 
modes, so the modes are then the the rightful sub-systems to be considered as 
being entangled or not. In this approach where the sub-systems are modes, 
situations where there are differing numbers of identical particles are treated as 
different quantum states, not as differing physical systems, and the symmetri- 
sation principle required of quantum states for identical particle systems will 
be satisfied by using Fock states to describe the states. 

Furthermore, it is argued that the overall quantum states should conform 
to the superselection rule that excludes quantum superposition states of the 
form (I99|) as allowed quantum states for systems of identical particles - massive 
or otherwise. Although the justication of the SSR in terms of observers and 
their reference frames formulated by other authors has also been presented for 
completeness, a number of fairly straightforward reasons were given for why it 
is appropriate to apply this superselection rule for massive bosons, which may 
be summarised as: 1. No way is known for creating such states; 2 No way 
is known for measuring all the properties of such states, even if they existed; 
and 3. There is no need to invoke the existence of such states in order to 
understand coherence and interference effects. Invoking the existence of states 
that as far as we know cannot be made or measured, and for which there are no 
known physical effects that require their presence seems a rather unnecessary 
feature to add to the non-relativistic quantum physics of many body systems, 
and considerations based on the general principle of simplicity (Occam’s razor) 
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would suggest not doing this until a clear physical justification for including 
them is found. As two mode fermionic systems are restricted to states with at 
most two fermions, the focus of the paper is then on bosonic systems, where 
large numbers of bosons can occupy two mode systems. 

However, although there is related work involving local particle number 
super-selection rules, this paper differs from a number of others by extending 
the super-selection rule to also apply to the density operators p^, ... for 
the mode sub-systems A, B, ... that occur in the definition ([3]) of a general 
non-entangled state for systems of identical particles. Hence it follows that the 
definition of entangled states will differ in this paper from that which would 
apply if density operators p^, p^, ... allowed for coherent superpositions of 
number states within each mode. In fact more states are regarded as entangled 
in terms of the definition in the present paper. Indeed, if further restrictions 
are placed on the sub-system density operators - such as requiring them to spec¬ 
ify a fixed number of bosons - the set of entangled states is further enlarged. 
The simple justification for our viewpoint on applying the local particle number 
super-selection rule has three aspects. Firstly, since experimental arrangements 
in which only one bosonic mode is involved can be created, the same reasons 
(see last paragraph) justify applying the super-selection rule to this mode sys¬ 
tem as applied for the system as a whole. Secondly, measurements can be 
carried out on the separate modes, and the joint probability for the outcomes 
of these measurements determined. For a non-entangled state the joint proba¬ 
bility (1231) for these measurements depends on all the density operators p^, p^, 
... for the mode sub-systems as well as the probability Pr for the product state 
... occuring when the general mixed non-entangled state is prepared, 
which can be accomplished by local preparations and classical communication. 
For the non-entangled state the form of the joint probability PAB..{i,j, •■•) for 
measurements on all the sub-systems is given by the products of the individual 
sub-system probabilities P^ii) = etc that measurements on the 

sub-systems A,B, ...yield the outcomes etc when the sub-systems are in 
states Pfi, Pfi, ... , the overall products being weighted by the probability Pr 
that a particular product state is prepared. If p^, p^, ... did not represent 
allowed quantum states then the interpretation of the joint probability as this 
statistical average would be unphysical Thirdly, attempts to allow the density 
operators Pr, Pr, ... for the mode sub-systems to violate the super-selection 
rule provided that the reduced density operators p^, pg for the separate modes 
are consistent with it are shown not to be possible in general. 

As well as the above justifications for applying the super-selection rule to 
both the overall multi-mode state for systems of identical particles and the 
separate sub-system states in the definition of non-entangled states, a more 
sophisticated justification based on considering SSR to be the consequence of 
describing the quantum state by an observer (Charlie) whose phase reference is 
unknown has also been presented in detail in Appendix [10] for completeness. For 
the sub-systems local reference frames are involved. The SSR is seen as a special 
case of a general SSR which forbids quantum states from exhibiting coherences 
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between states associated with irreducible representations of the transformation 
group that relates reference frames, and which may be the symmetry group for 
the system. 

In regard to entanglement measures, we discussed the particle entanglement 
measure of Wiseman et al [31], [H] and found that a non-zero result for the 
particle entanglement measure shows that the state must be entangled. How¬ 
ever, as for other entanglement measures the problem with using the particle 
entanglement measure to detect entangled states is that there is no obvious way 
to measure it experimentally. On the other hand, as will be seen in the accom¬ 
panying Paper II, the quantities involved in entanglement tests can be measured 
experimentally. 
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6 Appendix 1 - Projective Measurements and 
Conditional Probabilities 

6.0.3 Projective Measurements 

For simplicity, we will only consider projective (or von Neumann) measurements 
rather than more general measurements involving positive operator measure¬ 
ments (POM). If O is a physical quantity associated with the system, with 
eigenvalues \i and with 11^ the projector onto the subspace with eigenvalue Xi 
then the probability P{i) that measurement of O leads to the value Xi is given 
by EH 

P{i) = Tr(n,p) (146) 

For projective measurements Hi = 11^ = IlJ and = 1, together with 

i 

ofii = fiiO = Aifii. 

Following the measurement which leads to the value Xi the density operator 
is different and given by 

Pcond(^o) = (n.pno/P(*) (147) 

This is known as the reduction of the wave function, and can be viewed in 
two ways. From an ontological point of view a quantum projective measure¬ 
ment changes the quantum state significantly because the interaction with the 
measurement system is not just a small perturbation, as it can be in classical 
physics. From the epistemological point of view we know what value the physi¬ 
cal quantity 11 now has, so if measurement of H were to be repeated immediately 
it would be expected - with a probability of unity - that the value would be Xi. 
The new density operator Pcondi^p) satisfies this requirement. It also satisfies 
the standard requirements of Hermitiancy, unit trace, positivity - as is easily 
shown. 

To show this formally we have for the mean value for 11 following the mea¬ 
surement 

= TrljiXondi^P)) 

= Tr(fl (n,pn,))/P(i) 

= A,Tr(n,p)/P(i) 

= A, (148) 

whilst for the variance 

(aII^)^ = Tr((ll-(fl)/p_^(f],*)) 

= Tr(ll2p,,^(*))-(fl)' 

= >H->H 

= 0 (149) 
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which is zero as expected. 

If following the measurement of the results of the measurement were dis¬ 
carded then the density operator after the measurement is 

Pcondi^) = Pcondi^, *) = XI (150) 

i i 

which is the sum of the 'Pcondi^i 0 ®^ch weighted by the probability P{i) of the 
result Ai occuring. Note that the expression for 'Pcondi^) is not the same as the 
original density operator p. This is to be expected from both the epistimological 
and ontological points of view, since although we do not know what value Xi has 
occurred, it is known that a definite value for has been found, or that mea¬ 
surement process has destroyed any coherences that previously existed between 
different eigenstates of We note that 'Pcondi^) ^iso satisfies the standard 
requirements of Hermitiancy, unit trace, positivity - as is easily shown. 

6.0.4 Conditional Probabilities 

Suppose we follow the measurement of Q resulting in eigenvalue Xi with a mea¬ 
surement of A resulting in eigenvalue where the projector associated with 
the latter measurement is . Then the conditional probabiltity of measuring 
A resulting in eigenvalue following the measurement of that resulted in 
eigenvalue Xi would be 

P{j\i) = 

= Tr{E,i%p%))/Pit) 

= Tri{E,%)p{flMJ))/Pi^) (151) 

where the cyclic properties of the trace and the idempotent property of the 
projector have been used. If the measurements had taken place in the reverse 
order the conditional probabiltity of measuring 17 resulting in eigenvalue Xi 
following the measurement of A that resulted in eigenvalue p^ would be 

P{i\j) = Tr((n,S)p(g,nO)/P(j) (152) 

We note that the actual probability of measuring Xi then p^ would be the 
joint probability 

P{j after i) = P{j\i) P{i) = Tr{(E.jfii)p (fliEj)) (153) 

whilst the actual probability of measuring p^ then Xi would be the joint prob¬ 
ability 

P{t after j) = Pm PU) = Tr((n,S,) p (S.n,)) (154) 

and we note that in general these two joint probabilities are different. 

If however, the two physical quantities commute, then there are a complete 
set of simultaneous eigenvectors \ Xi,Pj) for 17 and A. It is then straightforward 
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to show that Hi'E.j = in which case P{j after i) = P{i after j) = P{i,j), 

so it does not matter which order the measurements are carried out. The overall 
result 


= P{M P{^) = pm P{j) 

= TrifiiEjpEjfii) 

= Tr{uMjp) (155) 

is an expression of Bayes theorem. 

A case of particular importance where this occurs is in situations involving 
two or more distinct sub-systems, in which the operators 17 and A are associated 
with different sub-sj^stems. For two sub-systems A and B the operators 17 and 
A are of the form 17^ and 17^, or more strictly 17^ (8) Is and 1a ® 17b. It is 
easy to see that (17^ ® 1b){1a ® 17b) = 17 0 17b = (1a 0 17b)(17a 0 1b), so 
the operators commute and results such as in Bayes theorem (I155p apply. 


6.0.5 Conditional Mean and Variance 

To determine the eonditioned mean value of A after measurement of 17 has 
led to the eigenvalue Xi we use rather than p in the mean formula 

^A^ = Tr(Ap). Hence 


= Tr(Ap^„„^(17,i)) 


= Tr(A (n,pn,))/P(i) 

(156) 

A = ^ 

3 

(157) 

^)i = 


3 

= ^^3 Tm3A{pXifEj)/P{i) 

j 



(158) 


j 


using Ej = E^, the cyclic trace properties and Eg. (11511) . Hence the conditional 
mean value is as expected, with the conditional probability P{j\i) replacing 
P{j) in the averaging process. 

For the conditioned variance of A after measurement of 17 has led to the 
eigenvalue Xi we use Pgo„(;(17,7) rather than 'p and the conditioned mean (A) 
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rather than 

in the variance formula = Tr{{A — (^Ap^'p). 

Hence 


(aA^)^ = rr((A-(A)/^_,(H,*)) 



= rr((A-(A)/(n,pno)/P(*) 

(159) 

Now 

(A-(A)/=y:(,,-(A)/H, 

j 

(160) 

so that 

j 




(161) 


j 


using the same steps as for the conditioned mean. Hence the conditional variance 
is as expected, with the conditional probability P{j\i) replacing P{j) in the 
averaging process. 

6.1 Detailed Inequalities for EPR Situation 

For separable states the conditional probability that measurement of pA on sub¬ 
system A leads to eigenvalue pa given that measurement of ps on sub-system 
B leads to eigenvalue pb is obtained from Eq. dsni) as 

P{pa-,Pa\pb,Pb) = ^PrPa{pa-,Pa)Pb{pb-,Pb)/^PrPb{pb,Pb) (162) 

R R 


where 


P^{pa,Pa) = TrAiU^ji) pB{pb,Pb) = TrsiU^J^) (163) 

are the probabilities for position measurements in the separate sub-systems. 
The probability that measurement of ps on sub-system B leads to eigenvalue 
Pb is 

P{pb,Pb) = ^Pr Pb(pb,Pb) (164) 

R 

The mean result for measurement of pA for this conditional measurement is 
from Eq. m 




^PaP{pa,Pa\pb,Pb) 

PA 

^Pr {pa)rPb{pb,Pb)/P{pb,Pb) 

R 


(165) 
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where 


(166) 


^Pa)r = '^PaPa{pa,Pa) 

PA 

is the mean result for measurement of pA when the sub-system is in state 

The conditional variance for measurement of pA for the conditional mea¬ 
surement of pb on sub-system B which led to eigenvalue ps is from Eq. (uni) 

= Y^iPA-{PA)p^^p^f P{pa.Pa\pb,Pb) 

PA 

= ^^P^a)’^b,pb^b{pb,Pb)IP{pb,Pb) ( 167 ) 

R 

where 

(^Pa)pb,pb = Y^Pa - ^Pa)pb,pb)'' Pa{pa.Pa) 

PA 

is a variance for measurement of pA for when the sub-system is in state but 
now with the fluctuation about the mean for measurements condi¬ 

tional on measuring ps ■ 

However, for each sub-system state R the quantity (Ap^)^ is minimised 
if {pa)pj^ is replaced by the unconditioned mean {pa)r just determined from 
Pfl. Thus we have an inequality 

{^p\)Ypp, > ( Api )^ ( 168 ) 

where 

{^Pa)^ = Y^P^ ~ Pa{pa,Pa) ( 169 ) 

PA 

is the normal variance for measurement of pA for when the sub-system is in 
state p^. 

Now if the measurements of ps are unrecorded then the conditioned variance 
is 




T.{^Pa)^p^^,^P{xb,xb) 

XB 

EY^R i^PA)sB,XB^B{xB,XB) 

XB n 


which in view of inequality (15(11) satisfies 

{^Pa)p^ > J^Y^b (Apa)^ Pb{pb,Pb) 

PB R 

= YPBi^fAf 

R 


( 170 ) 


( 171 ) 
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using ^ Pb{pb,Pb) = 1- Thus the variance for measurement of momentum pA 

PB 

conditioned on unrecorded measurements for momentum ps satisfies an inequal¬ 
ity that only depends on the variances for measurements of pA in the possible 
sub-system A states 
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7 Appendix 2 - Inequalities 

These inequalities are examples of Schwarz inequalities. 

7.1 Integral Inequality 

If C'(A), D{X) are real, positive functions of A and P(A) is another real, positive 
function then we can show that 

J dXP{\)C{\). J dXP{\)D{X) > (/ dXP{X)y^CiX)D{X)y (172) 

To show this write x = f dXP{X)C{X) and y = J dX P{X)D{X). Then 

xy = j dXP{X)C{X) j dyiP{y)D{yi) 

= U dXdyiP{X)P{y)C{X)D{y) 

= fdX P(XfC(X)D(X) + ffdXd^(l- S(X - y)) P{X)P{y)C{X)D{y) 

(173) 

Also, write z = dX P{X)y^C{X)D{X)^ . Then 

2 = / dX P(A) VC{X)DiX) J dy P{y) sJC{y)D{y) 

= JJdXdy PiX)P{y) VC{X)D{X) ^/C{y)D{y) 

= fdX P{XfC{X)D{X) +ffdXdy(l- 6{X - y)) P{X)P{y) y'C{X)D{X)^/C{y)D{y) 

(174) 

so that 

xy-z = UdXdy{l-5{X-y)) P{X)P{y) (c{X)D{y) - ^C{X)D{XWC{y)D{y)) 

= ]^UdXdy{l-5{X-y))P{X)P{y) (c{X)D{y) + C{y)D{X) - 2^C{X)D{y)^C{y)D{X)) 

= If fdXdy(l-S(X-y)) P(X)P(y) (VC(X)D(y) - y/C(y)D(X)) ' 

> 0 (175) 

which proves the result. 

For the special case where D(X) = 1 and where J dXP{X) = 1 we get the 
simpler result 

fdXP(X)C(X) > (/ dXP(X)^/C(X)J (176) 


7.2 Sum Inequality 

If Cr and Dr are real, positive quantities for various R and Pr is another real, 
positive quantity then we can show that 

Pr^rJ2 Pr Dr > {'Ej Pr'J CrDr\ ( 177 ) 

R R \ R J 
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To prove this write x = Y^ Pr Cr and y = Y1 Pr Pr Then 
R R 


xy = PrP-rJ2 PsDs 

R S 

= J2J2 PrPs CrDs 

R S 

= Pr CrDr + X] S (1 “ ^Rs) Pr Ps CrDs 

R R S 

Also, write z = Pr ^/CrDr^ . Then 

2 = (ePr VCrDr) (j:Ps VCsDs] 


(178) 


= J2J2 PrPs \JGrDr\JC'sDs 


R S 


= J2 Pr CrDr + X) Z) (1 “ ^Rs) Pr Ps \JCrDr\JCrDs (179) 

R R S 


SO that 
xy — z = 


Z) Z) -Pr -Ps (1 “ ^Rs) (CrDs - sjC rDrsJOSDs'] 

R S 

2 ZZ PrPs(1-<Jrs) (CrDs + CrDr - 2\JCrDssJCrDr^ 
^ZZ PsPr{1-5rs) [Vc^-y^C^y 


> 0 


(180) 


which proves the result. 

For the special case where Dr = 1 and where Z Pr = 1 we get the simpler 

R 

result 


EPrCr > (ZPrV^ 


(181) 


This inequality is used in [54]. 
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8 Appendix 3 - Spin EPR Paradox 

8.1 Local Spin Operators 

For two sub-systems 1 and 2 there are numerous possibilities for defining sep¬ 
arate commuting spin operators for the two systems. One situation of interest 
is where each sub-system is associated with two bosonic modes, the standard 
annihilation operators being oi and bi for system 1 and 02 and 62 for system 2. 
The local spin operators for each sub-system can be defined as 

Si = {b[ai + d\bi)/2 Si = {b\di - a\bi)/2i Sl = {V[bi - d\ai)/2 

si = (^02-I-0262)72 S'y = (6^02 - 0262)721 S'f = (6|62 - 0202)72 

(182) 


These satisfy the usual angular momentum commutation rules and those or the 
different sub-systems commute. The squares of the local vector spin c^erators 
are related to the total number operators Ni = b\bi -|-ojoi and N 2 = 'h\b 2 -1-0202 
as E(^a)" = ( 6 Vi72)(1Vi 72 + 1) and E(^a)" = ( N 2 / 2 ){N 2/2 + 1). The total 


spin operators are 

Sa = sl + sl a = x,y,z 

and these satisfy the usual angular momentum commutation rules. 


(183) 


8.2 Conditional Variances 

The question is whether the conditional variances (^ASh'^ ^ for measuring 

Sxi for sub-system 1 having measured Sx 2 for sub-system 2, and for 

measuring Syi for sub-system 1 having measured Sy2 for sub-system 2 violate 
the Heisenberg Uncertainty Principle. 

7 (*■)'“ (‘84) 

where the measurements on sub-system 2 are left unrecorded. If this inequality 
holds we have an EPR violation. 

For separable states the conditional probability that measurement of Sxi on 
sub-system 1 leads to eigenvalue Sxi given that measurement of Sx2 on sub¬ 
system 2 leads to eigenvalue Sx2 is obtained from Ea. IpO)) as 

P{Sxl, Sxl\Sx2: — E Pr P^{Sxl, Sxi)P^{Sx 2, 8x2 )/^ Pr P2^{Sx 2,8x2) 

R R 

(185) 

where 

Pi^{Sxi, 8 xi) = Tr.iUl^j,) P^{Sx 2 . 8 x 2 ) = Tr 2 if(l^Jil) (186) 
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are the probabilities for position measurements in the separate sub-systems. 
The probability that measurement of Sx 2 on sub-system 2 leads to eigenvalue 
Sx2 _ 

P^Sx2-,^x2^ — E PrP^{Sx2,Sx2 ) (187) 

R 

The mean result for measurement of Sxi for this conditional measurement 
is from Eg. dTOl) 

{SxlJ^ — / ^ Sxl P{Sxl, Sxl\Sx2^ ^x2^ 

' ' Sx 2 

^xl 

= Pi^iSx2,Sx2)/P{Sx2,Sx2) (188) 

R '' 'R 

where 

( 5 , 1 )^ = ^ SxiP^SxuSxi) (189) 

^xl 

is the mean result for measurement of Sxi when the sub-system is in state 
The conditional variance for measurement of Sxi for the conditional mea¬ 
surement of Sx 2 on sub-system 2 which led to eigenvalue Sx 2 is from Eq. o 


^ V(s.l-(^.l). fP(Sxl,Sxl\Sx 2 ,Sx 2 ) 

\ / Sx 2 i^x 2 ' ' Sx 2 i^x 2 


= P2{Sx2.Sx2)/P{Sx2,Sx2){I^Q) 

o'' Sx2,Sx2 
ri 


where 


^S'tiU =y^{Sxl-{Sxl)^ fP^(Sxl,Sxl) 

Sx2,Sx2 ' ' Sx2^Sx2 

^xl 


is a variance for measurement of Sxi for when the sub-system is in state 


but now with the fluctuation about the mean (Sxi)^ for measurements 

^ ' Sx2-,Sx2 

conditional on measuring Sx 2 - 

However, for each sub-system state R the quantity (AS'^^V is min- 

' ' Sx2iSx2 


imised if ( Sxi ) ^ is replaced by the unconditioned mean ( Sxi ) just de- 

' ' Sx 2 -,^x 2 ' ' ^ 

termined from Thus we have an inequality 


/ ^ \ R 

\ / Sx 2 .Sx 




where 


(ASl,y = ^(s,i - {SxiyPi^{Sxi,Sxi) 


(191) 


(192) 
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is the normal variance for measurement of Sxi for when the sub-system is in 
state 

Now if the measurements of Sx 2 are unrecorded - as would be the case from 
the point of view of the experimenter on spatially well-separated sub-system 1 
when measurements on this sub-system take place at the same time - then the 
conditioned variance is 



P{Sx 2 ,Sx 2 ) 

P^{Sx2.^.-^ (193) 

Sx2 o ' ' Sx2.Sx2 


which in view of inequality (15(11) satisfies 

Sx2 ^ 

= (194) 

R 

using P'^{Sx 2 : Sx 2 ) = 1- Thus the variance for measurement of spin Sxi 

Sx2 

conditioned on unrecorded measurements for spin Sx 2 satishes an inequality 
that only depends on the variances for measurements of Sxi in the possible 
sub-system 1 states 

Now exactly the same treatment can be carried out for the variance of spin 
Syi also conditioned on unrecorded measurements of measurements for momen¬ 
tum Sy2. Details are given in Appendix [5] We have with 




Sy2 ,Sy2 



YUsI,) PiSy2,Sy2) 

Sy2 ' ' Oy2,Psy2 

— (Syi\ ) P{Sy2: Sy2\Sy2: Sy2) 

\ / Oy2,Sy2 

Syl 

y ^ Syl P{Syl, Syl|S'y2 5 Sy2) 

Syl 


the inequality 

R 

with 

{AS^yl^ = Y^iSyl - {Syi)r Pl^{Syl,Syl) ( 196 ) 

Syl 


is the normal variance for measurement of Syi for when the sub-system is in 
state 
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We now multiply the two conditional variances, which it is important to 
note were associated with two different conditioned states based on two different 
measurements - spin Sx 2 and spin Sy 2 - carried out on sub-system 2. 

However, from the general inequality in Eq.([T77|) 

PrDr p (yI, Pr \/ CrDb\ (198) 

R R \ R J 

we then have 

i (e Ph\/(aS^.>''(a%)'‘) (199) 

But we know from the HUP that for any given state that ( ) — 

\ \ \Szi ) P, so for the conditioned variances associated with a separable state 

(Ag,).jAsi.)._^ > 1(ep»i(&>"i)' 

= ( 200 ) 


showing that for a separable state the conditioned variances ( AS. 


^xl 


and 




Sy2 


still satisfy the HUP. It is important to note that these variances 


were associated with two different conditioned states based on two different 
measurements - spin Sx 2 and spin Sy2 - carried out on sub-system 2, the results 
of which the observer for sub-system 1 would be unaware of. Thus if the EPR 
violations as dehned in Eg. (11841) are to occur then the state must be entangled. 
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9 Appendix 4 - Extracting Entanglement due to 
Symmetrisation 


9.1 Two Particle Case - Bosons 

The approach of Killoran et al [53] can be first applied to the simple case 
of = 2 bosons initially in the A modes aO and al and were discussed in 
Subsection 13.1.11 Here we present the detailed derivation of the results. The B 
modes bO and 51 are initially unoccupied. 

The occupied state is 

1$^) = ^{|aO(l)) |al(2)) + |a0(2)) |al(l))} (201) 

in first quantisation. This is regarded by Killoran et al [95] as an entangled state 
for sub-systems consisting of particle 1 and particle 2. In second quantisation 
the occupied state |<i>A) and the unoccupied state |$b) are given by 


|1>a) = |$s) = |0),o|0),i (202) 

These are regarded as separable states for the A modes aO and ol and separable 
states for the B modes 50 and 51. 

In second quantisation we consider the effect of the beam splitter on an input 
state 

|^zn) = |$A)®|$B) (203) 

The effect is to produce an output state given by 




(r5j -b tal) {rb\ + ta\) 




^/T 


|0).o |0)al ® |0),o |0)w 


^"(|0)aO I0)al ® Il)w Il)w) + ^^(I0)a0 |l)al ® |1)60 |0),1 + |l)a0 |0)al ® |0),o |l)w) 
+i'(|l)ao|l)al®|0)M|0)w) (204) 


Measurements can then be done on the output state based on projecting the 
state onto eigenstates for the number operators for the A and B mode-based 
sub-systems. The projectors n'^(iVA) for sub-system A onto eigenstates with 
Na = 0,1, 2 bosons are given by 

n^(o) = |o)„olo)„i(oLo(oLi 

n^(l) = (|l)ao| 0 )al(l|ao( 0 |al + |OU|l)al(OL(l|al) 

n^(2) = (|2)„„ |0)„, (2|„„ (OL, + |i),o |I),, (iLo (IL, + |0),o |2)^, (OLo (2^,) 

(205) 
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with corresponding expressions for projectors II^{Nb) for sub-system B. 

To demonstrate entanglement extraction for particle based sub-systems with 
particle 1 in one sub-system, and particle 2 in the other sub-system we choose 
projectors corresponding to there being one particle in the A modal sub-system 
and one particle being in the B modal sub-system. Thus the output state is 
projected onto the states with Na = 1 and Nb = 1 and we get after normalising 

|^o.t(i,i)) = AA(n^(i)®n^(i)) 

= ■^(|l)o0 l*^)ol ® I*^)b0 1^)61 ■*" I^)a0 l^)al ® |1)60 |0)&l) 

(206) 

This is still a bipartite entangled state of the of two modal sub-systems, A and 
B. 

If we construct a mathematical correspondence of the form 

|a0(l)) ^ |l),ol0)„i |al(2))^ |0),o|l),i 

|al(l)) ^ |0),o|l)ai |a0(2))^ |l),o|0),i (207) 

we see that the projected output state given in (12061) as a bipartite entangled 
state of the of two modal sub-systems, A and B, has the same mathematical form 
as the bipartite entangled state of the of two particle sub-systems containing 
particle 1 and particle 2. respectively. 

9.2 Two Particle Case - Fermions 

Here the details for the simple case of = 2 fermions initially in the C modes cO 
and cl are presented, following the same approach as in the previous SubSection. 
The D modes dO and dl are initially unoccupied. Fermion modes are denoted 
c and d to distinguish them from bosonic modes a and b. 

The occupied state is 

|ci>c) = Icl(2)) - |c0(2)) |cl(l))} (208) 

in first quantisation. This is regarded by Killoran et al |95] as an entangled state 
for sub-systems consisting of particle 1 and particle 2. In second quantisation 
the occupied state |<i)c) and the unoccupied state |<i)D) are given by 

I'J’c) = |l).o|l)cl \‘^D) = \0),o\0)dl 

I'J’c) = \^d) = \0),o\0),, (209) 

These are regarded as separable states for the C modes cO and cl and separable 
states for the D modes dO and dl. 
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In second quantisation we consider the effect of the beam splitter on an input 
state 

|$zn) = |$c)®|$D) (210) 

The effect is to produce an output state given by 

|$out) = U\^in) 

_ (rdj) + tCg) {rd[ + tcj) 


yr 


yr 


|0)eo|0).l®|0),o|0),i 


— ^^(|0)cO l®)cl ® |l)dO |l)dl) + |0)cO |l)cl ® |l)dO |0)dl + |l)cO l®)cl ® |0)do |l)dl) 


+i"(|l)coll)ci®|0)dolO)di) 


( 211 ) 


Note the minus sign in the second term - this is due to the fermion creation 
operators anti-commuting dgcj = —c(dg. 

Measurements can then be done on the output state based on projecting the 
state onto eigenstates for thejrumber operators for the C and D mode-based 
sub-systems. The projectors n'^(N"c) for sub-system C onto eigenstates with 
Nc = 0,1, 2 bosons are given by 

n''(o) = |0),o|0),i(0|,g(0|,, 

n^(i) = (|i)colo)ci(iU(oU + |o).oli).i(oU(iU) 

n^( 2 ) = (| 2 )^g |0),1 ( 2 |^g (0|,1 + |l),g |1),1 (l|,g (1|,1 + | 0 ),g | 2 )^, (0 | ,g ( 2 |^J 

( 212 ) 

with corresponding expressions for projectors Ii^{Nio) for sub-system D. 

To demonstrate entanglement extraction for particle based sub-systems with 
particle 1 in one sub-system, and particle 2 in the other sub-system we choose 
projectors corresponding to there being one particle in the C modal sub-system 
and one particle being in the D modal sub-system. Thus the output state is 
projected onto the states with Nc = 1 and No = 1 and we get after normalising 

|$o.t(i,i)) = A^(n^(i)®n^(i)) 

= ^ ~ l®)dl) 

(213) 


This is still a bipartite entangled state of the of two modal sub-systems, C and 
D. 

If we construct a mathematical correspondence of the form 


|c0(l)) ^ |l),o|0)el |cl(2))^|0),g|l),i 

|cl(l)) ^ |0),g|l),i |c0(2))^|l),g|0),i (214) 

we see that the projected output state given in (12131) as a bipartite entangled 
state of the of two modal sub-systems, C and H, has the same mathematical form 
as the bipartite entangled state of the of two particle sub-systems containing 
particle 1 and particle 2. respectively - even down to the correct minus sign in 
the second term. 
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9.3 Three Particle Case - Bosons 

The key ideas in the approach by Killoran et al [53] are more fully illustrated 
by considering one of their specific cases, namely a quantum state with iV = 3 
identical bosons for a system with four modes - two A modes aO and al and two 
B modes 50 and 5f - in which there are two bosons in mode aO and one boson in 
mode of. The other modes 50 and 5f are initially un-occupied. The modes aO 
and 50 could be two different spatial modes for a bosonic atom in one hyperfine 
state, and af and 51 could be two different spatial modes for a bosonic atom in 
another hyperfine state. With particles labelled 1, 2 and 3 the quantum state in 
terms of first quantisation is given by 

1$^) = ^{|a0(l))|a0(2))|al(3)) + |a0(l))|a0(3))|al(2)) 

+ |a0(2)) |a0(l)) |al(3)) + |a0(2)) |a0(3)) |al(l)) 

+ |a0(3)) |a0(l)) |al(2)) + |a0(3)) |a0(2)) |al(l))} (215) 

In first quantisation there is no state for the B modes, since the vacuum state 
is not recognised as a quantum state of anything. Following the approach of 
regarding labelled identical particles as sub-systems we consider a bipartite di¬ 
vision of the three particle system with the first sub-system as consisting of 
particle 1 and 2 and the second sub-system consisting of particle 3. The same 
state |$y 4 ) can be written as 

1$^) = ^{|a0(l))|a0(2))}|al(3)) 

+ y|^{|a0(l)) |al(2)) + |a0(2)) |al(l))} |a0(3)) (216) 

In this form the state appears to be an entangled state for the two sub-systems. 
The Hrst term (which has amplitude I/a/S) represents a state for the sub-system 
of particle 1 and 2 with both particles in single particle state |a0) and a state for 
the sub-system of particle 3 with this particles in single particle state |al). The 
second term (which has amplitude y^2/3) represents a state for the sub-system 
of particles 1 and 2 with one particles in single particle state |a0) and the other 
in single particle state |al), and a state for the sub-system of particle 3 with 
this particle in single particle state |a0). It is this entanglement which Killoran 
et al [95] wish to extract by applying a beam splitter to the state |$^), the beam 
splitter being associated with a unitary operator U whose effect is to transform 
the single particle states \ak) and \bk) into linear combinations of each other 
involving reflection and transmission coefficients r, t as follows 

U\ak) = r\bk)+t\ak) fc = 0,1 

U\bk) = t\bk)-r\ak) k = 0,l (217) 

For simplicity r, t are assumed to be real with r'^ +t'^ = 1. The beam splitter is 
just assumed to couple spatial modes of the same hyperfine state. This unitary 
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operator applies irrespective of the particular particle occupying the one particle 
states. 

In second quantisation the occupied state |$^) and the unoccupied state 


3 ) are given by 




I^a) = 

|2)a0 l^)al 

= |0),o |0),i 


\^ a ) = 

(aj)^ (al) 

V2 Vi 

|0)ao|0)al l<i>s) = |0),o|0)bl 

(218) 


where the Fock states are also written in terms of mode creation operators and 
vacuum states for the modes. In second quantisation it is clear that |$a) and 
Id’s) themselves are respectively separable states for the A and B modes. In 
second quantisation the effect of the unitary operator associated with the beam 
splitter follows from (I217|) noting that \ak) = |0) and \hk) = 6 ). |0) and is 

given by 

UalU-^ =rVf^ + tal ral fc = 0,1 (219) 

In paper II we show that two mode beam splitters can indeed be described by 
equations analogous to (12191) . In second quantisation we consider the effect of 
the beam splitter on an input state 


l^m) = \<^a) 0 |$b) 


( 220 ) 


The effect is to produce an output state given by 




(rfej +talf {rV^ 


ta\) 


V2 


^/T 


|0)aO I0)al ® \0)W |0) 


/bl 


’’ (I*^)a0 l*^)al ® |2)b0 |l)&l) + ^ ^(I®)o0 l^)ol ® I^)b0 ’^(|l)a0 l*^)al ® |1)&0 |l)bl) 

+rt^maO |0)al ® |0)bO |l)bl + ^(|l)a0 |l)al ® |l)bO |0)bl) + ^"(12)^0 |l)al ® |0)bO |0)bl) 

( 221 ) 


Note this state is normalised to unity as = (r^+t^)^ = 1. The input 

state is a bipartite separable state of two modal sub-systems, one containing the 
two A modes aO and al and the other the two B modes 60 and 61. The output 
state terms each are eigenstates of number operators Na = ajao + a\ai and 
Nb = 6960 -I- 6 | 6 i with eigenvalues Na = 0, Nb = 3 for the term, Na = 1, 
Nb = 2 for the r'^t term, Na = 2, Nb = 1 for the rt^ term, Na = 3, Nb = 0 
for the t^ term. The same result as in (12211) can also be obtained using (12171) in 
conjunction with the first quantisation form of the input state given by (12151) 
though the algebra is more complex. In contrast to the input state, the output 
state is a bipartite entangled state of two modal sub-systems, one containing 
the two A modes aO and al and the other the two B modes 60 and 61. Both 
input and output states are states with the same total of = 3 bosons. 


no 














Measurements can then be done on the output state based on projecting the 
state onto eigenstates for the number operators for the A and B mode-based 
sub-systems. The projectors 11'^(iV^) for sub-system A onto eigenstates with 
Na = 0,1, 2, 3 bosons are given by 


n^(o) = |o)„olo),i(oLo(oLi 

n^(l) = (|l)ao|0)al(l|ao(0|al + |OU|l)al(OL(lLl) 

n^(2) = (|2)„o |0)al (2Lo (OLi + |l),o |l)al (iLo (iLl + |0).o |2)al (OLo (2Ll) 

n^(3) = (|3)„o |0)„, (3^0 (OL, + |2)^o |1)^, (2^0 (IL, + |l)^o |2)^, (2^, + |0),o |3),, (3Li) 

( 222 ) 


with corresponding expressions for projectors 11^ (TVs) for sub-system B. 

To demonstrate entanglement extraction for particle based sub-systems with 
particles 1 and 2 in one sub-system, and particle 3 in the other sub-system we 
choose projectors corresponding to there being two particles in the A modal sub¬ 
system and one particle being on the B modal sub-system. Thus the output 
state is projected onto the states with Na = 2 and Nb = 1 and we get after 
normalising 


l$o.t(2,l)) 


Ar(n^(2)®n^(i)) 

|2)a0 l*^)ol ® |0)bO |l)hl + Y^(|l)a0 |l)ol ® |l)b0 |0)bl 

(223) 


This is still a bipartite entangled state of the of two modal sub-systems, A and 
B. 

If we construct a mathematical correspondence of the form 

|a0(l)) |a0(2)) ^ |2),o |0)„i |al(3)) ^ |0),Jl),i 

^{|a0(l))|al(2)) + |a0(2))|al(l))} ^ |l),o |l)ai |a0(3)) ^ |1),J0),, 

(224) 

we see that the projected output state given in (12231) as a bipartite entangled 
state of the of two modal sub-systems, A and B, has the same mathematical form 
as the bipartite entangled state of the of two particle sub-systems containing 
particles 1 and 2 and particle 3. respectively. This type of result is proved in 
more general cases in [SS] - here we have exhibited the key features of their 
approach in a particular case. 

It is on this basis that Killoran et al [53] assert that the action of the beam 
splitter is to extract the entanglement due to symmetrisation that was present in 
the quantum state |$y 4 ) for the particle sub-systems containing particles 1 and 
2 and particle 3. respectively. It is of course not their ingeneous mathematical 
derivation that is in dispute - it is the interpretation. From the point of view 
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of sub-systems being modes, not particles the action of the beam splitter is 
to create an entangled state of the two modal sub-systems, A and B from a 
state that was separable. That this entangled output state can be projected 
onto eigenstates of the number operators for the two modal sub-systems, A and 
B which have the same mathematical form as the presumed entangled initial 
state for the particle sub-systems containing particles 1 and 2 and particle 3. 
respectively is of course interesting, but it does not show that labeled identical 
particles can be regarded as physically accessible sub-systems. Apart from the 
logical issue that sub-systems must be both distinguishable from each other via 
physical measurements, it is noteworthy that the approach of Killoran et al 
[55] rested on physical processes that involved coupling modes, not identified 
indistinguishable particles. 
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10 Appendix 5 - Reference Frames and Super- 
Selection Rules 

Several papers such as [Sg, ED, [S3], |SS], [33], [S7], [55] explain the link between 
reference frames and super-selection rules (SSR). In this Appendix we present 
the key ideas involved. 

10.1 Two Observers with Different Reference Frames 

The first point to appreciate is that there are two observers - Alice and Char¬ 
lie - who are involved in describing the same quantum system, which has been 
prepared via some physical process We will refer to Charlie as the external ob¬ 
server, Alice the internal observer. The system could be a multi-mode system 
involving identical particles, it could just be a single mode system or it could 
even be a single particle with or without spin. Alice and Charlie each describe 
quantum states in terms of their own reference frames, which might be a set of 
coordinate axes for the case of the spin or position states for the single particle 
system, or it could be a large quantum system with a well-defined reference phase 
in the case of multi-mode or single mode systems involving identical particles. 
Alice and Charlie may each choose from a set of possible reference frames - for 
the single particle case there are an infinite number of difference choices of coor¬ 
dinate axes for example, related to each other via rotations and/or translations. 
In Situation A - which is not associated with SSR - Alice and Charlie do know 
the relationship between their two reference frames (and can communicate this 
relationship via classical communications) - such as in the case of the single 
particle system when the relative orientation of their two different coordinate 
axes are known. In Situation B - which is associated with SSR - Alice and 
Charlie do not know the relationship between their two reference frames - such 
as in the multi-mode or single mode system involving identical particles when 
the relative phase between their two large quantum phase reference systems is 
not known. Alice and Charlie describe the same system via density operators 
a and p, and the key question is the relationship between these two operators 
in situations A and B and for various types of reference frames. In terms of the 
notation in m p ^ a and p —>■ p. In some situations the assumption that Alice 
even possesses a well-defined reference frame may be invalid, in which case it is 
important to realise that it is Charlie’s quantum state which is of most interest 
for describing the system. This description may differ from what hypothetical 
observer Alice would regard as the quantum state. 

10.2 Symmetry Groups 

A particular relationship going from Alice’s to Charlie’s reference frame is speci¬ 
fied by the parameter g, which in turn defines a unitary transformation operator 
T{g) that acts in the system space. Particular examples will be listed below. 
If there was a third observer - Donald - and the relationship going from Char¬ 
lie’s to Donald’s reference frame is specified by the parameter h, which in turn 
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defines a unitary operator T{h), then if we symbolise the relationship going 
from Alice’s to Donald’s reference frame by the parameter hg, it follows that 
T{hg) = T{h)T{g). This shows that the unitary operators satisfy one of the 
requirements to constitute a group, referred to generally as the transformation 
group. The other requirements are easily confirmed. The unitary operator 
T(0) = 1 corresponding to the case where no change of reference frame oc¬ 
curs (specified by the parameter 0) exists, and satisfies the requirement that 
T{0g) = T{0)T{g) = T{g0) = T{g)T{0). The unitary operator T{g~^) = T{g)‘^ 
corresponding to the relationship specified as that converts Charlie’s ref¬ 
erence frame back to that of Alice exists, and satisfies the requirement that 
T(0) = T{g~^)T{g) = T{g)T{g~^). Hence all the group properties are satisfied. 

A few examples are as follows: 

1. Translation group - single spinless particle system, with p, x.the momen¬ 
tum, position vector operators. Here is a vector giving the translation of 
Charlie’s cartesian axes reference frame from that of Alice, thus g =.j^. The 

unitary translation operator is T{j^) = exp(ip • a^/h). 

2. Rotation group - single particle system, with J the angular momentum 

vector operators. Here is a unit vector giving the axis and rotation angle 
(f) for rotating Alice’s cartesian axes reference frame into that of Charlie, thus 
g =.j^., The unitary rotation operator is (f) = exp(i(/)J • 

3. Particle number U(l) group - single mode bosonic system, with a the 
mode annihilation operator and Na = the mode number operator. Here 9a 
is the phase change Alice’s to Charlie’s reference frame. The unitary operator 
is f{9a) = e^p{iNa9a). 

4. Particle number U(l) group - multi-mode bosonic system, with a as a 

typical mode annihilation operator and N = of a the total number operator. 

a 

Here 9 is the phase change from Alice’s to Charlie’s reference frame. The unitary 
operator is T{9) = exp(iiV0). 

In these examples the system operators p, J, Na, N etc are the generators 
of the respective groups. In many situations the generators commute with the 
Hamiltonian for the system (or more generally with the evolution operator that 
describes time evolution of the quantum state), in which case the group of 
unitary operators T{g) is the symmetry group, and the generators are conserved 
physical quantities. 

10.3 Relationships - Situation A 

In Situation A, where the relationship between the reference frames for Alice 
and Charlie is known and specified by a single parameter g, Alice’s description 
of the state a is related to Charlie’s description p for the same state via the 
unitary transformation 

p = fig)afig)-^ (225) 
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Note that this is a passive transformation - no change of state is involved, just 
the same state being described by two different observers. 

As an example, consider the spinless particle and the translation group. If 
is a position eigenstate then T(_^) = A pure quantum po¬ 

sition eigenstate described by Alice as ct = 1^1’) (‘hI with state vector |$) = 

would be described by Charlie as p = |d>) (dt| but now with 
which is also a pure quantum position eigenstate but with eigenvalue — a^. 
This is as expected since Alices’s cartesian axes have been translated by _c^ to 
the origin of Charlie’s axes without change of orientation. In the case of mo¬ 
mentum eigenstates 


p^j we have T{a^) p^j = exp(ij^ • j^/h) so a pure 

quantum momentum eigenstate described by Alice with | $) = p ^ would be 

described by Charlie with |d>) = exp(jp • j^/h p which is also a pure mo¬ 
mentum eigenstate with the same eigenvalue p . Alice and Charlie describe the 

pure momentum eigernstate with the same density operator 'p = a, the phase 
factor cancels. 

For more general pure states, consider a quantum state described by Alice 
as a = 1$) ($1 with state vector |$) = j States of this form 

can represent localised states when is only significant in confined spatial 

regions, or they can represent delocalised states, such as momentum eigenstates 
^ when (^(_^) = exp(i^ • We see that Charlie also describes 

a pure quantum state p = |dt) (^'1 but now with |dt) = T(_c^) |$) = / d^((){^ + 

4 ) 4 ) = / ^4 ^>( 4 ) I 4 ), so the wavefunction is now + ^). 

Note that if Alice’s state vector was written in terms of momentum eigen¬ 
states |$)=/dp4p) pV then Charlie’s state vector l^-) = [ dp ^{p) 


has a momentum wave function = exp(z^ • j^/h) 4>{p^ related to that 

of Alice by a phase factor. Note that a state which is a quantum superposition 
of momentum eigenstates as described by Alice is also described as a quantum 
superposition of momentum eigenstates by Charlie. A similar feature applies in 
all situation A cases, and is related to SSR not applying in situation A. 

The case of the particle with spin and the rotation group is outlined in Ref. 


10.4 Relationships - Situation B 

In Situation B, where on the other hand the relationship between frames is 
completely unknown, all possible transformations g must be given equal weight, 
and hence the relationship between Alice’s and Charlie’s description of the same 
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state becomes 


p = J w{g)dgf{g)af{g) ^ 

= g [a] (226) 

where J w{g)dg is a. symbolic integral over the parameter g, which includes a 

weight factor w{g) so that J w{g)dg = 1. This linear process connecting a to'p 

is the ”t/- twirling'’’ operation. Again, this is a paassive transformation. 

It is straightforward to show that for any fixed parameter h that 

f{h)pf{h)-^=p (227) 

showing that Charlie’s density operator is Q invariant under the transformation 
group - unlike the case for Situation A. 

As an example, consider the single mode bosonic system and the U(l) group. 
If \na) is a Fock state then T(9a) \na) = exp(mo0a) \na)- Consider a pure quan¬ 
tum state described by Alice as the Glauber coherent state a = |$) ($| with state 
vector |$(/3)) = ^ C{na,(i) |na), where (^(na,/3) = exp(-|/3p/2)/3”“ /^J{na)\. 

Ua 

It is straightforward to show that 

f(0,) |$(/3)) = |$(/3exp(*0,))) (228) 

so that the Glauber coherent state is transformed into another Glauber coherent 
state, but with /3 changed via a phase factor to /3exp(i0a)- The quantum state 
described by Charlie is given by 

P = y ^ l$(/3exp(i6»o))) (^’(/3exp(i6»a))| (229) 

= C{na, f3) C{ma, py f{da) \na) {ma\ f{6ay 

Ua TTla 

= EE C{na,P)C{ma,Py \na) {rrial j ^ exp(I[na - ma]0a) 

Ua rria 

= \C{na,P)\‘^ \na) {Ual 

Ua 

= E V \na){na\ (230) 

which is a mixed state consisting of a Poisson distribution of Fock states with 
mean occupation number ria = |/3p- In view of the first expression for 'p it 
can also be thought of as a mixed state consisting of Glauber coherent states 
each with the same amplitude |/3| = y/W^, but with all phases (arg/3 -|- 0a) 
equally probable. Thus, whereas Alice describes the state as a pure state that 
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is a quantum superposition of Fock states with differing occupancy numbers, 
Charlie describes the same state as a mixed state involving a statistical mixture 
of number states. The former violates the SSR whereas the latter does not. A 
similar feature applies in all situation B cases, and is related to SSR applying 
in Situation B. Whether Alice could ever prepare such a state in the first place 
is controversial - see the discussion presented above in SubSections 13.21 and 13.41 
However, assuming she could, the quantum state as described by Charlie is a 
mixed state. 

The situation just studied relates of course to the debate [102] regarding 
whether the quantum state for a single mode laser operating well above thresh¬ 
old should be described by a Glauber coherent state or as a Poisson statistical 
mixture of photon number states. The first viewpoint (Alice) describes the 
state from the point of view of an internal observer with a reference frame, 
the second (Charlie) describes the same state from the point of view of an ex¬ 
ternal observer for whose reference frame relationship to that of the internal 
observer is unknown. The debate is regarded by [53] as settled on the basis that 
both viewpoints are valid, they are just at cross purposes because they refer to 
descriptions of the same quantum state by two different observers. 

It should not be thought however that the quantum state would always be de¬ 
scribed in such a fundamentally different manner for all Situation B cases. As an 
example, consider the multi-mode bosonic system and the U(l) group. Consider 
the pure quantum state described by Alice as the multi-mode N boson Fock state 
a = 1$) ($1 with state vector |$(Af)) = \nin 2 ---na...\N) — n \ni) Inz) .. \na) ..., 

a 

where N = ^ Uq. We have T{d) |nin 2 ...na...; N) = exp(iA^0) |nin 2 ...na...; N), 

a 

so that the same state would be described by Charlie as p = (dtj and with 
l^t) = \nin 2 ...na---; N). This is also a multi-mode N boson Fock state with 
exactly the same occupancies. The product exp{iN9) exp{—iN9) of phase fac¬ 
tors averages out to unity and here p = ct, so Alice and Charlie both describe 
the multi-mode Fock states in the same way. Another example for two mode 
bosonic systems and the U(l) group is provided by the one boson Bell states 
(the BS notation used here is non-conventional). These are entangled two mode 
states that Alice would describe via the state vectors |$^) = (|10) ± |01))/-\/2. 
We have T{9) = exp(i0) |$=*=), so that the same state would be described 

by Charlie with |4'=*=) = (|10) ± |01))/-\/2. Again the product of phase factors 
averages to unity and p = tr, so Alice and Charlie both describe the quantum 
states as Bell states, and in the same form. 

10.5 Dynamical and Measnrement Considerations 

Discussions of the relationship between equations governing the dynamical be¬ 
haviour of Alice’s and Charlie’s density operators depend on whether the evolu¬ 
tion is just governed by a Hamiltonian or whether master equations describing 
evolution affected by interactions with an external environment are involved. 
Such matters will not be treated in detail here, nor will the issue of relating 
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Alice’s and Charlie’s measurements. The latter issue is dealt with in |53) . 

However, in the case where Alice describes the Hamiltonian evolution of her 
density operator via the Liouville - von-Neumann equation 

d 

i/i—n = [i?, ct] (231) 

where in Alice’s frame the Hamiltonian is 77, and where in addition the trans¬ 
formation group is also the symmetry group so that T{g)HT{g)~^ = H for all 
g, it is easy to see that for both Situations A and B, Charlie’s density operator 
will evolve via the same LVN equation 

zh^^p=[H,p] (232) 

Thus both Alice and Charlie will describe the same dynamical evolution, though 
of course the initial (and hence evolved) states may differ in the two cases. 

10.6 Nature of Reference Frames 

Reference frames of differing types are involved for the various transformation 
groups. The common feature is that they are thought of as actual physical sys¬ 
tems themselves which are either macroscopic classical systems or macroscopic 
quantum systems in states associated with the classical limit. They are intended 
to be essentially unaffected by the presence of the systems for which they are 
acting as reference frames. In some cases relatively uncontroversial examples 
exist, such as for the cartesian axes associated with the translation and rotation 
groups associated with the single particle system. The physical reference system 
may be a large magnet whose magnetic field points in a well defined direction 
and defines a z axis, combined with an electrostatic generator whose electric 
field is in another well defined direction at right angles that defines an x axis. 
In other cases the existence of suitable reference frames is less clear. 

In this Subsection we will describe possible phase reference frames as if they 
are entirely separated (or uncorrelated) with the system of interest. In terms 
of the treatment by Bartlett et al [64j . [53) these are non-implicated reference 
frames. In the next SubSection and in the next Appendix phase reference frames 
that are correlated with the system of interest will be described - these are the 
so-called implicated reference frames of Bartlett et al. 

For the large quantum system with a well-defined reference phase associated 
with the U(l) group in the case of multi-mode or single mode systems involving 
identical particles, the usual choice is a single mode bosonic system such as a 
single mode BEC or a laser with a large mean occupancy, and which is thought 
of as being prepared in a Glauber coherent state |4)(q;)) in order to provide 
the phase reference frame, the reference phase being arg a. Whether such a 
reference frame really exists is controversial. The discussion presented above in 
Subsections 13.21 and 13.41 raises the question of whether such a phase reference 
state could ever be prepared, so this choice of a physical phase reference is rather 
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unsatisfactory. However, from the point of view of this presentation we assume 
it does, so that - as in the previous example - Alice can describe the reference 
state as another coherent state. Again, whether Alice could ever prepare such 
a state is questionable. 

Another possibility for a physical phase reference is a macroscopic low fre¬ 
quency harmonic oscillator, whose quantum energy eigenstates |n) - with n = 
0,1, ..., riniax and energies nhuj can be used to construct phase eigenstates \6p) 
with p = 0,1, ..., Umax and dp = p x 27r/(umax + 1), and which are defined by 

? 

^ ^max 

|0p) = exp(m6»p) \n) (233) 

These states are orthonormal. The separation between the equally spaced phase 
angles A0 = 27r/(niaax + 1) can be made very small if Umax is large enough. 
Under the effect of the harmonic oscillator Hamiltonian H = fuoN, where N is 
the number operator, the phase state \9p) evolves into \9p — ujAt) during a time 
interval At, so if the time intervals are chosen so that wAt = 27r/(nmax +1), the 
phase angle 9p changes into 0p_i. Thus the system behaves like a backwards 
running clock [?], the phase angles 9p defining the positions of the hands. If the 
clock initially has phase 9p the probability of finding the clock to have phase 9q 
after a time interval At is given by 


P{9q,9p,At) 


1 sin^((nniax + 1)A/2) 
(?T-max + l)^ sin^(A/2) 


(234) 


where A = 9p — 9q — wAt. For times At such that ojAt <C 27r/(nniax + 1) the 
probability of the phase remaining as 9p is close to unity. Thus if the phase state 
\9p) is used as a phase reference, it will remain stable for a time At satisfying 
the last inequality. For At ^ lOO^s and Umax 10^ so that phase is defined to 
^ 10“^ radians, an oscillator frequency w ^ 10° s“^ would sufhce for this phase 
reference standard. Such macroscopic oscillators do exist, though the process 
to prepare them in the phase reference quantum state \9p) would be technically 
difficult. Whether such a system would be useful as a phase reference for optical 
fields or a BEG is another issue 


10.7 Relational Description of Phase References 

In this Subsection phase reference frames that are correlated with the system 
of interest will be described - these are the so-called implicated reference frames 
of Bartlett et al [51] , [S3] . 

One such approach to describing phase references in the U(l) group case 
is via the concept of maps. For simplicity consider a one mode system S, the 
basis vectors for which are Fock states \m)g, where it is sufficient to restrict 
m = 0,1, ... ,mjniLx. The reference system R, will also be a one mode system 
with Fock states |n)^, where n is large. Product states \m)g 0 |n)^ for the 
combined modes exist in the Hilbert space Hs O Hfi and are eigenstates of the 
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various number operators, including the total number operator Nt = Ns + 

- where the eigenvalue is I = m + n. The product states may be listed via 
m = 0,1, ..., runiax and n = 0,1, ... or m = 0,1, ..., mmax and / = m, m + 1, .... 
Here we will describe how a coherent superpostion of number states, such as a 
Glauber coherent state can be represented. 

In the so-called internalisation or quantisation of the reference frame the 
product state Ito)^ G |n)^ is mapped onto the product state |m)g. G |n — m)^ 
where n > mmax- Thus 

\m) g. G |n)^ —>■ \m) g. G \n — (235) 

Hence for a linear combination of system states given by 

^max 

|$)g. = Cm |TO)g (236) 

m—0 

we have for the state |$)g 0 |n)^ in Hs G Hr 

ITlmax '^max 

Cm\m)s(^\n-m)i^=\'Sr,)Rs 

m—0 m—0 

(237) 

The mapping results in an entangled state where there are n bosons distributed 
betweeen the two modes. This state |'I'n)i{5 ^ state which is compatible 

with the SSR and is in one-one correspondence with the original system state 
|<I))g. Note that to create this state the reference state |n)^ must have more 
bosons in it than mmax- The density operator for the original pure system 
S state would be as = |d>)g (d>|g, and we note that this state violates the 
SSR. The state would be essentially a Glauber coherent state if Cm = 
exp(—|Q;p/2)a’"/(-\/^), with Wmax ^ |ap. However, for the mapped state 
|'I'„)^g. the reduced density operator ps is given by 

^max 

= \Cm\ \i^)si''^\s (238) 

m—0 

This is a mixed state and is compatible with the SSR. For the Glauber coherent 
state |$)g. this is the Poisson distribution of number states. Hence the original 
SSR violating superposition of number states for system S is mapped onto 
a state in the combined system for which the reduced density operator is a 
statistical mixture and is consistent with the SSR. as would correspond to 
Alice’s description of the state, 'ps to Gharlie’s. 

In the alternative so-called externalisation of the reference frame the map¬ 
ping is between product states, and is the reverse of the previous mapping. The 
product state \m)g G \n)R is mapped onto the product state \m)g \m + n) 
in the Hilbert space Hs G Hr where the former is spanned by vectors \m)g and 
the latter by vectors |m -I- n)^, and where n > mmax- Thus 

\m)g ^ \n)j^ ^ \m)g ^ \m + n)g (239) 


120 



The mapping of the Hs 0 Hji state then is 

^max 

m—0 

^max / ^max \ 

—> Cm, |w)g. 0 |n-)^ = I Cm |w)g I 0 = l“n)i{S 

(240) 

The mapping results in a non-entangled state which is incompatible with the 
SSR. The state in the subspace Hs is a coherent superposition of number states, 
whilst that in Hji is a Fock state. The reduced density operator in Hs is af 
given by 

^max ^max 

= E E \^)s {k\s (241) 

m—0 k—0 

which is the same as as = |4))g, ($|g and involves coherences between different 
number states in contradiction to the SSR. Clearly this second mapping just 
reverses the first one. 

Of these two treatments of phase reference frames, the internalisation ver¬ 
sion has a closer link to physics in that the pure state can in prin¬ 

ciple be created and does lead to a way of creating a state that is in one- 
one correspondence with any SSR violating pure state |d>)g, though it is in 
the form of an entangled state of the S, R sub-systems rather than just S 
alone. This is an important point to note - the original SSR violating state 
does not exist as a state of a separate system, all that exists is an SSR com¬ 
patible entangled state that is in one-one correspondence with it. However, 
the general process for creating a state such as is not explained. For 

simple cases such as |<i))g = (|0)g -I- \ l)g)/^/2 the creation of the required state 
= (l^)s ® \'’^)r + |l)s C) |n — l)^)/-\/2, where n > 1 would seem feasible 
via the ejection of one boson from a BEC in a Fock state \n)j^ into a previously 
unoccupied mode. . 

10.8 Irreducible Matrix Representations and Super-selection 
Rnles 

If \i) {i = 1,2, ...) are a set of orthonormal basis vectors in the system state 
space, then the group of unitary operators T{g) is represented by a group of 
unitary matrices D(g) 

ng)\i) = Y,Dj.{g)\j) (242) 

3 

with elements Dji{g), and such that D{hg) = D{h)D{g) etc corresponding to 
the group properties of the operators. This is a matrix representation of the 
transformation group. 
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The theory of such group representations and their application to quan¬ 
tum systems is well established, following the pioneering work of Wigner in the 
1930s. We can just use the results here. A key concept is that of irreducible 
representations. Within the system state space we can in general choose so- 
called irreducible sub-spaces, denoted as Tq, of dimension da and spanned by 
new orthonormal basis vectors iTaA) (A = 1,2, ... ,da) such that 

da 

f(g) |r„A) = (g) |r„/i) (243) 

fi=i 

For each irreducible sub-space Fq, there is no smaller sub-space for which the 
operation of all T{g) just leads to linear combinations of vectors within that 
sub-space. The da x da matrices D°‘{g) then form an irreducible matrix repre¬ 
sentation for the transformation group. For different a the representations are 
said to be inequivalent. 

The irreducible matrices satisfy the so-called great orthogonality theorem 

ma 

J w{g)dgD’^^{g)Df^{g)* = ^SapSf,^5xr (244) 

The proof of this result is based on Schur’s lemma. 

The importance of the irreducible representations and the consequent or¬ 
thogonality theorem lies in its application to Situation B cases, where we have 
seen that Charlie’s density operator p is invariant under any of the transforma¬ 
tions T{h)'pT{h)~^ = p. Suppose we represent p in terms of the basis vectors 
iFcA) associated with the irreducible representations 

p = |r„A)(rpr| (245) 

aA /3t 

where R will be a Hermitian, positive definite matrix with unit trace since it 
represents a density operator. Applying the transformation gives 

f{h)pf{h)-^ = Y.Y.^"rD;x{h)\rati){rpi\Dl^ihr 

= p (246) 

Averaging over h and using the great orthogonality theorem gives 

p = E E (E J ir^p) (r„p| (247) 

This is in the form of a mixed state involving irreducible state vectors iF^p) 
each occuring with a probability P“ given by 

= E ;r^AA = (248) 

A 
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which is the same for all /i associated with a given irreducible representation 
htj. This is clearly a positive real quantity and since 


ESC 

Oi fl 


a fi X ^ a X 

Tr'p = 1 


(249) 


the probabilities sum to unity as required. 

The final result for Charlie’s density operator 

p = EE |rc/r) (r«^| (250) 

a 

demonstrates the presence of a super-selection rule. In Charlie’s description 
of the quantum state there are no coherences between states iTa/r) associated 
with differing irreducible representations of the transformation group. This 
represents the general form of the SSR for all transformation groups in Situation 
B cases. 

As an example, consider the U(l) group and the single mode bosonic system. 
Since the Fock states satisfy T{0a) \na) = exp(ma0a) |«a) they form the basis 
for the irreducible representations of the U(l) group, the occupation number 
Uq specifying the irreducible representation and the 1x1 matrices exp(ma6*a) 
being the unitary matrices. Hence Charlie will describe the quantum state as 

d ='^P{'na) \na) {na\ (251) 

ria 


which is a statistical mixture of Fock states with no coherences between different 
Fock states. This result is of the same form as in Eg. (11411) and is in accord with 
the SSR on boson number. 

As another example, consider the U(l) group and the multi-mode bosonic 
system. Here sums of products of Fock states 

\n 1 n 2 -na-\N) ='^\ni)\n2) ■■\na) - N = '^na (252) 


such that the total occupancy is A^ = ^ Ua can be used to form irreducible 

a 

representations for the transformation group in terms of linear combinations of 
the products with the same N. Writing these linear combinations as 

E (253) 

{nin2...na---} 

we have since T{9) \nin 2 ...na...; N) = exp{iN0) \nin 2 ...na...', N) we see that 
T{6) = exp{iN9) |'I')(,) also, so the \'^%) define the irreducible basis states. 

The total occupancy N specifies the irreducible representation, but here there 
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are many irreducible representations with the same N depending on the various 
/i. In this case Charlie will describe the state as 

P = ( 254 ) 

N M 


which is a statistical mixture of multi-mode states all with the same total 
occupancy N. Although there are coherence terms between individual modal 
Fock states, there are no coherences between states with different total occu¬ 
pancy. This result is of the same form as in Ea. dlOOp and again is an example of 
a super-selection rule operating in terms of Charlie’s description of the quantum 
state. 

Finally, we note that in situation A where the relationship between the 
frames is known and there is no invariance for Charlie’s density operator, we do 
not have SSR applying. For the single particle case and the translation group 


the momentum states 


PJ 


define the irreducible representations, each specified 


by p , and as we saw Charlie’s description of the quantum state involved linear 


combinations of these irreducible basis vectors, in contradiction to the SSR. 


10.9 Non-Entangled States 

The essential feature of an non-entangled or separable state is that the sub¬ 
systems are considered to be unrelated to each other. Hence, both for Alice 
and Charlie there will be separate reference frames for each sub-system, with 
transformation groups - TA{ga) for sub-system A, Tsigb) for sub-system B, etc 
which relate the reference systems of Alice to those of Charlie. The transforma¬ 
tions Pa, Pb, ■■ are different. The overall transformation operator would be of 
the form T{ga,gb, ■■■) = TA{ga) <8 TB(gb) <8 ... • Alice would describe a general 
non-entangled state as having a density operator 

o’= ^ Pr (g) 0 ... (255) 

R 


It then follows for Situation B where the reference frames for Alice and Char¬ 
lie are unrelated, that Charlie would describe the same state via the density 
operator 

P = ^PrPb‘S>Pr<E)Pr<E)... (256) 

R 


where 


ffR = j w{gf)dgMgc)dlfc{gc)-^ C = A,B,... (257) 

Note that separate twirl operations are applied to the different sub-systems, as 
explicitly shown in the papers by Vaccaro et al [56] (see Section IIIA, Eqn. 3.3 
therein) and Paterek et al [58| (see Section 6). This leads for general trans- 
formation groups to the local group super-selection rule, where the Pr involve 
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no coherences between states associated with differing irreducible representa¬ 
tions of the transformation group. We see that Charlie also describes a non- 
entangled state and with the same mixture probability Pr as for Alice. Thus 
non-entanglement or separability is a feature that is the same for both Alice 
and Charlie, as ought to be the case. 

In the context of sub-systems consisting of modes (or sets of modes) occupied 
by identical bosons, the case of interest is Situation B, with each transformation 
group being U(l). Here the relationship between Charlie’s and Alice’s phase 
reference frames are unknown. Hence irrespective of Alice’s description of the 
sub-system states CTr, ... we see from the previous section that Charlie will 
describe the separate sub-system states Pr, Pr, as statistical mixtures of number 
states for the separate modes (or total number states for the sets of modes in 
each sub-system). Thus from Charlie’s point of view the separate mode density 
operators will satisfy the SSR. Thus we see that the introduction of reference 
frames and two observers - Charlie being the external one whose description of 
the quantum states is of primary interest - leads to the same SSR outcome as 
the simpler considerations set out in SubSections 13.21 and 13.41 Essentially the 
same considerations have been used in m, lis] and the other papers to justify 
the local photon number superselection rule. 
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11 Appendix 6 - Super-Selection Rule Violations 

? 

11.1 Preparation of Coherent Snperposition of an Atom 
and a Molecule ? 

A key paper dealing with the coherent superposition of an atom and a molecule 
is that by Dowling et al m, entitled “Observing a coherent superposition of an 
atom and a molecule”. Essentially the process involves one atom A interacting 
with a BEG of different atoms B leading to the creation of one molecule AB, 
with the BEG being depleted by one B atom. 


11.1.1 Hamiltonian 

The Hamiltonian is given by 


H = hlOAb^A^A + + ^(PnJ)Ab2 + 


(258) 


where bA,bM and 62 are standard bosonic annihilation operators for the atom, 
molecule and BEG modes respectively, uja,(^m and W 2 are the corresponding 
mode frequencies and k defines the interaction strength for the process where a 
molecule is created or destroyed from/to an atom A and a BEG atom B. A is 
the frequency difference between the molecular state AB and the two separate 
states for atoms A and B - this is zero on Eeshbach resonance - and is given by 


A = UJm — UJA — UJ2 


(259) 


The Hamiltonian commutes with the total number operator Ntot-, where 

Vot = ‘^'b'^M^M + b\bA +^62 (260) 

where the molecule number operator is multipled by two. 


11.1.2 Initial State 

Initially the state of the system is given by the density operator Eqs (10) and 
( 11 ) in the paper 

WoL = y ^exp(-iiVtot6»)|^')o^(^'|o^exp(+iiVtot6») (261) 

= IAI/3) (262) 

where |A) is a state with one atom A and |/3) is a Glauber coherent state for 
the BEG of atoms B.The super-operator acting on the pure state |'l')oi 
called the twirling operator, the group of unitary operators exp{—iNtotQ) depend 
on a phase variable 9 and are a unitary representation of D(l), the generator 
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being Ntot- These operators act as a symmetry group for the system and leave 
the Hamiltonian invariant. The initial state is also given by 


WoL 

— PA-m{^) ® P2{^) 

(263) 

Pa-m{^) 

= l^)MI 

(264) 

?2(0) 

f dO 

= / ^ exp(-m 26 ') (8) (/3| exp(+in20) 

(265) 


= '^Pn{<n>)\n) {n\ 

n 

(266) 


= y ^ l/3exp(-i6»)) (/3exp(-'i6»)| 

(267) 


where n 2 = 6262 is the number operator for the BEC mode and Pn{< n >) = 
{exp(— < n >) < n >" /n!} is a Poisson distribution, whose mean is < n > = 
|/3p. Initially then there is one atom A and the BEC is in a statistical mixture of 
number states with a Poisson distribution, which is mathematically equivalent 
to a statistical mixture of Glauber coherent states |/3exp(—i0)) with the same 
amplitude \/< n > but with all phases (arg /3 + 0) being equally weighted. 

11.1.3 Implicated Reference Frame 

In the paper by Dowling et al m the BEC is acting as an implicated phase 
reference frame (see [64] . m)- The state of the reference frame as described 
by Charlie is given by 


Pref = P2(0) = y ^ exp(-m 26 ») |/3) (/?| exp(+m 26 ») (268) 

and from Eq. (I258|) , there is an interaction between the reference BEC and the 
separate atom A and molecule M systems. However, because < n >= |/3p is 
very large, the BEC is essentially unchanged during the process, as reflected in 
the use of approximations in eqs (27), (28) of the paper. Another implicated 
phase reference frame situation, but involving a two mode reference frame is 
discussed in the paper by Paterek et al 

Overall, in terms of the discussion in Appendix ITOl Wnr would be Charlie's 
description of the initial state, whereas Alice would describe it as |'I')p^ ('^'Iol- 
Presumably in the paper by Dowling et al m what is referred to as the 
’’state of the laboratory” be Charlie’s reference frame, and what they refer to as 
the ’’internal reference frame” would refer to that of Alice. However, whether 
Alice could actually prepare such a state as i® controversial - see 

Subsections 13.21 and 13.41 though here this is assumed to be possible. 

11.1.4 Process - Alice and Charlie Descriptions 

There are three stages in the process, the first being with the interaction that 
turns separate atoms A and B into the molecule AB turned on at Feshbach 
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resonance for a time t = tt/(2k < n >), the second being free evolution at large 
Feshbach detuning A for a time r leading to a phase factor (j) = At, the third 
being again with the interaction turned on at Feshbach resonance for a further 
time t = Trj{2K < n >). The typical initial state |^')q£, given by |A) |/3) (eq 
(11)) evolves into |'I') 3 ^ given by (see eq. (32) of paper) 



(269) 


using approximations set out in eqs (27), (28) of the paper that depend on 
< n > being large. Here \M) is a state with one molecule AB. Thus it looks 
like a coherent superposition of an atom state |A) and a molecule state \M) 
has been prepared, the atom plus molecule system being disentangled from the 
BEC. Alice would describe the final state of the system as |^') 3 j;^ from 

her point of view a coherent superposition of an atom and a molecule has been 
prepared. 

However, for Charlie the final state of the system is described by a den¬ 
sity operator which is reconstructed by applying the twirling operator to 
l^) 3 L (^IsL • Noting that 



Tr2{\l3exp{—i6)) {fiexp{—i6)\) = (/3exp(—z0)|/3exp(—id)) = 1 (271) 

we see that Charlie’s final reduced density operator for the atom-molecule system 
is 


PA-Aii^) = Tr2W3L 



sin2(^) \A){A\+cos^{^) \M) (M| 


(272) 


Thus the coherence terms like |A) {M\ and \M) (A| do not appear in the final 
density operator when the average over 9 (not /3) is carried out. 

For Charlie the density operator for the atom and molecule is of course a 
statistical mixture of a state with one atom and no molecule and a state with 
no atom and one molecule. The authors of |111) actually point this out in the 
paragraph after eq (35) where (presumably for the case f = 7r/4) it is stated “the 
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state is found to be ... an incoherent mixture of an atom and a molecule.”. The 
probabilities for detecting an atom A or a molecule AB are as in eq (33) of the 
paper. In terms of Charlie’s description, the density operator at the end of the 
preparation process does not signify the existence of a coherent superposition of 
an atom and a molecule, as the title to the paper might be taken to imply. The 
existence of such a coherent superposition would of course be present in Alice’s 
description, but it is Charlie’s (laboratory) description that is more relevant. 

11.1.5 Interference Effects Without SSR Violation 

Note that interference effects are still present since the atom or molecule de¬ 
tection probabilities depend on the phase (j) associated with the free evolution 
stage of the process. However, as in many other instances, the presence of co¬ 
herence effects does not require the existence of coherent superposition states 
that violate the super-selection rule. The authors actually point this out in the 
paragraph after eq (35), where it is stated “we have clearly predicted the stan¬ 
dard operational signature of coherence, namely Ramsey type fringes, but the 
coherence is not present in our mathematical description of the system.” What 
they are referring to is Charlie’s description of the final state - which indeed 
shows no such coherence, but the belief that coherent superposition states are 
needed to predict coherence effects is mistaken. 

To drive this point home, the process can be treated with the initial state for 
the BEC being given as a Fock state \N). With the interaction being given as in 
Ea. (|258|) (eq (14) in the paper) the state vector is a simple linear combination 
of two terms 

|^(t)) = A{t) |A) \N) + Bit) \M) |iV - 1) (273) 

This is of course an entangled state. Coupled equations for the two amplitudes 
Af) and Bf) can easily be obtained and simple solutions obtained for stages 
where the Feshbach detuning is either zero or large. The state vector is continu¬ 
ous from one stage to the next , and the reduced density operator at the end of 
the three stage process for the atom plus molecule sub-system can be obtained. 
It is of the form 


Pa-m{3) = Tr2(|'k(3))(vl/(3)|) 

= sin^it)\A){A\+cos^i^)\M){M\ (274) 

which is of course a statistical mixture of a state with one atom and no molecule 
and a state with no atom and one molecule - and is exactly the same result as 
obtained in the paper by Dowling et al. m- Note that coherence effects in 
regard to the interferometric dependence on f for measurements on the final 
state has been found without invoking either the description of the BEC via 
Glauber coherent states or the presence of a coherent superposition of an atomic 
and a molecular state. The result can easily be extended for the case where the 
BEC is initially in a statistical mixture of Fock states with differing N occuring 
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with a probability Pn- Each initial state |A) |A^) evolves as in Eq. (12731) . We 
then would have 


N 

= EPn \A){A\+cos^i^) |M) (M|) 

= sin\^)\A){A\+cos^i^)\M){M\ (275) 

which is the same as before. Allowing for a statistical mixture of Fock states 
makes no difference to the interferometric result. 

11.1.6 Conclusion 

Dowling et al |lllj state in their abstract that “we demonstrate that it is pos¬ 
sible to perform a Ramsey-type interference experiment to exhibit a coherent 
superposition of a single atom and a diatomic molecule” . However the inter¬ 
ferometric effects (involving the dependence on cj)) cannot be said to exhibit the 
existence of such a coherent superposition, since the same interferometric re¬ 
sults can be obtained without ever introducing such a quantum state. There is 
not a convincing case that quantum states that violate the super-selection rule 
forbidding the creation of coherent superpositions of Fock states with differing 
particle numbers can be created, even in Alice’s reference system. The fact that 
an SSR violating state created in Alice’s reference system is not 

surprising, because in the process considered the initial state |/3) for the BEC 
was assumed as a factor in Alice’s initial state, and this was itself inconsistent 
with the SSR. Furthermore, such SSR violating states are not needed to de¬ 
scribe coherence and interference effects, so that justification for their physical 
existence also fails. 

11.2 Detection of Coherent Snperposition of a Vacuum 
and a One-Boson State ? 

Whether such super-selection rule violating states can be detected has also not 
been justified. For example, consider the state given by a superposition of 
a one boson state and the vacuum state (as discussed in [112j 'l. Consider an 
interferometric process in which one mode A for a two mode BEC interferometer 
is initially in the state a |0) -I-/3 |1), and the other mode B is initially in the state 
|0)-thus |d>(i)) = (a |0)-|-/3 |1 ))a®|0 )^ in the usual occupancy number notation, 
where |ap-|-|/3p = 1. Modes A, B could refer to two different hyperhne states of 
a bosonic atom with non-relativistic energies /iw^.and hujB, mode annihilation 
operators a, b. The modes are first coupled by a beam splitter, which could be 
a resonant microwave pulse that causes transitions between the two hyperhne 
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states and which can be described via a unitary operator Ubs such that 

® |0)b) = (|1)a ® |0)_B “ * I®)a ® I^)_b)/''^ 
%s(|0)^®|i)b) = H|i)a®|0)b + |0)^®|i)b)/^ 
^bs(|0)^®|0)5) = (|0)^®|0)^). (276) 

After passing through the beam splitter the system is allowed to evolve freely 
for a time r, the Hamiltonian being Hfree = {me? +hu}A)o?a-\- {mc^ + fiuJB)b^b - 
where collisional effects have been ignored and the rest mass energy included for 
completeness. Following the free evolution stage, the modes are then coupled 
again via a beam splitter, and the probability of an atom being found in modes 
A, B then being measured. A straightforward treatment of the evolution shows 
that the final state is given by 

!'!'(/)) = «(|0)^®|0)s) 

+/3 exp(—i{mc ^/h + oja}t) 

X 810),)-,i±i52^(|0),8ID,)) 

(277) 

where A = wb — to a is the detuning. The probabilities of finding one atom in 
modes A, B respectively are 

Pio = W sin2(AT/2) Poi = W cos2(At/2) (278) 

Thus whilst coherence effects occur depending on the phase difference (j) = At 
associated with the interferometric process, the overall detection probabilities 
only depend on the initial state via |/3p. There is no dependence on the rela¬ 
tive phase between a and /3, as would be required if the superposition state 
a|0) + /3|1) is to be specified from the measurement results. Exactly the 
same detection probabilities are obtained if the initial state is the mixed state 
p{i) = |ap(|0)^ (0|^® |0)s (0 |b) + |/3p(|l)^ (1^® |0 )b (0|b): in which the vac¬ 
uum state for mode A occurs with a probability |ap and the one boson state 
for mode A occurs with a probability |/3p. In this example the coherent super¬ 
position associated with the super-selection rule violating state would not be 
detected in the interferometric process. The paper by Dunningham et al m 
considers first a detection process that involves using a Glauber coherent state 
as one of the input states. Similar interference effects as in Eq. (12781) are ob¬ 
tained. A second detection process in which the single term Glauber coherent 
state is replaced by a statistical mixture with all phases equally weighted in 
considered next, leading to the same interference effects. This again confirms 
that it is not necessary to invoke the existence of coherent superpositions of 
number states in order to demonstate interference effects. 
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